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Untersuchungen ttber die Luftwiderstande 
von kérnigen und staubférmigen Giitern im Luftstrom* 


Von W. Rausch 


1. Einleitung; Zweck und Ziel der Untersuchung. Fiir das Studium der Strémungsvorgange 
von gasférmigen Medien mit Festteilchen, die heute auf vielen Gebieten der Technik, insbesondere 
in der Verfahrenstechnik, von groBer Bedeutung sind, ist die Kenntnis des Luftwiderstandes der 
Festteilchen von gréBter Bedeutung. 

Bei der Bewegung von Festteilchen in Gasen spielt neben der Masse der Teilchen vor allem der 
Luftwiderstand, der von der Strémung auf die Teilchen ausgeiibt wird, eine entscheidende Rolle. 
Der Luftwiderstand der Teilchen ist von der Form, der Gré8e und der Zahl der Teilchen in einem 
bestimmten Raum sowie einer Reihe anderer GréBen abhangig. Wenn man den Luftwiderstand 
der Teilchen und die Strémungsverhaltnisse kennt, kann man die Bahnen der Teilchen mit Hilfe 
von Bewegungsgleichungen' ** errechnen. Die in der Literatur angegebenen Bewegungsgleichun- 
gen sind jedoch nur beschrankt anwendbar und gelten nur unter der Voraussetzung, da die Gut- 
beladung so gering ist, daB die Strémungsverhaltnisse durch die Festteilchen nicht beeinfluBt 
werden. 

Bei derartigen Untersuchungen wurde also bisher in der Regel angenommen, daB die einzelnen 
Teilchen sich gegenseitig nicht beeinflussen. Diese Annahme trifft in den meisten Fallen nicht zu. 
Es ist eine altbekannte Tatsache, die schon von vielen Forschern beobachtet worden ist, daB ein 
Gemisch aus vielen Teilchen in einem ruhenden Gas schneller nach unten fallt als ein einzelnes 
Teilchen. Weiter ist zu beachten, daB in der Praxis die Gemische meistens aus Teilchen verschie- 
dener Gro8e und unregelmaBiger Form bestehen. Man kann daher die Beziehungen zur Berechnung 
des Luftwiderstandes von Einzelkérpern nicht ohne weiteres auf Gemische iibertragen und muB 
versuchen, die Theorie des Widerstandes eines Gemisches den hier vorliegenden Verhdltnissen an- 
zupassen. 

Die vorliegende Arbeit befaBt sich daher mit der theoretischen und der experimentellen Be- 
stimmung des Luftwiderstandes von kérnigen und staubférmigen Giitern im Luftstrom. Da der 
Luftwiderstand von kérnigen und staubférmigen Giitern mit zunehmender Feinheit immer gréBer 
wird, kann von dem Luftwiderstand auch auf die Feinheit des Gutes und umgekehrt von der Fein- 
heit des Gutes auch auf den Luftwiderstand geschlossen werden. Die Bestimmung des Luftwider- 
standes ist daher gleichbedeutend mit einer Feinheitsbestimmung. 


2. Theorie des Luftwiderstandes von kérnigen und staubformigen Giitern. a) Das Einzel- 
teilchen. Die Luftwiderstande von einzelnen einfachen Kérpern, z. B. der Kugel, der Kreis- 
platte, sind bereits seit langer Zeit Gegenstand weitgehender experimenteller und rechnerischer 
Untersuchungen gewesen” *. Dabei kann diese Aufgabe fiir die Kugel als gelést betrachtet werden 
und es ist heute mdglich, z. B. den Widerstand einer Kugel in irgendeinem strémenden Medium 
mit unbedingter Zuverlassigkeit zu berechnen. Dies gilt auch fiir ganz kleine Kugeln von der 
GréBenordnung der Teilchen eines Staub—Luft-Gemisches. 

Die Berechnung des Luftwiderstandes W von Einzelkérpern erfolgt nach der Beziehung 


Weg rae, (1) 
28 


* Gekiirzte Fassung der gleichnamigen von der Fakultat fiir Maschinenwesen genehmigten Dissertation 
des Verfassers, Technische Hochschule Karlsruhe, 1957. 

1 W. Barth, Z. angew. Math. Mech. 30 (1950) S. 255. 

2 W. Barth, Ing.-Arch. 16 (1948) 5. 147. 

3 G. Weidner, Grundsatzliche Untersuchungen iiber den pneumatischen Fordervorgang, insbesondere tiber 
die Verhiltnisse bei Beschleunigung und Umlenkung. Diss. Karlsruhe 1954 und Forschung auf dem Gebiete 
des Ingenieurwesens 21 (1955) S. 145. 

4 L. Prandtl, Fiihrer durch die Stromungslehre, $. 177. Braunschweig 1949. 

5 E. Eiffel, Recherches expérimentales sur la résistence de l’air exécutées — la Tour Eiffel, S. 78, Paris 1907. 
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Hierbei bedeuten F' die Stirnflache des Kérpers in Anstrémrichtung, w die mittlere Anstrém- 
geschwindigkeit, y, das spezifische Gewicht des strémenden Mediums und « einen Widerstands- 
beiwert, der von der Reynoldsschen Kennzahl 

dy, w 


Re = 
y 


abhangt; d, ist eine charakteristische Kérperlange, » die kinematische Zahigkeit des strémenden 
Mediums, g die Erdbeschleunigung. 

Fiir die Kugel ist in Abb. 1 die GréBe a als Funktion der Reynoldsschen Kennzahl Re dar- 
gestellt. Bei Reynoldsschen Zahlen Re<1 ist « umgekehrt proportional zu Re. Das bedeutet, 
daB nur die Zahigkeitskrafte des umstrémenden Mediums den Widerstand bedingen. Bei Re> 2000 
wird der Wert « nahezu 
konstant. In diesem Fall 
bestimmen ausschlieBlich 
die Tragheitskrafte den 
Luftwiderstand des Kér- 
pers. Dazwischen liegt 
ein Ubergangsgebiet, in 
dem beide Einfliisse wirk- 
sam sind. Alle wbrigen 
Korper zeigen grundsatz- 
lich ein ahnliches Ver- 
halten, wenngleich sich 
die absoluten Werte er- 
heblich —_unterscheiden 
} kénnen. 

b) DasGemisch. Im 


Gegensatz zu der groBen 


a Zahl der Untersuchungen 

: : E zo | &  iiber den Luftwiderstand 
10 10 0 0 0 0 0 0 ow von Einzelkérpern fehlen 
Abb. 1. Abhangigkeit des Kugelwiderstandsbeiwertes « yon der Reynoldsschen Zahl Re. bis heute fast vollig Un- 


terlagen itiber den Wider- 
stand von Gemischen von Feststoffen in Gasen und Fliissigkeiten. Zum Aufbau einer Theorie 
des Luftwiderstandes von Gemischen miissen daher zunichst einige grundlegende Betrachtungen 
angestellt werden. 

Man wird dabei zweckmaBigerweise so vorgehen, daB8 man versucht, die gleichen Uberlegungen 
wie sie fiir die Kinzelkugel gemacht wurden, sinngemaB auf ein Gemisch zu iibertragen. Die Form 
der angegebenen Beziehung (1) diirfte jedoch fiir die Bestimmung des Luftwiderstandes von Gut- 
gemischen nicht geeignet sein. Das Gemisch besteht im allgemeinen aus kleinen Kérperchen, deren 
Gestalt nicht genau angegeben werden kann. Es ist daher nicht méglich, den Luftwiderstand wie 
beim Einzelkérper auf die Stirnflache zu beziehen, ganz abgesehen davon, daB jedes Einzelteilchen 
eine andere Stirnflache besitzt. Daher scheint es zweckmaBig zu sein, bei einem Korngemisch den 
Luftwiderstand nicht auf eine geometrische GréBe zu beziehen. Im Gegensatz zur Stirnflache ist 
das Gewicht des Gemisches in der Regel bekannt oder kann leicht bestimmt werden. Das Problem 
vereinfacht sich dann auf die Bestimmung einer GréBe n, die wie folgt definiert ist: 


W _ Luftwiderstand 


are Q = Gutgewicht ~ (2) 


Der Wert n ist eine dimensionslose Zahl, die angibt, wievielmal gréBer der Luftwiderstand W 
des Gemisches als das Gewicht Q des Gemisches ist. Beim freien Fall im ruhenden Gas wird n=1. 
Nachstehend soll nun versucht werden die Beziehung (1) so umzuformen, daB sie auch auf 
Gemische angewendet werden kann. Es werde die Annahme gemacht, daB das Gutgemisch aus 
Kugeln vom Durchmesser d,, und vom spezifischen Gewicht y, bestehe, und daB sich die Teilchen 
gegenseitig nicht beeinflussen mégen. Der Gesamtwiderstand des Gemisches ist damit gleich der 
Summe der Luftwiderstande der einzelnen Teilchen. Beachtet man, daB das Gewicht Q einer Kugel 
ze dé 


Q=n5 
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ist, dann ergibt sich aus (1) und (2) folgender Ausdruck: 
37, w 
ay a dee (3) 
Diese Gleichung wurde ausgewertet und in Abb. 2 iiber der Abszisse ny;,/Yo dann die relative 
Anstrémgeschwindigkeit w mit dem Kugeldurchmesser d,, als Parameter aufgetragen, wobei Vo 
das spezifische Gewicht des Wassers von 1000 kg/m? bei 4° C ist. Der Mittelwert der Anstrom- 
geschwindigkeit w allerTeilchen 
des Gemisches errechnet sich 
aus der Differenz der Luft- 
geschwindigkeit v und der Gut- 
geschwindigkeit ¢ zu 
w=vu—ce. 


7o{m/sek| 


Das Diagramm Abb.2 gilt nur 
fir Luft von 20° C und 760 mm 
Hg. Um dieses Diagramm auch 
bei anderen Betriebszustanden 
anwenden zu kénnen, wird die 
Gleichung (3) wie folgt umge- 


=| 
formt. Man setzt 
v Re 
————— 4 
aeaies (4) fi 
git , ee ED 
und 10 0 19 10° 0° 10° uy 
7 ?, 
Co —— a — ut " (5) Abb. 2, Zusammenhang zwischen Luftwiderstand, Anstrémgeschwindigkeit und Durchmesser 
Q Y kugeliger Teilchen, 
Ws, gelig 
. . 5 TU: 1 
Hierbei bedeuten 7 die absolute cae 


Zahigkeit des Férdermediums und 
o die Dichte des Férdermediums. 
Die Beziehungen (4) und (5) er- 
geben 


Fis (6) 
wy, 


Erweitert man nun (3) auf beiden 


Seiten mit ’”Z und ersetzt d,, 


04 
nach Beziehung (6), dann ergibt 
sich 


Brute: 3 Ee a 


igh = 
Yon? Agro 


7 w “2 dr, 
) g 
(7) 
Ue Yr Pang 
In (7) kann (n ne } als die ab- % 5 gt 7) Pr 7 oF Hee 
; Yr. ; Abb. 3. Zusammenhang zwischen Luftwiderstand, Anstrémgeschwindigkeit und 
hangige und |w a als die unab- Durchmesser kugeliger Teilchen. 
hangige Verdnderliche angesehen ; 
1s : Ve Vr ral ee 
werden. Durch Auswerten der Beziehung (7) und Auftragen von (n ee als Abszisse sowie | w a 
0 


als Ordinate entsteht Abb. 3. Paraméter der Kurvenschar ist wieder der Kugeldurchmesser d,,. 
Dieses Diagramm ist fiir alle Medien anwendbar unter der Voraussetzung, daB das spezifische 
Gewicht des gasférmigen oder fliissigen Mediums klein gegeniiber dem spezifischen Gewicht des 
Materials ist. Fir den Fali, daB die spezifischen Gewichte y, und y,; von der gleichen GréBen- 
ordnung sind, bleibt das Diagramm das gleiche, nur hat die Einheit der Abszisse dann den Wert 


C52) Ft 
(ee Peale 
Yo 1] 
Diese Darstellung eignet sich vor allem dazu, eine Verbindung des Luftwiderstandes eines 
Einzelkérpers mit dem Luftwiderstand eines Gemisches herzustellen. Sind die Werte n und w 
228 
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fiir ein Gemisch beispielsweise durch den Versuch bestimmt, so kann mit Hilfe der Abb. 2 oder 
Abb. 3 der Durchmesser einer Kugel eines Kugelgemisches angegeben werden, das unter der Vor- 
aussetzung, da sich die Kugeln gegenseitig nicht beeinflussen, den gleichen Widerstand wie das 
Gemisch haben wiirde. Dieser Durchmesser mige als gleichwertiger Kugeldurchmesser d,, des 
Gemisches definiert werden. Durch Angabe des gleichwertigen Kugeldurchmessers, des spezi- 
fischen Gewichtes des Materials, des spezifischen Gewichtes und der Zahigkeit des Férdermediums, 
der Gasgeschwindigkeit und der Gutgeschwindigkeit ist somit der Widerstand des Gutgemisches 
eindeutig festgelegt. Dabei ist allerdings zu beachten, dafs der gleichwertige Kugeldurchmesser 
im allgemeinen keine Konstante ist, sondern daB er von den Betriebsverhaltnissen, insbesondere 
von der Materialbeladung bzw. der Materialdichte und anderen GréBen wie Kornform, Korn- 
zusammensetzung usw. abhangig sein kann. 

Statt des gleichwertigen Kugeldurchmessers d, kann es in manchen Fallen auch zweckmabig 
sein die dazugehérige Sinkgeschwindigkeit w, als MaB fiir den Luftwiderstand des Gemisches ein- 
zufiihren. Die Sinkgeschwindigkeit kann aus (3) leicht errechnet werden, bzw. aus den Abb. 2 
oder 3 entnommen werden, wenn man n=1 setzt. 


c) Ahnlichkeitsttberlegungen. Fiir ein homogenes Medium sind die Ahnlichkeitsbedin- 
gungen allgemein bekannt, die eingehalten werden miissen, sobald man von einem Strémungs- _ 
vorgang auf einen anderen schlieBen will. Auf den Fall jedoch, daB sich in dem strémenden Medium 
Festteilchen bewegen oder daB es sich um die Strémungen von Mehrstoffgemischen handelt, sind 
die bisher durchgefiihrten Uberlegungen im allgemeinen nicht anwendbar. 

Barth? hat versucht, auch fiir diesen Fall Ahnlichkeitsbedingungen aufzustellen und kommt 
zu folgenden Ergebnissen: Wenn zwei Strémungsvorgange ahnlich verlaufen sollen, miissen fol- 
gende Bedingungen erfiillt sein: 


1. Geometrische Ahnlichkeit zwischen Modell und GroSausfihrung, 


Tragheitskrafte vdy vd 


* Zahigkeitskrafte oe ZN yo BE UIE, 


Tragheitskrafte v 


Se panE epse ae Vda = Fr = konst. (Froudesche Kennzahl). 


Bei Mehrstoffgemischstrémungen ist auBerdem noch den weiteren Forderungen 


ee konst., 
V « 

5. —| = konst. 
Ny 


zu geniigen. Neben den bereits bekannten Bezeichnungen bedeutet hier d eine charakteristische 
Lange und v eine charakteristische Geschwindigkeit. Durch einen Stern als Index soll eine zweite 
charakteristische GréBe des gleichen Strémungsfeldes gekennzeichnet werden. 


Die gleichzeitige Erfillung aller dieser Gesetze fiir eine Mehrstoffgemischstrémung ist in den 
meisten Fallen unméglich'. Man muf deshalb haufig einen Ausweg in der Weise suchen, daB man 
an Stelle der strengen Ahnlichkeitsbedingungen Ersatzbedingungen verwendet. Diese Ersatz- 
bedingungen verzichten auf die Einhaltung der geometrischen Ahnlichkeit hinsichtlich der Teilchen 
im Strémungsfeld. Dagegen muB die Bedingung erfiillt sein, da® die auf die Teilchen wirkenden 
Krafte im gleichen Verhaltnis zu den iibrigen Kraften in der Strémung stehen. 


In einer eindimensionalen Strémung gilt fiir ein Teilchen mit der Masse m 
W—Q=m—. 


Driickt man hierin W durch (1), sowie Q durch die Fallgeschwindigkeit nach Abschnitt b) aus, 
und legt man dem Widerstandsheiwert a die GesetzmaBigkeit 


__ konst. konst. 


ue (ean 


J 


1 W. Barth, ZVDI, 92 (1950) S. 105. 
2 W. Barth, Chemie-Ingenieur-Technik, 26 (1954) S. 29. 
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zugrunde, so erhalt man die Bewegungsdifferentialgleichung 
g (oe)? — gute a Fyn, 0) 


Hieraus lassen sich die folgenden Ahnlichkeitsbedingungen ableiten: 


wy my 
rae —= konst.. (10) 
we 2—-x 
(“*) = konst. (11) 


Die Bedingungen (10) und (11) gelten fiir ein Einzelteilchen. Fiir ein Gemisch tritt als weitere 
Bedingung hinzu, da die Materialbeladung yw, d.i. das Verhaltnis des geférderten Gutgewichts 
je Zeiteinheit G zum geférderten Luftgewicht je Zeiteinheit Vy,;, konstant bleibt, also 


G 
= Vy, = konst. (12) 


Der gleichwertige Korndurchmesser d,, der als Maf fiir den Luftwiderstand des Gemisches 
aufgefaBt werden kann, muf daher als Funktion dieser AhnlichkeitsgréBen dargestellt werden 
kénnen. Um dimensionslose Werte zu erhalten, wird man daher das Verhaltnis d,/d einfiihren. 


Es mu8 sein 
dk ce we ™ 
Fo \ es |. (13) 


v 


Wenn der Luftwiderstand eines Gemisches bei einem Betriebszustand bekannt ist, kénnen mit 
Hilfe der oben angegebenen Beziehungen fiir eine Reihe anderer Betriebszustande die Luftwider- 
stande angegeben werden. Ks sind zwei ausgezeichnete Falle zu unterscheiden: 


I) x =1. Hier wird 


w,.U Ww di, 
e/a inet ae one eee Ley ed illles 
ae = konst., = = konst., 1 = Konstsr, i= 0 |u ie 


Diese Bedingungen gelten fiir den Fall, daB das Stokessche Widerstandsgesetz anwendbar ist. 
II) x = 0. Hier wird 


dg v? d (4 ; dg’ ov 
Diese Bedingungen gelten fiir den Fall, daB das quadratische Widerstandsgesetz erfiillt ist. 


Fir die beiden Falle gilt auch v?/gd = konst., wie aus den Gleichungen leicht abgeleitet werden 
kann, und dies ist nichts anderes als die Froudesche Kennzahl. 


w? w? wz w 
-/ — konst., + — konst., fs = konst. , aie @ f t ; 


3. Die Bestimmung des Luftwiderstandes von kérnigen und staubformigen Giitern. Eine rein 
rechnerische Ermittlung des Luftwiderstandes von Gemischen ist im allgemeinen genau so wenig 
méglich, wie es bei Einzelkérpern der Fall ist. Erst wenn fiir das Gemisch eine geeignete Dar- 
stellung des Luftwiderstandes vorhanden ist, kénnen tiber die Ahnlichkeitsbeziehungen hinweg 
Aussagen fiir einen gegebenen Betriebszustand gemacht werden. Es kann daher nur eine experi- 
mentelle Ermittlung des Luftwiderstandes in Frage kommen. Zu seiner Bestimmung kénnen 


folgende Verfahren Verwendung finden. 


a) Windsichtung?? und Sedimentation’. Von H. W. Gonell wurde ein Gerat beschrieben, 
das zur Widerstands- und Feinheitsbestimmung von Gemischen Verwendung finden kann. Mit 
diesem Gerdt kann die Sinkgeschwindigkeit der einzelnen Gutteilchen des Gemisches in einem 
stromenden Medium bestimmt und die Zusammensetzung des Gemisches hinsichtlich der Sink- 
geschwindigkeit der einzelnen Gutteilchen ermittelt werden. In Abb. 4 ist ein Gerat dargestellt, 
das nach diesem Prinzip arbeitet. 

Die Sichtung des untersuchten Gutes erfolgt dadurch, da unter den konischen Ansaugstutzen 
eine mit Gut gefiillte Schale gehalten wird. Die Ansaugluft reift alle Teilchen mit sich, deren 


1 W. Barth, Fette und Seifen 55 (1953) S. 617. 

2 VDI-Richtlinien, Leistungsversuche an Entstaubern, 2. Ausgabe 1949. Deutscher Ingenieur-Verlag 
GmbH. Diisseldorf. 
3 H. W.Gonell, ZVDI 27 (1928) S. 945. 


394 W. Rausch: Luftwiderstande von kérnigen und staubformigen Giitern im Luftstrom —_Ingenieur-Archiv 


Gewicht bei der eingestellten Luftgeschwindigkeit kleiner als ihr Luftwiderstand ist. Alle ibrigen 
Teilchen bleiben in der Schale oder fallen in sie zuriick. Fiir jede eingestellte Luftgeschwindigkeit 
wird der Gutriickstand in der Schale gemessen. Aus der Luftgeschwindigkeit kann auf den Luft- 
widerstand bzw. die KorngréBe geschlossen werden. 

Bei sehr feinkérnigem Gut treten bei der Windsichtung Schwierigkeiten auf. Die Luftgeschwin- 
digkeiten werden sehr klein und damit schwerer mefbbar. Ein Festhaften der untersuchten Kérner 
an der Wand des Rohres und Zusammenbacken des Gutes verfalschen oft das Ergebnis. In manchen 
Fallen verwendet man deshalb zweckmafigerweise statt der Windsichtung die Sedimentation. Am 
besten hat sich die Methode nach Andreasen! bewahrt. Das Gerat nach Andreasen zeigt Abb. 5. 

Mit Hilfe dieser beiden Verfahren kénnen nur Aussagen iiber den Luftwiderstand bei einer An- 
stromgeschwindigkeit, d. h. bei der Sinkgeschwindigkeit (n = 1) gemacht werden. Es besteht keine 
Méglichkeit die Zusammenhange zwischen Luftwiderstand, 
Materialbeladung und Anstrémgeschwindigkeit aufzuzeigen. 
Gerade diese Zusammenhange sind aber fiir den Ablauf von 
Strémungsvorgangen bei Gasen und Festteilchen von grund- 
legender Bedeutung. 


— 
——— 
—— 
— 


aS — 


eas utgeschwindigkelt c 


ie Vi Abb. 5. Apparat 
Abb. 4, Windsichtapparat. nach Andreasen. Abb. 6. Geschwindigkeits- und Druckverlauf in der MeBdiise. 


b) Verfahren nach Barth. Die beschriebenen Nachteile der Sedimentation und der Wind- 
sichtung werden ausgeschaltet in einem Verfahren nach Barth?. Hier wird der Luftwiderstand des 
Gutes im freien Luftstrom gemessen. Damit kénnen beliebige Anstrémverhdltnisse erreicht werden 
und auBerdem kann der Einflu8 der Materialdichte leicht studiert werden. Wegen seiner Bedeue 
samkeit fiir die gestellte Aufgabe soll dieses Verfahren eingehender betrachtet und die auf diese 
Weise gewonnenen MeBwerte den Untersuchungen zugrunde gelegt werden. 

1. Beschreibung des Verfahrens. Die Grundlagen des Verfahrens sollen an Hand von Abb. 6 
studiert werden. In eine Rohrleitung vom Durchmesser dy wird eine MeBdiise eingebaut, die aus 
einer zylindrischen Bohrung vom Durchmesser d, baw. d, und einer gewissen Lange 1 ae einem 
konischen Ubergangsstiick am vorderen Rohrende (in Strémungsrichtung gesehen) besteht. 

Durch die Rohrleitung mége ein Staub—Luft-Gemisch gesaugt und die statischen Driicke an den 
MaBebenen 0—0, 1—1 und 2—2 bestimmt werden. Der Einfachheit halber soll angenommen werden 
daB alle Teilchen des Gemisches die gleiche Gré®e haben. Hinter der MeBdiise wird man nonaales 


1 Siehe FuGBnote 1 auf S. 323. 
2 Siehe FuBnote 2 von Seite 319. 
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weise zunachst ein Filter und anschlieBend eine MeBblende anordnen, um die durchgehenden Luft- 
und Gutmengen bestimmen zu kénnen. In der Diise selbst stellt sich grundsatzlich ein Geschwin- 
digkeits- und Druckverlauf ein, wie er in Abb. 6 aufgezeichnet ist. 

Betrachten wir ein Teilchen beim Durchlaufen des Rohres und der anschlieBenden Diise nach 
Abb. 6. Vor Eintritt in die Diise mégen Luftgeschwindigkeit v, und Teilchengeschwindigkeit c, 
annahernd gleich grof sein. Beim Kintritt in die Diise steigt die Luftgeschwindigkeit sehr schnell 
auf den Wert v, an, wahrend das Teilchen infolge seiner Tragheit zunachst seine Geschwindigkeit 
beizubehalten sucht. Dabei entsteht eine hohe Relativgeschwindigkeit w—=v—c, die eine starke 
Beschleunigung bewirkt. Eine quantitative Betrachtung zeigt, daf® in diesem Falle der Einflu8 
des Gewichtes und der Wandreibung vernachlissigbar klein sind. Erst allmahlich beim Durch- 
laufen der Bohrung erhéht sich die Teilchengeschwindigkeit auf c,. Die relative Anstrémgeschwin- 
digkeit, und damit auch die Beschleunigung, wird kleiner, so daB der EinfluB von Gewicht und 
Wandreibung immer mehr an Bedeutung gewinnt. Bei geniigend langer zylindrischer Diisenlange | 
stellt sich ein Beharrungszustand so ein, daB w=v—c=konst. wird und der Luftwiderstand nur 
noch von Gewichtseinfliissen und Wandreibungsverlusten verursacht wird. 

Zur Beschleunigung des Gutes mu auf die Rohrlange dl ein bestimmtes, zusatzliches Druck- 
gefalle dp, der Luft aufgewendet werden. Bezeichnet man den Querschnitt des zylindrischen Diisen- 
teils mit f, und macht man die Annahme, dafs die von den Teilchen auf die Wand ausgeiibten 
Reibungskrafte und das Gutgewicht vernachlassigbar klein sind, so ergibt die Anwendung des 
Impulssatzes 


dW = dm = dp.f, (14) 
de dc 
dm ae dm=e¢ = Opes 
Nun ist aber 


ip ne 
CAE: 


und damit wird 


G 
dpifis des 
Pf te 


Fihrt man die Materialbeladung y ein und integriert, so ergibt sich 
ip 
(ae ie 0 
ines 4 


4, 


Zu den aus Abb. 6 ersehbaren Bezeichnungen miissen fiir die folgenden Ausfiihrungen noch einige 
weitere eingefiihrt werden: 


Ap! 


2 
i) Wie 


ie. 


Druckverlustbeiwert des Diiseneinlaufs, 


le : Druckverlustheiwert der gesamten MeBdiise. 


Mit Hilfe der Beziehung (15) sollen nun die Druckgefalle Ap!! und Ap!, sowie die dimensionslosen 
Druckverlustbeiwerte &!! und &! ermittelt werden. Fir den Fall v = v, = konst. findet man 


, ey. Yr 
Apt = Ap? — Api =— = "s (a) = 87 Pa uO) 
i) Ope i 
aes) 
und ferner : 
L Yr Us Yr 
Ap) = Ap, — Ap, = Z Vg (ey — &) = Ee “Oe” (17) 


I ce 2 I pm 
= 2/2 aE 
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Die Annahme v = v, = konst. wurde der Einfachheit wegen auch fiir die konische Verengung 
(Bereich zwischen MeBebene 0—0 und 1—1) zugrunde gelegt, obwohl diese Voraussetzung in der 
konischen Verengung nicht erfiillt ist. Der dadurch entstehende Fehler kann vernachlassigt werden. 

Die Druckgefalle Ap!! und Ap! werden durch Messung bestimmt. Fiir die einzelnen MeBpunkte 
werden die KenngréBen &!, €!' und « bestimmt und €! baw. €!' in Abhangigkeit von wu graphisch 
aufgetragen. Aus dieser Auftragung kénnen €! und él fir 4 = konst. (z. B. uw = 2, 4, 6) entnommen 
und auf Grund der nachstehenden Beziehungen Schliisse auf die Geschwindigkeit des Gutes in der 
Diise gezogen werden. Es muf sein 


Cae at Co 
ge me oe, . nape? ie 
fre yer: 
Ups ae CPs é (19) 


Um eine Aussage iiber die GréBe von c, und c, machen zu kénnen, muB cy, bekannt sein. Beziiglich 
der Geschwindigkeit c) kann man sagen, daB vy > c) > 0 sein muh. Es ist dabei zu beriicksichtigen, 
daB ein Teil der Teilchen auf die Wand 
des konischen Einlaufs zwischen den MeB- 
ebenen 0—0 und 1—1 (Abb. 6) auftreffen 
wird, wobei durch den Aufprall ein Teil 
der Geschwindigkeit vernichtet wird. Ks 
besteht daher die Méglichkeit, daB die 
Geschwindigkeit c, wesentlich kleiner als 
v, wird, auch wenn in der Robrleitung 
vor der Diise c, und vy praktisch gleich 
groB angenommen werden kénnen. Um 
auf alle Falle klare Verhaltnisse zu er- 
halten, wurde in vielen Fallen ein Einlauf 
nach Abb. 7 verwendet. Bei dieser An- 
ordnung wurde vor der zylindrischen 
Diise ein Prallkegel angeordnet, auf den 
Jii\!|\) die ankommenden Gutteilchen auftreffen. 
Sse SCDabei wird die Geschwindigkeitsenergie 
der Gutteilchen vernichtet, so daB cy prak- 

tisch gleich 0 gesetzt werden kann. 


= Der Verlauf der Gutgeschwindigkeit 
eu i in der zylindrischen Bohrung zwischen 
Filter i den MeBebenen 1—1 und 2—2 kann da- 


durch bestimmt werden, daB man die 
Druckdifferenzen Ap!" fir Diisen mit 
verschiedenen Lingen / ermittelt. Man 
kann die GréBe c, am Diisenende aus (15) 
bzw. (18) und (19) bestimmen und damit den Verlauf der Gutgeschwindigkeit in Abhangig- 
keit von der Diisenlange festlegen. 

Zur Bestimmung der Werte n und w wird nun die Diisenlange in einzelne Intervalle von der 
Lange Al eingeteilt und hierfiir die GréBen n und w errechnet. Diese Intervalle werden zweck- 
maRigerweise so gewahlt, daB sie der Langendifferenz zweier in ihrer Lange aufeinanderfolgender 
MeBdiisen entsprechen. Fiir die Berechnung des Wertes n kann man von (2) ausgehen. Unter den 


getroffenen Voraussetzungen der Vernachlassigbarkeit von Wandreibungskraften und Gutgewicht 
bekommt man mit 


16° 
Abb. 7. Anordnung der MeBdiise. 


W =m— und Oe Tine 


den Wert 


_lid _1lddl Cmae 
pdb = edie (20) 


Fir die numerische Berechnung geht man zu Differentialquotienten iiber und erhalt 


ee ie Som 
eg Al gAl ee 


n 


(21) 
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Der Wert c,,,/v, fiir ein Intervall wird dadurch ermittelt, da8 man, ausgehend von der angenommenen 
Geschwindigkeit c) am Anfang der Diise (bei Einbau des Kegels nach Abb. 7 wird cy = 0) die Ge- 
schwindigkeitszunahmen in den einzelnen Intervallen addiert und fiir das untersuchte Intervall 
eine Zunahme im Mittel von + Ac annimmt. Aus den Werten Ac/v, und c,,/v, wird dann (n V4/Vo)s 
w und die Materialdichte w® errechnet. Die GréBe ® ist definiert als 


Oye urs 
pe = w—, 


und stellt ein Maf fiir den gegenseitigen Teilchenabstand dar. Die gefundenen MeBwerte werden 
in ein (n y;/7o, w)-Diagramm eingetragen, aus dem der Widerstand des Gemisches, ausgedriickt 
durch die GréBe des gleichwertigen Kugeldurchmessers 


Hh d,, abgelesen werden kann. Durch Auftragen der ein- 
aie zelnen MeBwerte kann man ein Bild iiber den Luft- 
TAA widerstand des Gutes unter den verschiedenen Betriebs- 


bedingungen gewinnen. Dabei ist zu beachten, daB 
durch die Aufteilung in einzelne Intervalle gewisse 
Fehler bei der Auswertung entstehen kénnen, die aber 
aus der Lage der einzelnen Punkte zu erkennen sind 
und gegebenenfalls leicht korrigiert werden kénnen. 


ay 
i 
= 
WEY 
| | & [Mtertopr 
imal 
2-28 fll, —o¢ 
, iit . we 
SZ = 
Bees 
atts DruckmeBstellen\ = = 
ore 
Uy 
UT] 


Abb. 8. Gerat zur Bestimmung des Luftwiderstandes Abb. 9. MeSblende zur Luftmengenmessung. 
von Staéuben und Mehlen. 


2. Beschreibung des Mefgerdtes. Zur Durchfiihrung der beschriebenen Untersuchungen wurde 
das in Abb. 8 dargestellte Gerat entwickelt. Das zu untersuchende Gut wird in den Aufgabe- 
behalter J eingefiillt. Den Boden dieses Aufgabebehilters bildet ein Drehteller 2, der durch einen 
KurzschluBlaufermotor angetrieben wird. Die Drehzahl des Drehtellers kann durch Anderung 
der Ubersetzung mit Hilfe von Stufenscheiben verandert werden. Durch Verschieben der Behilter- 
wand in vertikaler Richtung wird ein Spalt freigegeben, durch den das Material aus dem Aufgabe- 
behalter J ausflieBt. Um einen gleichmaBigen AusfluB zu erreichen, wurde auf dem Drehteller ein 
Konus mit Seitenarmen befestigt, der sich mit dem Teller dreht und in dem Behalter gewissermaBen 
als Rithrwerk wirkt. Das Material wird durch einen zur Feinregulierung der Aufgabemenge ver- 
stellbaren Abstreifer 3 vom Drehteller abgestreift und fallt in den Trichter 4, der um eine auBerhalb 
der Symmetrielinie liegende Achse geschwenkt werden kann. Das Material fallt entweder in das 
AuffanggefaB 5 oder in das Zufihrungsrohr 6 und wird damit der MefBdiise zugefiihrt. Dieses 
Zufiihrungsrohr 6 zur MeBdiise hat an einem Ende eine konische Erweiterung, durch die Luft aus 
‘der Atmosphare angesaugt wird und sich mit den Festteilchen mischt. Das Gemisch strémt durch 
die MeBdiise 7 und gelangt in ein Stoffilter, in dem das Material zuriickgehalten wird. Die ge- 
reinigte Luft durchflie8t die MeBblende 8, anschlieSend das Geblase 9 und gelangt dann ins Freie. 
Das Filter ist in den Filterkasten 10 eingebaut. Durch Lésen einiger Fliigelschrauben kann der 
Deckel des Filterkastens entfernt und das Filter ausgebaut werden. Die MefBdiise ist in Abb. 7 
dargestellt. Sie sitzt direkt auf dem Filterdeckel, der auch eine MeBbobrung zur Messung des 
Druckes hinter der Diise enthalt. In die MeBdiise ist aus den auf Seite 326 dargelegten Griinden 
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ein Kegel vor der Diise eingebaut worden, der durch Entfernen des Einsatzes a Abb. 7 ausgebaut 
werden kann, so da die Diise auch ohne Kegel betrieben werden kann. Die Bohrung fiir die Druck- 
an (nms) meBstelle vor der Diise miindet in den 
ae Raum zwischen Diise und Zufiihrrohr 6 zur 
ye MeBdiise. Die Mef&blende 8 ist in Abb. 9 

of dargestellt. Sie wurde wie die Mefdiise 7 

a eid a BA sleichfalls in den Filterkasten eingebaut. 
Die eine Druckentnahmestelle fiir die MeB- 
blende miindet in den Filterkasten, die 
andere in das an die MeBblende anschlie- 
Bende Rohr. Das Gerat wurde so ausge- 
fiihrt, daB alle Teile leicht ausgebaut bzw. 
ausgewechselt werden kénnen. Die Diffe- 
renzdriicke bzw. die Wirkdriicke der MeB- 
diise 7 und der MeBblende 8 werden mit 
wassergefillten U-Rohren gemessen. Die 
MeBblende 8 wurde mit einer Gasuhr ge- 
eicht. Die Einstellung der Luftmenge er- 
folgt durch Drosselung des Regulierschie- 
bers 12 oder durch Anderung der Drehzahl 
des Geblases 9. Das Geblase 9 wird von 
einem Universalmotor angetrieben; die 
Geblasedrehzahl kann deshalb durch Vor- 
schalten eines Ohmschen Widerstandes ein- 
fach geregelt werden. Das Geblase ist zwei- 
stufig und hat eine Férderhéhe von etwa 
500mm W.S. bei einem Durchsatz von 
6l/sec. bei Normallast und maximaler 


Drehzahl von 3000 U/min. 
[Pat Tent ane vee as ee ee) 3. Versuchsdurchftihrung und Versuchs- 


400 


Diise A, l=20mm 


0 7 Zi 3 . i 
ergebnisse. Die Durchfiihrung der Messun- 
y 8 
Abb, 10, Druckverlust 4p, in Abhingigkeit von der Gutmenge G. gen erfolgt in der Weise. daB man zunachst 
> 
i mit dem Gerat reine Luft ansaugt. Es wird 
a7 dann die Aufgabevorrichtung in Betrieb 


genommen und mit Hilfe des Abstreifers 

eine bestimmte Gutmenge eingestellt. Durch 

prea Anderung des Widerstandes im Stromkreis 
5 des Motors stimmt man die Drehzahl des 


oe | ese Absaugegeblases so ab, da8B eine bestimm- 
i Se gam | ae te, gewiinschte Luftmenge angesaugt wird. 

— vréiybek YO ee Durch kurzzeitiges Schwenken des Trich- 
43 3 : : Lisentinge 


oy, 
Vp 352; 43m/sek 1-0 he Me 


7 


Lufigeschwindigkeit v,=362m|sek | __— Vn 
U,=35,2m./sek 0,3 
x i : _} Disenlange arte 
Z-Omm, 
02 ee ; 
Soe I tel eli | 


) 6 7 8 
ie 
Luftgeschwindigheit or Jo —== 
© V,= 352m/Sek x 
x v,=43 m/sek Z{mm| 


4 U,=51 m/sek 0 10 20 30 40 50 60 
Abb, 11. Druckverlustbeiwert ¢, in Abhingigkeit yon der Material- Abb. 12. Geschwindigkeitsverhiltnis c/v, in Abhangigkeit yon 
beladung 4. Material: Glaskugeln. der Diisenlinge 1. Material: Glaskugeln. 


ters 4 fallt das Gut in das AuffanggefaB 5. Mit einer Stoppuhr wird die Aufgabezeit gemessen 
und die aufgefangene _Gutmenge durch Abwiegen des AuffanggefaBes 5 bestimmt. Diese Gut- 
mengenbestimmung wird vor und nach jeder Ablesung des Wirkdruckes an der MeBblende und 


TU 
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an der MeBdiise vorgenommen. Es ware natiirlich grundsatzlich auch méglich, durch Auswiegen 
des Filters im Filterkasten 10 die aufgegebene Gutmenge zu bestimmen, jedoch hat das andere 
Verfahren, insbesondere bei vielen Messungen, den Vorteil der Einfachheit und kurzen Zeit- 
dauer. Die Messungen erfolgen in der Regel in der Weise, daB fiir eine konstante Gasmenge 
bzw. fiir einen konstanten Blendenwirkdruck Ap,, die Differenzdriicke Ap! baw. Ap" bei ver- 
schiedenen Gutmengen bestimmt werden. Die gemessenen Werte Ap! und Ap" werden dann 
zweckmaBigerweise in Abhangigkeit von G aufgetragen. Diese Messungen werden fiir Diisen mit 
verschiedenen Langen, z. B. 1 = 5, 10, 20, 40 und 80 mm, durchgefihrt. 


{m/sek| 


gig fir. =60 


gilltig fur u*=30 


-gilig Fr Bf 


_— 


Die Zahlen be/ den eingehragenen Punkten 
geben die Grobe von Lu* an 


o{a| a] v,=51 m/sek 
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200 500 400 500 600 700 800 500 1000 
Abb. 13. Diagramm nach Abb. 2 mit aus Versuchen gewonnenen Punkten, Material: Glaskugeln (Durchmesser 90 bis 140 1). 


Zur Nachpriifung der durchgefiihrten Uberlegungen wurden Gemische zweier verschiedener 


-Stoffe, namlich Glaskugeln und Aluminiumsilikat (Kontaktstoff), in dem beschriebenen MeBgerat 


untersucht. 

«) Glaskugeln. Die Glaskugeln sind im Handel erhaltlich. Sie haben einen Durchmesser 
zwischen 90 und 140 4 und ein spezifisches Gewicht y, = 2600 kg/m’. Bei der Verwendung dieser 
Glaskugeln bestand die Gewahr, daB sich das bali 


Gemisch aus Teilchen etwa gleicher GréBe zu- 

sammensetzt. AuBerdem ist auch die Form und zy Zs 

der Luftwiderstand des einzelnen Teilchens ge- A aliailte 

nau bekannt. Die Voraussetzungen der Rech- ELEIE 

nung sind damit im vorliegenden Fall erfiillt, ; 8 8 8 

so daB ein Vergleich der Rechnung mit den wpe 

Versuchsergebnissen méglich ist. es ee 
Die Untersuchungen wurden mit der Mef- [ v,-a52m/sek 

diise (Abb. 7) mit Kegeleinsatz durchgefihrt. v,-43 m/sek 

Es wurden Mefdiisen mit den Lingen | von 0, | 


10, 20 und 40 mm verwendet. In Abb. 10 * wurde 
der Druckverlust Ap, in Abhangigkeit von G fiir 
v = konst. fiir die verschiedenen Diisen auf- i ae 30 40 50 60 0 80 
getragen. Es zeigt sich, daB die Kurven nahezu Abb. 14. Ereiolncrtiges Eine rid rolnesiee dj, in Abhingigkeit von der 
linear mit G. ansteigen. Aus den MeBwerten wutr- Materialdichte p®, Material: Glaskugeln von 90 bis 140 yw. 

de é! baw. é! ermittelt und in Abhangigkeit von yw in Abb. 11* dargestellt. Aus dieser Abbildung 
wurden dann die &! und £!!-Werte fiir 1 = 2, 4 und 6, sowie fiir die Geschwindigkeiten v, = 35,2, 43 


* Diese Abbildung gibt als Beispiel den charakteristischen Verlauf wieder. Kurven fiir weitere Luft- 
geschwindigkeiten und vollstandige Zahlentafeln siehe die Dissertation des Verfassers. 
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und 51 m/sek entnommen und unter der Annahme, daB cy = 0 und cites =, gesetzt werden kann, 
cy = + &} bestimmt. In Abb. 12 * wurde der Verlauf des Geschwindigkeitsverhaltnisses c/v, tiber 
der Diisenlange aufgetragen. Er ergibt sich aus der Summe der Geschwindigkeitssteigerungen Ac 
in den einzelnen Intervallen. 

Die Lage der einzelnen MeBpunkte im (n, w)-Diagramm zeigt Abb. 13. Es ist verstandlich, dafs 
die MeBpunkte eine gewisse Streuung aufweisen werden. Der Luftwiderstand des Gemisches liegt 
zwischen einem gleichwertigen Korndurchmesser d, = 110 4 und d, = 200 yu. Versucht man die 
Punkte gleicher Materialdichte miteinander zu verbinden, so zeigt sich, daB mit zunehmender 
Materialdichte der Luftwiderstand abnimmt. Um diese Abhangigkeit besser erkennen zu kénnen, 
wurde in Abb. 14 der gleichwertige Kugeldurchmesser d, in Abhangigkeit von der Materialdichte pe 
aufgetragen. Bei der Interpolation auf ~® = 0 liegt der gleichwertige Kugeldurchmesser bei etwa 
110 w. Eine Nachpriifung der Kugeln ergab, dafs dieser Wert etwa dem wirklichen mittleren Kugel- 
durchmesser des Gemisches entspricht. 


7o{m/sek| 
50 Va | 
s: J 
Ze ( 
ve | 
‘His BPS 
Ye giiltig ture" =20 | ae a 
| a 
WT) Ko ih anes hey Hy =a 
vis ees hl i Agility tir. =10 
4 : ee Bie 
pi 268 
rae er oe 
 36y ee 
30 e ie 
ae ey 676 Noes, 
7 6 pe eae firgs*=5 
6,95 
: -giilig fir p*=2.5) 
20 ¥ E 
dy 25h / = 
4 8 
if Ary ae 8 
(iy thn nya £8 = 
10 - Serr, ol D2= 41 m/sek 
ae 
. Ke 
7 d= 64h Le 2 v,=20 m/sek 
a pate a Die Zahlen bei den eingetragenen Punkten | _, de 
7 dy = 5G geben die Griife von pan i 
100 200 300 400 500 600 700 800 900 


Abb. 15. Diagramm nach Abb.2 mit aus Versuchen gewonnenen Punkten. Material: Kontaktstoff (Aluminiumsilikat). 


Eine Abhangigkeit des gleichwertigen Kugeldurchmessers d, von der Luftgeschwindigkeit v 
14Bt sich nicht feststellen. Dies bedeutet, daf die fiir die Einzelkugel giiltigen Abhangigkeiten in 
bezug auf die Anstrémgeschwindigkeiten auch fiir das Gemisch gelten. Beziiglich des Luftwider- 
standes des Gemisches in den einzelnen Intervallen ist festzustellen, daB in dem Langenintervall 
von 0 bis 10 mm und von 20 bis 40 mm etwa die gleichen Werte ermittelt wurden, wahrend sich 
fiir das Intervall von 10 bis 20 mm etwas geringere Werte ergaben. Diese Unterschiede sind wahr- 
scheinlich auf MeBungenauigkeiten bzw. auf kleine Unterschiede in den Abmessungen der Diisen 
zurickzufiihren. Ks erscheint daher zulassig zu sein, Mittelwerte zu verwenden, zumal die Korrektur 


eines Wertes infolge des Zusammenhanges der einzelnen MeBgréBen untereinander auch die anderen 
Werte beeinflussen muf. 


6B) Kontaktstoff (Aluminiumsilikat). Die gleichen Untersuchungen wie mit den Glas- 
kugeln wurden mit dem Kontaktstoff (Aluminiumsilikat), ausgefiihrt. Sein spezifisches Gewicht 
ist y, = 1500 kg/m*. Das Gemisch besteht aus Kérnern verschiedener GréBe. Die Messung und 
Auszahlung einer geniigend groBen Anzahl unter dem Mikroskop ergab einen mittleren Durch- 
messer des Gutes von etwa 70 uw. Die Messungen wurden fiir drei verschiedene Geschwindigkeiten 
(v, = 41, 31,5 und 20 m/s) durchgefiihrt. Wie bei den Versuchen mit den Glaskugeln wurde auch 


hier die MeBdiise nach Abb. 7 mit dem Prallkegeleinsatz verwendet. Die Diisenlingen betrugen 0, 
10, 20, 40 und 80 mm. 


* Siehe FuBnote * auf S. 329. 
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Die Auswertung erfolgt in der gleichen Weise wie bei den Messungen mit den Glaskugeln. Die 
Kurvenblatter 15 und 16* stellen Versuchsergebnisse dar. Aus der Lage der MeBpunkte ergibt 
sich, daf$ in diesem Fall auch eine Abhangigkeit des gleichwertigen Kugeldurchmessers von der 
Luftgeschwindigkeit besteht. Fiir alle MeBpunkte, die zur gleichen Luftgeschwindigkeit gehéren, 
wurde ein Mittelwert gebildet und in Abb. 16 eingetragen. Man erkennt, da8 mit abnehmender 
Luftgeschwindigkeit v der Luftwiderstand abnimmt. Die auftretenden Streuungen der MeBwerte 
bei den einzelnen Luftgeschwindigkeiten diirften auf Ungenauigkeiten der Messungen und der 
Me8einrichtungen zuriickzufiihren sein. 


c) Sonstige Verfahren. Dieser Abschnitt soll einen Uberblick iiber die in der Literatur ver- 
éffentlichten Verfahren geben, die auf ahnlichen Uberlegungen und Versuchsanordnungen wie die 
bisher beschriebenen Methoden beruhen. 

1. ImpulsmeBverfahren nach Barth!. Im 1ex(4| 


Verfahren nach Abschnitt 2 konnte man die 2 ve 
Gutgeschwindigkeit c in einem zylindrischen 2w 5 

Rohr ohne weiteres errechnen, sobald der 0 5 v,=%5m/sek 
Druckverlauf iiber diesen Weg bekannt war. et eell ce le 
Voraussetzung fiir diese Rechnung war, daB 5: ELELELE 
auBer der Beschleunigungskraft alle anderen abl g 2 2 A 
Krafte vernachlassigt werden durften. Dies ‘ oP 
ist jedoch nicht immer der Fall. Die Beschleu- ee ae be 
nigung des Gutes nimmt mit zunehmender Gut- Z, [y-¥ m/sek 2 
geschwindigkeit ab und wird im Beharrungs- ® | v,-af5m/sek : 

zustand zu Null, wenn das Material seine End- | %=20 m/sek 
geschwindigkeit erreicht hat. Der Luftwider- Be 
stand des Gutes ist dann genau so groB wie ? 10 ee e 


die Kraft, mit der es an der Rohrwand reibt, Abb. 16. Gleichwertiger Kugeldurchmesser dj, in Abhangigkeit von der 
zuziiglich des Gutgewichtes. Fir diesen Fall Materialdichte ,®. Material: Kontaktstoff. 

laBt sich die Gutgeschwindigkeit nicht mehr einfach errechnen, da in dem Reibungsfaktor eine 
weitere Unbekannte aufgetreten ist. Mit dem ImpulsmeBgerat nach Barth kann die Material- 
geschwindigkeit gemessen werden. 

Das Material—Luft-Gemisch wird hierbei, aus der waagerechten Férderleitung kommend, in 
einen reibungslos aufgehangten unten offenen Zylinder geleitet und in ihm um 90° senkrecht nach 
unten abgefiihrt. Auf diese Weise wird der gesamte Impuls des Gemisches in Richtung der Férder- 
leitung auf den Zylinder iibertragen und kann mit Waage und Gewichten abgewogen werden. Eben- 
so wie der Impuls des Gemisches kann der der reinen Luft bestimmt werden. Die Differenz zwischen 
beiden ergibt den Impulsanteil des Materials. Uber den Impulssatz lat sich hieraus leicht die 
Gutgeschwindigkeit bestimmen. Zusammen mit den MeBwerten der Luftgeschwindigkeit, der ge- 
férderten Gutmenge je Zeiteinheit und des Druckes kénnen jetzt wieder die den Luftwiderstand 
kennzeichnenden Werte w und d, erhalten werden. 


2. Widerstandsbestimmung nach Blaine?. Bei diesem Verfahren wird vom Widerstand, den 
eine Materialsaule einem durchstrémenden Gas entgegensetzt, auf die Kornfeinheit geschlossen. 
Dies ist nur méglich, wenn das Zwischenraumvolumen der untersuchten Materialsdule immer das 
gleiche ist. Um dies zu erreichen, wird eine jeweils genau bestimmte Gewichtsmenge, d. h. ein 
genau bestimmtes Volumen des Materials in den zylindrischen Hohlraum des von Blaine vorge- 
schlagenen MeBgerates eingefiillt. Aus der DurchfluBmenge und dem Druckverlust des durch das 
Material strémenden Gases kann die KorngréBe bestimmt werden. Gegeniiber den bisher beschrie- 
benen Verfahren hat die Methode den Nachteil, daB nicht im freien Luftstrom gemessen wird und 
deshalb der Einflu8 der Materialschichtung nicht ausgeschaltet werden kann. 


3. Verfahren nach H. M. Carlson, P. M. Frazier und R. B. Engdahl?. Die Verfasser haben in 
eine Rohrleitung eine Diise und dahinter eine Blende eingebaut. Bei Férderung eines Gas—Material- 
Gemisches beniitzen sie die Wirkdruckanzeige an der Blende zur Luftmengenbestimmung und die 


* Diese Abbildung gibt als Beispiel den charakteristischen Verlauf wieder. Kurven fir die weiteren Luft- 
geschwindigkeiten und vollstandige Zahlentafeln siehe die Dissertation des Verfassers. 
1 W. Barth, Strémungsvorgange beim Transport von Festteilchen und Flissigkeitsteilchen in Gasen. 


Chemie-Ing.-Technik 30 (1958), S. 171. ; 
2 F, Gille, Die Priifung der Mahlfeinheit mit dem Gerat von Blaine, Schriftenreihe der Zementindustrie 


Heft 7, 1951. 
3 H. M. Carlson, P. M. Frazier u. R. B. Engdahl, Trans. Am. Soc. Mech. Engs 70 (1948) S. 65. 
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Wirkdruckanzeige an der Diise zur Materialmengenbestimmung. Da jedoch alle theoretischen 
Uberlegungen von der Verwendung der DurchfluBgleichung fir Diisen und Blenden ausgehen 
und die Einfliisse der Materialdichte, der Froudeschen Zahl usw. nur als Korrektur der Konstanten 
der DurchfluBgleichung erscheinen, kann das Verfahren nur nach Eichung der gegebenen Anlage 
verwendet werden. 

4. Verfahren nach L. Farbar!, Dieser baute zylindrische MeBdiisen in pneumatische Férder- 
leitungen ein und bestimmte den Druckabfall langs dieser Diisen. Die Messungen zeigen den be- 
kannten Verlauf des Druckabfalls von Gas—Material-Gemischen in zylindrischen Diisen. Die 
theoretischen Uberlegungen gehen von der Anwendung der Energiegleichung aus. Dies ist jedoch 
nicht zulassig, da mit Energiebetrachtungen die Materialreibung nicht erfaBt werden kann. Da- 
gegen ist die Verwendung des Impulssatzes exakt, denn bei ihm miissen nur dufBere Krafte, die 
gemessen werden kénnen, betrachtet werden. 

5. GutmengenmeBgerdt nach W. Barth?. Das Gerat ist so ausgefihrt, daB es in pneumatische 
Forderleitungen eingebaut werden kann. Die charakteristischen Teile der MeBeinrichtung sind ein 
zylindrisches und ein darauffolgendes, diffusorartig erweitertes Rohr. Durch geeignete Mab- 
nahmen kann erreicht werden, daB die Druckdifferenzen im zylindrischen Teil nur abhangig von 
der Materialmenge und die Druckdifferenzen im Diffusor nur abhangig von der Luftmenge werden. 
Aus den gemessenen Druckunterschieden zwischen Anfang und Ende des zylindrischen Rohres 
bzw. Anfang und Ende des Diffusors kann mit Hilfe einer EKichkurve die Gutmenge bestimmt 


werden. 


4, Schlu&bemerkungen. Ziel der Untersuchungen war es, eine méglichst einfache und allgemein 
gultige Darstellungsweise des Widerstandes eines Feststoff—Luft-Gemisches zu finden. Unter 
diesem Gesichtspunkt wurden die bekanntesten Verfahren der Widerstandsbestimmung betrachtet. 
Fiir manch wertvollen Hinweis zu dieser Arbeit méchte der Verfasser Herrn Prof. Dr.-Ing. W. Barth 
aufrichtig danken. 


(Eingegangen am 22. November 1957.) 
Anschrift des Verfassers: Dipl.-Ing. Walter Rausch, Hof/Saale, Kulmbacherstr. 29. 


' L. Farbar, Industrial and Engineering Chemistry 44 
> y 44 (1952) S. 2947. 
 W. Barth, R. Nagel u. K. van Waveren, Chemie-Ingenieur-Technik 29 (1957) S. 599. 
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Uber eine zweckmafige Walzenform fiir das Auswalzen diinner Folien 
Herrn Prof. Dr. Rudolf Sonntag zum 65. Geburtstag gewidmet 
Von W. Roth 


1. Einleitung. Zur Auswalzung diinner Folien werden Walzen verwendet, deren Kontur 
nicht zylindrisch ist. In der vorliegenden Arbeit wird ein Verfahren angegeben, welches es ge- 
stattet, die erforderliche Kontur bis zu einer gewuinschten Genauigkeit zu ermitteln. Wenn z. B. 
Aluminiumfolien mit einer Starke von 0,02 mm ausgewalzt werden, dann bedingt ein so geringes 
MaB eine auferordentlich sorgfaltige Ausfiihrung der Walzen. Die Walzen miissen nicht nur auf 
Hochglanz poliert sein (Rauhigkeit etwa 0,2 4), sondern die Kontur mu8 auch nach einem ge- 
nauen Gesetz von der Zylinderform abweichen. Waren die Walzen namlich genau zylindrisch, 
so wirde infolge der durch die Krafte P (Abb. 1) bedingten Durchbiegungen die Starke der Folie 


P 


Abb. 3. Belastete Walze mit faBformiger Kontur. 


Abb, 1. Zylindrische Walzen in deformiertem Zustand. 


Abb. 2. Walze mit faBformiger Kontur. Abb. 4. Belastungen an der undeformierten Walze. 


zur Mitte hin zunehmen. Die Walzen wiirden sich wesentlich an den Randern abstiitzen, was 
einen relativ starken Fehler zur Folge hatte. Um eine solche ungleichmaBige Folienstarke zu 
vermeiden, miissen die Walzen faBformig ausgefiihrt werden (Abb. 2). Die genaue Kontur der 
Walze ergibt sich dabei aus der Bedingung, dafB nach Aufbringung der Krafte P auf die Zapfen 
sich eine konstante Pressung p in der 4uBeren Druckfaser der Walze einstellen mu8, wobei noch 
gefordert wird, daB diese Druckfaser in eine gerade Linie tihergeht (Abb. 3). 


2. Lésung des Problems. Eine Lésung der Aufgabe wird iiber die Differentialgleichung der 
elastischen Linie gewonnen. In Abb. 4 sind an der undeformierten Walze die wirkenden Krafte 
eingezeichnet. Das Biegemoment an der Stelle x ist 


M(x) = P(a+2)—py, 


wobei p = P/I gesetzt werden kann, wie die Gleich- 
gewichtsbedingung gegen Verschiebung in verti- 
kaler Richtung aussagt. Die Differentialgleichung 
der elastischen Linie ist somit 


2 
x 
Ee t= 21 Abb. 5. ZweckmaBige Wahl des Koordinatensystems zur 
(1) Beschreibung der Kontur. 


Nee ag 


Bei der Integration dieser Gleichung muf nun darauf geachtet werden, daB das Tragheitsmo- 
ment J nicht konstant ist. Wegen der jedoch nur geringen Abweichung von der Zylinderform 
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darf in z (Abb. 5) linearisiert werden, so daf 
1 1 Z\—4 1 3 
Sy a —\ y—/l—4— 
Ts el aa 
ist, mit Jy = r4/4 (Abb. 5). Gleichung (1) hat damit die Form 
Hine Re Zz x? 
y =#3,(! 1 Fy (a x zi): (la) 


Soll die Druckfaser der Walze in eine gerade Linie titbergehen, so muB (Abb. 5) die geometrische 
Bedingung 


Sf = 


oder 


erfillt sein. Da die Gleichung der elastischen Linie selbst weniger interessiert, wird y’’ aus (la) 
eleminiert, und es liegt so die Differentialgleichung fiir die Kontur der Walze 


Path ( 4 -)(« bx zi) = ° (2) 


| 
gine ph 


Yo 
vor. Aus Symmetriegriinden empfiehlt sich die Substitution 
x—l=u, 

und aus (2) erhalt man fiir z = z(u) die lineare Differentialgleichung 

dz 4 P oly aca oe eee (2a) 

du2 Ein" Z| zi EJo ! z)— 3° e 

Die allgemeine Lésung einer linearen Differentialgleichung zweiter Ordnung muB zwei Inte- 

grationskonstanten enthalten, welche durch die Randbedingungen bestimmt werden. Bei dem 
vorliegenden Problem muf dz/du = 0 an der Stelle u = 0 und z = 0 an den Stellen u = + 1 


sein. Da die gesuchte Funktion z eine gerade Funktion beziiglich der Stelle u = 0 sein muB, 
wird der Ansatz 


ye O2n we (3) 


n=0 


auf die Lésung von Gleichung (2a) fiihren. Die Randbedingung dz/du = 0 bei u = 0 ist mit 


beliebigen x2, (n = 0, 1,...) erfiillt; aus diesem Grund wird also in (3) nur noch eine Integra- 
tionskonstante enthalten sein. Durch Einsetzen von (3) in (2a) findet man fiir die «2, das System 
ie lege? PL domes 
a PR Papin eo" 
ae Os dA ek i is hs 
POE a ret ET re OEE OT. 
2 Pb 2 P (4) 
Og yc tea ee 
Pl 2 le 2 P 
by = Wis 


ots ETF, e Bo Utes 


wobei 2a +1 = b gesetzt ist. Die Koeffizienten 2, lassen sich also rekursiv in Abhangigkeit 
von a, berechnen, und die Konstante a, tibernimmt demnach die Rolle der noch zu bestimmen- 
den Integrationskonstanten, welche dann durch die Randbedingung 


co 


(a oe O2n Le (5) 


n=0 


festgelegt ist. Die maximale Abweichung z,,,, von der Zylinderform liefert Gleichung (3) mit u = 0 


max 
Cees == OE OG. 0 


Die Rechnung gestaltet sich tibersichtlich, wenn die Substitution 
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durchgefiihrt wird. Es ist f der Biegepfeil eines eingespannten Tragers von der Linge I, der 
Biegesteifigkeit E Jj, belastet durch die Endkraft P und die Streckenlast p = P/l (Abb. 6). Die 


Gleichungen (4) kénnen nun auch in der Form 


48 f b 12 b 
+) PA See SS aage 
cy aeraie Has) aie 
48 f b 48 f 12 
 J4 : 2 q — __ 
3-40, 1 5 To! 5 7 00 ie 
4 
BY Gaelic seo, E +2L,,2=0, (4a) 
0 
7-80,2 SL 5,684 22 o,h=0, 
0 


angegeben werden. 
Eine praktische Auswertung wird sich nun so gestalten, daB man zunachst einmal annimmt, 
die Kontur der Walze werde durch die Parabel 
2 =X + Ow? 
hinreichend genau beschrieben, also alle a2, (n = 2,3,...) waren Null. Die Erste der Glei- 
chungen (4a) liefert dann 


Re 
Xo IE = oT (430 SS5 f) 
und (5) damit 
6 6b 
epee, ie ae 2 sit yet : 
a ee an 24 6 be , Abb. 6. Zur Veranschaulichung 
cs ae ir, der gewahlten Substitution f. 


Betrachtet man diesen Wert e, als einen ersten Naherungswert, so wird man zu einer zweiten 
Naherung kommen, wenn von Ansatz (3) 


By = Ogg + Ooo U + Oy, ut (6) 


als Approximationsfunktion unterteilt wird. Die Zweite der Gleichungen (4a) liefert 


oy, 4 = 7 FF yt | (4 f be f) 


yao 
WO OAT bef | 
tere 


und aus (5) folgt 


b 
l 

Zmax2 = 2 = %o2 = pe ee 
ji + : (6 


SchlieBlich sei noch fiir den Ansatz 
%3 = os + Og3 U® + ys Ut + Orgs US (6a) 


der nachste Term aus der Dritten der Gleichungen (4a) 


8 b f [24/b\2f Ae 
teal =e 7 [5 (T) (45590 f) 


angegeben, mit dem der Naherungswert 


sitar rye le) Gl 


Zmax3 — 3 = %o3 — 4 ia f 2m) Ti i f\, 7168 /b \8 hi 3 
f er: (6 I 7 Ts T 795 7 (15 I (5) as (7 al 


erreicht wird. ; 
Es liegt in der Art der Problemstellung und der hier angewandten Lésungsmethode begriindet, 
daB die Folge («om — e) eine Nullfolge ist, also lim oo m = e gilt. Die Folge der wom konvergiert 


m— oo 
sehr rasch gegen den gesuchten Grenzwert e, wie eine Betrachtung der oben angegebenen Werte 
e,, e, und e, zeigt. Fat man die e; (t = 1, 2,.. .) als Funktionen des Quotienten f/ry auf, der als 
23 
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auRerordentlich klein anzusehen ist, so stimmen e, und e; in den konstanten Gliedern bereits 
iiberein. Obwohl eine Beriicksichtigung selbst noch des in f/r, linearen Gliedes tiberfliissig scheint, 
soll die linearisierte Form von e, angegeben werden, welche dann als eine den praktischen An- 
forderungen ausreichende exakte Richtformel bezeichnet werden kann. Es ist also 


ers (6 1 ! ‘A Be |S (7) 


Aber auch ein anderer Gedankengang la$t die schnelle Konvergenz der «om erkennen und 
beweist damit die Brauchbarkeit der entwickelten Methode. Eine Naherungslésung fiir die ge- 
suchte Kontur der Walze liefert namlich schon die elementarste Biegelehre. Abb. 7 zeigt einen 
eingespannten Trager von der Lange /, belastet durch die Endkraft P, das Endmoment P a und 
die Streckenlast p = P/l. Der Querschnitt des Tragers sei konstant und damit auch die Biege- 
steifigkeit E J = E Jy. Unter Einflu8 der Belastungen wird der Trager sich durchbiegen. Der 
Endquerschnitt senkt sich dabei um die Strecke 


b IPI 1 b 
joe (6 I ET = 5 (0 i uf. 
Dieser Biegepfeil kann als Naherungswert fiir e angesehen werden, so 
da8 dann aber auch fiir die Kontur die Funktion 


a= fgo—y 
Abb. 7. Erste Niherung durch Ersatz- als Naherung gilt, wobei y die Biegelinie des Tragers (Abb. 7) ist. 
triger mit Binzellast, Streckenlast 1). @ eine GréBenordnung von nur einigen hundertstel Millimetern 
und Endmoment. $ % 
besitzt, praktisch eine zylindrische Walze vorliegt, ist die Annahme, 
mit konstantem Querschnitt zu rechnen, zutreffend. Man kénnte nun fortfahren und von der so 
gefundenen ersten Kontur mit einem zweiten Iterationsschritt eine noch genauere Funktion 
fiir die Walzenkontur herleiten, doch mite in diesem zweiten Rechnungsgang bereits die 
Veranderlichkeit des Traigheitsmomentes J beriicksichtigt werden. 


Allgemeiner ist es iiber die Differentialgleichung fiir die Kontur der Walze 
d?z M 
ae ey (8) 


mit 


fiir kleine z/ry méglich, ein Iterationsverfahren anzugeben, welches direkt von der Gleichung (8) 
ausgeht und die Form hat 


a M 
dx? E ee 
oder 
x2 
dz, P(a+2—3) l 7p! 
FES EJ, rami (9) 
wobei 


gesetzt ist und unter M das schon oben angegebene innere Biegungsmoment an der Stelle x der 
Walze zu verstehen ist. Die z, (n = 1,2,...) geniigen den Randbedingungen z,(0) = 0 und 
z,(!) = 0. Der erste Schritt der Iteration wird mit z, = 0 ausgefiihrt, so daB (9) die Form hat 


P( bz 
Ee a+e— 3) 
dx? 1M 


Man erhalt den Wert 


Pe ee eee (10) 
2E Jo Tals ag a 
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welcher mit x = 1 auf den Naherungswert 


Pe b 
¢® So aT ET, (6 I i 
fiihrt, der schon angegeben wurde. 


Der zweiten Iteration liegt die Gleichung 


dx* ET, 


ds, Pletal 


zu Grunde mit z, aus (10); doch ist die Ermittlung von z, schon so umstandlich, daB hier die Rech- 
nung bereits abgebrochen werden soll. Der Rechnungsgang verliert an Ubersichtlichkeit. 

Wenn gréBere Genauigkeit gefordert wird, ist das oben dargelegte Verfahren vorzuziehen. 
Kin Vergleich der ao m mit fy zeigt, daB in a». und x3 die konstanten Glieder bereits mit f, iiber- 
einstimmen, und dies ist ein Anhaltspunkt fiir die Genauigkeit der abgeleiteten Formel (7) fiir e. 

Es steht im Ermessen des Rechners, das Verfahren noch so zu verallgemeinern, daB der Aus- 
druck fiir das Tragheitsmoment 1/J um weitere Potenzen von z erweitert, so daB also 


1 1 Zz 2 
5=3% (1 12-465 wes 


gesetzt wird. Dies ware zwar nur fiir Walzen mit starkerer Kriimmung in der Meridianebene 
erforderlich, eine Erweiterung der Theorie also, welche sich fiir das vorgelegte zu lésende Problem 
nicht lohnt. Eine Folge davon ware auch, da die Gleichungen (4) bzw. (4a) dann komplizier- 
teren Bildungsgesetzen geniigen. 


3. Zahlenbeispiel. Ein Zahlenbeispiel mag die Verhaltnisse noch erhellen. Es sei 1 = 45 cm, 
¢—I2cm, ry, = 14cm, E = 2,2-10% kg/cm? und P = 3000kg. Mit J, = 3,02- 104 cm‘ ist 
dann f = 8,6-10-? mm und f, = 1,41: 10-2? mm. Fir den Quotienten f/r) erhalt man f/ry = 
= 6,14- 10-5, und Formel (7) liefert e = 1,409- 10-2 mm. Der Unterschied von e und fj ist also 
bereits so gering, daB er mit Rechenschiebergenauigkeit gerade noch feststellbar ist. Die Gleichung 
der Kontur ware dann nach (6) oder (6a) hinreichend genau beschrieben. 


4, Zusammenfassung. Es werden im Wesentlichen zwei Verfahren zur Ermittlung der 
Walzenkontur angegeben. Wahrend das erste Verfahren darauf beruht, daB von vornherein in 
der Differentialgleichung der elastischen Linie der Walze eine Veranderlichkeit des Tragheits- 
momentes J beriicksichtigt wird, entwickelt sich der Rechnungsgang iibersichtlich und so, da 
praktisch bis zu jeder gewiinschten Genauigkeit die notwendige Abweichung von der zylindri- 
schen Form der Walze ermittelt werden kann. Das zweite Verfahren verliert mit den ersten 
Iterationsschritten bereits an Ubersichtlichkeit und wird deshalb nicht weiter entwickelt. Es 
werden jedoch Vergleiche gezogen, die beweisen, dafi beide Verfahren sehr rasch konvergieren. 


Aus dem Lehrstuhl und Institut fiir Technische Mechanik (Prof. Dr.-Ing. Rudolf Sonntag) Karlsruhe. 


(Eingegangen am 22. November 1957.) 
Anschrift des Verfassers: Dipl.-Ing. Werner Roth, Karlsruhe, Kaiserstr. 12, Institut f. Techn. Mechanik. 
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Beitrag zur Theorie ebener Fachwerke 
Von H. Baldauf 


1. Einleitung. Die iibliche Theorie zur Berechnung der Stabkrafte in Fachwerken geht u. a. 
von der Voraussetzung aus, daB die Stabe in den Knoten durch reibungsfreie Gelenke mit einan- 
der verbunden sind, und da8 die Lasteintragung nur in den Knoten erfolgt. Unter dieser Voraus- 
setzung erhalten die Fachwerkstabe nur Normalkrafte. Dabei wird vernachlassigt, daB die Stabe 
in den Knoten biegesteif miteinander verbunden werden und durch ihr Eigengewicht Belastungen 
senkrecht zur Stabachse erhalten. Dadurch treten in Fachwerkstaben stets Biegemomente auf, 
die im allgemeinen zu kleinen Zusatzspannungen fiihren und als Nebenspannungen bezeichnet 
werden. Die exakte Beriicksichtigung der Dehnsteifigkeit und der Biegesteifigkeit von Fach- 
werkstaben fiihrt bei Anwendung der Kraftgré8enmethode auf einen hohen Grad statischer Un- 
bestimmtheit, so daB man sich zur Berechnung der Nebenspannungen meist mit Naherungs- 
lésungen begniigt!. Kennzeichnend ist dabei, daB zunachst die Berechnung in der tiblichen Form 
(Cremona-Plan, Rittersches Schnittverfahren) durchgefiihrt wird, und daB die Nebenspannungen 
in einem zweiten Rechnungsgang bestimmt werden. Bei dem zweiten Rechnungsgang wird meist 
nur die Biegesteifigkeit der Stabe beriicksichtigt. In der vorliegenden Arbeit wird ein Verfahren 
beschrieben, das in Anlehnung an das Reduktionsverfahren? gestattet, die Schnittkrafte eines 
Fachwerkes unter Beriicksichtigung der Dehnsteifigkeit und der Biegesteifigkeit in einem Rech- 
nungsgang zu bestimmen. 


Abb. 1. Definition der Belastung und der Randschnittkrafte Abb. 2. Definition der Randverformungen 
eines Fachwerkstabes. eines Fachwerkstabes. 


2. Grundlegende Beziehungen am Stab und an den Knoten. Es werden folgende Voraus- 
setzungen getroffen: Die Theorie erster Ordnung ist anwendbar. Das Hookesche Gesetz ist giiltig. 
Die Abmessungen der Knotenscheiben sind klein gegeniiber den Stablangen, so daB die Stab- 
langen gleich den Netzlangen des Fachwerksystems gesetzt werden kénnen. Das Tragheitsmo- 
ment J;,, die Flache F;, und das Kigengewicht g,, eines Stabes ist zwischen zwei Knoten i und k 
konstant, ebenso der Elastizitatsmodul E. Zwischen den Schnittkraften (Abb. 1) und Form- 
eae (Abb. 2) an den Stabenden i und k bestehen dann folgende Beziehungen (s. auch 


y a ee a ee eee , oe 
R HES FY  ES;, 6 EJ,, 24 me: 
1 ers is 13 
ees =) ee es Sik lik 
Pk 0 ae bse li. Ev 2 0 0 ESS Gu Pi 
2 
M, 0 0 =k aly 0 0 zs Sik ze M; 
aa 1 
Q:, 0 0 0 ap 0 0 ae li, Q; q 
uw 0 0 0 =] eee eee 
hk . LEE, Po re ms 
Ni, 0 0 0 0 0 pl — #34, by N; 
1 0 0 0 0 0 0 cig 1 


19) Hartmann, Stahlbriicken (Der Briickenbau 3. Band), Wien 1951, S. 6 
2 §. Falk, Ing.-Arch. 24 (1956) S. 216. and) ea at ool 


XXVI. Band 1958 H. Baldauf: Beitrag zur Theorie ebener Fachwerke 339 


oder kurz 

Me = Bin Ni « (La) 
Die von der tiblichen Form abweichende Definition der Formanderungen an den Stabenden ist 
im Hinblick auf die Transformation der geometrischen Randbedingungen an den Knoten (10) 
zweckmaBig. AuBerdem ist bei dieser Definition vorteilhaft, daB bei der Umkehrung der Be- 
rechnung 

= Lei Me (1b) 


die Matrix {,; mit Ausnahme der von der Belastung g;, abhangigen Glieder mit der Matrix &,, 
iibereinstimmt: 


i es Bi 5, = Bik 8 Mik 
i= Ui, fir gy, =. (le) Bi, = Sik SID Hj, 

: ne ee . < Lage a at 
Die von der Belastung g,, abhangigen Glieder andern Bri — Shi? *bi Fk 
. . . a - gl =— 
bei der Umkehrung der Berechnung ihr Vorzeichen Bi hi “hi Bik 

fo) to) 


(Abb. 3). 

Fiir die Zahlenrechnung empfiehlt es sich, mit 
dimensionslosen Zahlen zu rechnen. Auf Grund der 
folgenden Beziehungen (2) wird die Matrix (1) um- 
gerechnet. 


JS Ski 
7 


+ Pele — mt tele a5 
ray 13 Bie Or 1? Re pd u=—U ea) ? 
M=M*P.1, Q=Q*P,, N=N*P,, $(2) 

Ie OK; 
x fe 


(he Sie EG Abb. 3. Zerlegung des Eigengewichtes in die Lastkompo- 
le nenten parallel und senkrecht zur Stabachse. 


Dabei bedeuten P,, ], und J, beliebig wahlbare Kraft, Lange und Tragheitsmoment. Zweck- 
maBig werden sie so ausgewahlt, daB die Verhaltnisse J,/J;,, und 1;,/1, méglichst gleich 1 werden 


oder nahe bei | liegen. Wird (2) in (1) eingesetzt, so ergibt sich 


a err) teste) 9 Oe  \F 
OF ea ee a 
M 0 0 = tt 0 0 + iat, || my 
Oe} too 0 pf 0 0 +62, lor | 
uj 0 0 0 0 1 Fh = wt, || ut 
N; 0--.. 0 0 0 (ciate ae OL ve 
1 0 0 0 0 0 0 ay 1 
Dabei bedeuten (Abb. 3) 
ee od Ee aa | 
m= +8). aie) (4) 
eerie Nees (iz): | 


Sind an einem Stab die geometrischen Randbedingungen an den Enden 1 und k vorgeschrieben, 
so sind damit die Randschnittkrafte nach (3) bekannt. Die Matrix zerfallt in zwei Teile 


y 

7 c l; 2 1 Als l 2 i: i 

ae li site ae (7) v6 Su (7) + Ye \| 9% 
(")- ; pie Ms (5) 

“ Fs (lit 1 Je (lik)? = i 

je eee Tue —zIule) +) er 
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und A 
a Je [lik ~ x es (6) 
ve ( : F,B(t) mag Wah 
if 


Daraus folgen die Randkrafte im Punkt 1 


lik 
M* = es Sila] > 
y Ik Ske ENE Ip ental\ Uf 
zIult) 
lik 
* MLN pee 3 * 
? 20+ Ye) a (EF ve) Haat (7) 
Us 1-Je (lix\® Daal hae 
6 Jik (i) 
Aes uf tue | Sik [lik 
N; Sha ft Bie oa Tosh 
Fal.) 


Durch Vertauschen der Indices i und k erhalt man die entsprechenden Ausdriicke fiir die Schnitt- 
krafte im Punkte k. 


PPR 


My=My Pl 


we ken=P Keo! 


Abb. 5, Beziehungen zwischen den Verschie- 


Abb. 4, Belastung und Randschnittkrafte 
bungen zweier Stabenden an einem Knoten. 


am Knoten K. 


Die Gleichgewichtsbedingungen an einem Knoten K lauten (Abb. 4) 
Life “if 2 N,, C08 &, — D, Q,, Sin, aa 0, 
K K 


Py, + DN, sina, + FQ, cosa, =0, (8) 
K K 
== x oD DES 0 
K 


Sind an einem Knoten alle Stabendkrafte mit Ausnahme der Krafte am Stabende m bekannt, 
so ergeben sich diese Randschnittkrafte aus (8): 


OF + cos (x, —a,,) + sin (a, —«,,) 0\ /Q* 
Nx | =— > — sin (x, —«,,) + cos («, —<@,,) Orne 
Ms s 0 0 +1) \ms 
=“ COS sini 0 PRE 
—2 | silo cos xe 0 Pele (9) 


0 0 —1} \ Mg 
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Die senkrechte und waagerechte Verschiebung und die Verdrehung eines Knotens ist bestim- 
mend fir die geometrischen Randbedingungen aller Stabe, die in diesem Knoten angeschlossen 
sind. An einem Knoten bestehen zwischen den geometrischen Randbedingungen der Stab- 
endenn und m folgende Beziehungen (Abb. 5): 


‘y* + cos (a, — On) + sin («,, — &,,) a ee 
un | =| — sin (a, — Oh.) + cos («, — &,) 0 Un ey) 
pt 0 0 4 1 Pn 


3. Rechnungsgang. Bei einem zweistabigen Knoten (Abb. 6 Knoten A) wird an einem 


Doe (Stabende 0) der Vektor M mit den unbekannten Komponenten Ve, sed Oa. 


o> Vo angenommen: 
ye \ l 0 0 0 0 ( 0 0 
ge 0 1 0 0 0 0 0 
Ms 0 0 1 0 0 0 0 
to | OF {=| 01 ¥e +101 %e +).0| Mo] 1100 +] oles +|ojNi-+]o]1.. (1) 
us 0 0 0 0 il 0 0 
Ne 0 0 0 0 0 1 0 
i 0 ) 0 0 0 ) 0 I! 
Abb. 8. Streben- und Pfostenfachwerke. 
Abb. 7. Reihenfolge des Rechnungsganges. Abb. 9. K-Fachwerk und Rautenfachwerk. 


Mit Hilfe von (3) wird der Vektor 7, in Abhangigkeit der unbekannten Komponenten berech- 
net!. Auf Grund der Gleichgewichtsbedingung (9) und der Kontinuitatsbedingung (10) wird 
am Knoten A der Vektor 7, bestimmt und darauf mit (3) der Vektor 73. Mit der Kenntnis 
von 73 und 7, liegen iiber (10) die geometrischen Randbedingungen fiir den Stab 45 fest, so daB 
mit (7) die Randschnittkrafte am Stabende 4 berechnet werden kénnen und anschlieSend mit 
(3) die Schnittkrafte am Stabende 5. Am Knoten B kann nunmehr auf Grund der Gleich- 
gewichtsbedingung (9) und Kontinuitatsbedingung (10) mit 7; und 4, der Vektor 7, bestimmt 
werden, anschlieBend mit (3) der Vektor 7,. Damit liegen die geometrischen Randbedingungen 
des Knotens D fest, so daB sich nunmehr die Schnittkrafte des Stabes 89 nach (7) berechnen 
lassen. Dieser Rechnungsgang wird sinngemaB fortgesetzt. In Abb. 7 ist die Reihenfolge der 
Rechnungen angegeben. SchlieBlich ergeben sich die Vektoren 71) und 7.,, die die Bedingungen 
(9) und (10) erfiillen miissen. Dadurch erhalt man 6 Bedingungsgleichungen fiir die 6 unbe- 


1 Siehe FuBnote 2 von S. 338. 
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kannten Anfangswerte. Sind die Stiitzkrafte unbekannt (statisch unbestimmt gestitzte Fach- 
werke), so treten noch t unbekannte Stiitzkrafte hinzu, dehen aber t geometrische Stiitzen- 
bedingungen gegeniiber stehen. Insgesamt sind also zur Berechnung der 6 + t unbekannten 


Komponenten 6 + t Bedingungsgleichungen verfiigbar. 
Im betrachteten Beispiel (Abb. 6) sind die folgenden Stiitzenbedingungen vorgeschrieben: 


yo = 90, uy = 0, Yip =O. 


Dadurch enthalt der Anfangsvektor nur noch 4 unbekannte Komponenten. Werden die Stiitz- 
krifte als bekannt vorausgesetzt (statisch bestimmte Stiitzung), so sind von den 6 Bedingungen 
(9) und (10) am KnotenG nur noch 3 voneinander unabhangig. Die vierte unabhangige Be- 
dingung folgt aus der Randbedingung y*, = 0. Bei statisch unbestimmter Stiitazung gehen die 
Stiitzkrafte in den Gleichgewichtsbedingungen (9) als unbekannte GréBen ein und sind in der 
weiteren Rechnung entsprechend mitzufihren. 


4, Grenzen der Anwendharkeit, Erweiterungsméglichkeiten. Das angegebene Verfahren kann 
auf Strebenfachwerke mit und ohne Hilfspfosten und auf Pfostenfachwerke (Abb. 8) angewendet 
werden. Dabei beginnt die Berechnung an einem zweistabigen Knoten (A bzw. G). Die An- 
wendung bei K-Fachwerken und Rautenfachwerken (Abb. 9) ist in der beschriebenen Weise 
nicht méglich. 

Die Beschrankung auf konstante Flachen- und Tragheitsmomente und konstante Gewichte 
der einzelnen Stabe ist nicht nétig. Die Verdnderlichkeit dieser GréBen kann ohne weiteres 
beriicksichtigt werden, nur wird die Transformationsmatrix (1) komplizierter, auch ist dann 
die Beziehung (1c) nicht mehr erfillt. Ebenfalls kann das Verfahren auf Lastangriff von Einzel- 
kraften zwischen den Knoten erweitert werden. 

Soll auch der Einflu8 der Langskraft auf die Formanderungen senkrecht zur Stabachse 
beriicksichtigt werden, so muB von der Differentialgleichung 

d*y f 


d4 
E Jy, 57 — Nun ga = Bie (12) 


ausgegangen werden, mit den Lésungen (fiir J;, = konst., N,;, = konst., g’, = konst.) 


tk 


y=A+Bx+ Coosy,x« + Dsiny,, « — Bik 2 fu Nee 


2 Nik 
= ii 7. Pa ea Sik 2 ws 
yea Bx + COo) »;, « + D Gin v,, x aN fir NG 05 (13) 
mit 
oF pean Nik a eee), Nik 
Mik Ein’ Te Oe iee 


Auf diese Weise wire es méglich, die Knicksicherheit eines Fachwerkes fiir Knicken in der 
Fachwerkebene festzustellen und die Traglast eines Fachwerkes zu berechnen. Allerdings muB 
hierbei zunachst die Langskraft in den einzelnen Staben geschatzt werden (Cremonaplan, Ritter- 
scher Schnitt). 


5. Zusammenfassung. Es wurde ein Verfahren erlautert, mit dessen Hilfe die exakte Be- 
rechnung der Langskrafte, Querkrafte und Momente spezieller Fachwerkssysteme in einem 
Rechnungsgang durchgefiihrt werden kann. Auf mégliche Erweiterung zur Traglastberech- 
nung und Stabilitatsberechnung wurde hingewiesen. 


(Eingegangen am 26. November 1957.) 
Anschrift des Verfassers: Dipl.-Ing. Heinrich Baldauf, Dresden A 27, HalbkreisstraBe 71. 
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Ein Beitrag zur Berechnung regelmafig gegliederter Scheiben * 
Von H. Beck 


1. Einleitung. Gegliederte Scheiben stellen hochgradig statisch unbestimmte ebene Rahmen- 
tragwerke von regelmaBiger Struktur dar. Ihre Berechnung erfordert im Regelfall nach den 
tiblichen Methoden der Rahmenstatik (KraftgréBenmethode, Formanderungsgréienmethode) 
einen erheblichen numerischen Aufwand, und zwar auch dann, wenn man einen iterativen 
Lésungsweg benutzt. 

Der Grundgedanke des im Folgenden abgeleiteten Naherungsverfahrens besteht darin, viele 
unbekannte statisch iiberzdhlige EinzelgréBen zu méglichst wenigen unbekannten Funktionen 
zusammenzufassen. Anstelle eines linearen Gleichungssystems zur Bestimmung der vielen 
EinzelgréBen tritt dann ein Differential- a 
gleichungssystem zur Bestimmung der P(e) 
wenigen, die EinzelgréBen zusammen- reeeeeell TH 
fassenden Funktionen. Im giinstigsten ! 

Fall — der sich bei diesem Verfahren tat- l 
sachlich einstellt — reduziert es sich auf 
eine einzige Differentialgleichung zur Be- 
stimmung einer alle EKinzelgréBen festlegen- 
den Funktion. Dieser Gedanke wird durch 
den regelmaBigen Aufbau der hier zu be- 


Riegel 1 (I,4) 


i GeschoB 1 


handelnden Systeme nahegelegt, der zur 

Folge hat, daB bei Wahl gleichartiger sta- Riegel? (I) 
tisch Uberzahliger diese in Abhangigkeit 

des Ortes durch eine mathematisch ,,ver- 

ninftige“‘ Funktion erfaBbar sind. Die Er- Riegel m (I,m) 


gebnisse eines darauf basierenden Rechen- 
verfahrens sind um so besser, je gréBer die 
Zahl der durch eine Funktion zu ersetzen- 
den Uberzahligen ist, d.h. also, je um- 
fangreicher die iibliche Rahmenrechnung 
wiirde. 


; ne GeschoB m 


Riegel m+7 (I, met) 


Riege!r-1(L, n- 
2. Systemstruktur einer gegliederten gelr-t (Ens) 


Scheibe. Ein Beispiel fiir eine gegliederte 
Scheibe ist mit den gewahlten Bezeich- 
nungen in Abb. 1 dargestellt. Die Riegel 


Th 4 GeschoB r-1 


sind untereinander parallel und gleichlang, Riegelr (I>) 
und sie sind in jedem GeschoB mit der q 

gleichen Anzahl gleichabstandiger, zu ihnen D=Lp 
senkrechter Pfosten verbunden. Das Riegel- Abb. 1. Systembezeichnungen. 


tragheitsmoment ist fiir jeden Riegel iiber 

x konstant; die einzelnen Riegel kénnen jedoch ein voneinander verschiedenes Tragheitsmoment 
besitzen. Die Tragheitsmomente der Innenpfosten sind je GeschoB gleich, kénnen aber von 
GeschoB zu GeschoB verschieden sein. Ferner sei das Verhaltnis der Tragheitsmomente von 
Innen- und Endpfosten fiir alle Gescthosse gleich. Der Systemlinienabstand der Riegel, die 
lichte Weite zwischen den Riegeln und die relative Lage der Riegelsystemlinien zum lichten 
Riegelabstand kénnen von GeschoB zu Geschof variieren. Wir fiihren eine Numerierung der 
Riegel und der Geschosse ein, die wir gemaB der Darstellung in Abb. 1 von oben nach unten 
zihlen. Die Gesamtzahl der Riegel nennen wir r. Ein System mit r Riegeln hat stets r— 1 
Geschosse. Die in Abb. 1 und den folgenden Abbildungen gewahlten Darstellungen fiir die 
Lagerungsbedingungen bedeuten fiir die Ableitung des Verfahrens keine Spezialisierung. 


* Auszug Habilitationsschrift Darmstadt 1957; Referent: Professor Dr.-Ing. A. Mehmel, Korreferent: 
Professor Dr. C. Schmieden. 
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3. Struktur des Ersatzsystems. Im Sinne der Mechanik basiert das in seinen Grundgedanken 
bereits geschilderte Verfahren auf der Annahme eines kontinuierlichen Ersatzsystems fiir das 
gegebene diskontinuierliche System. Da die in dem gegebenen diskontinuierlichen System 
interessierenden mechanischen Gréfen erst durch die dem System eigene Diskontinuitat ent- 
stehen, darf das kontinuierliche Ersatzsystem die wesent- 
lichen Diskontinuitatseigenschaften des vorliegenden Sy- 
stems nicht verletzen. Die im Abstand a, vorhandenen 
Pfosten mit dem Tragheitsmoment I;,,, denken wir uns 
- : : durch kontinuierlich angeordnete Lamellen der Breite dx 

| fea, mit dem Tragheitsmoment I, ,, dx/a, ersetzt, die sich un- 
| mittelbar gegenseitig nicht an der Verformung hindern. 
em Die Tragheitsmomente J;,,, der Innenpfosten sind durch 

8 die Lamellentragheitsmomente vollstandig erfaBt. Falls die 
Tragheitsmomente der Endpfosten eines Systems von I, ,,/2 
verschieden sind (d.h. falls If, == I2,m/2 und bzw. oder 
T3 in += Io, m/2 ist), fiihren wir im Ersatzsystem noch zusatzliche Ersatzendpfosten mit den 
Tragheitsmomenten 


A A Loe. B B eae 
Ton = n(1 9 ): Tn = Ten(1 ype"); ismsr—l). (1) 


aye 


Abb. 2. Systembezeichnungen fiir nicht am 
Riegelende angeordnete AuBenpfosten. 


2,m 
ein. Wir setzen bei den weiteren Betrachtungen voraus, daB I4,, => 0 und IZ, > 0 ist. 


Wir haben in Abb. 1 zugrunde gelegt, daB die AuBenpfosten an den Riegelenden angeordnet 
sind. Fur den praktisch kaum vorkommenden Fall, da die AuBenpfosten vom Riegelende 


Py 
Riegel 1 
ae 
7 Lamellen (Inq -dx/atg) 
2 
Riegel 2 
Pr 
Riegel m 
8 
| | Liew 
Lamellen(Iz m,a2/a2) 
Pm+1 
= : = |= : Riegel m+1 


Riegel 7-1 


8 
' TEp-4 


Lamellen (Lp p.4°@2/a,) 


Riegelr 


Cs | 
D=L, L=Zp 


Abb. 3. Ersatzsystem mit Ersatzendpfosten. 


einen abetand aZ/2 und a®/2 haben (Abb. 2), ist den AuBenpfosten ein Einzugsgebiet von 
(a. + a#)/2 baw. (a, + a®)/2 zuzuweisen, so daB sich die Formeln (1) modifizieren zu 


A = 1 ad Ts im A 
Linn D) ( make |B 


B sa 1 aB UES os B 
Tis (1 aap |e 


(mae (1’) 
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’ Das Verhaltnis zwischen den Tragheitsmomenten der Endpfosten und der Innenpfosten 
wird als in allen Geschossen gleich vorausgesetzt : 


hig is i ny 
Moe a4 Teter limasr—l; — 
Ty in LE EY iW ee: : nites of Satter “i . 
Dann ist auch das Verhaltnis 
iG His 
n ,m 3 E,m 
eee bzw. Pas eek (lsm<r—l) (3) 


in allen Geschossen gleich. 

; In Abb. 3 ist das dem in Abb. 1 gegebenen System zugehérige Ersatzsystem dargestellt. 
Die auftretenden Lasten fassen wir in eine senkrecht zu den Riegelachsen wirkende Lastsumme 
p(x) und in r in Richtung der Riegelachsen wirkende Einzelkrafte H,, zusammen (Abb. 1). 


ek 
/ 


i - Riege/ m = 4+ fF} 


| Verformung : 
- ra : Riegel m+1 


a Pee 
CL) raumfeste Punkte L | ar eS eee ; in iste 
\ \| + \ 


Abb. 4. Verformung einer Lamelle. 


Fiir die rechnerische Belastung des Ersatzsystems teilen wir die gegebene Lastsumme p(x) 
entsprechend den Tragheitsmomenten der Riegel auf diese auf, so daB am Riegel m die senkrecht | 
zu seiner Achse angreifende Last 


Pa(*) = p(x) —= (4) 


ist. Sie ist in Abb. 3 symbolisch angegeben. Die RegelmaBigkeit im inneren Aufbau des Sy- 
stems wird auf diese Weise durch die RegelmaBigkeit im Lastangriff erganzt. Soweit die tat- 
sachlich gegebenen Lasten mit den rechnungsmabigen Lasten des Ersatzsystems in der Ver- 
teilung auf die einzelnen Riegel nicht identisch sind, entstehen in den Lamellen bzw. Pfosten 
zusatzliche Normalkrafte, die aus einer einfachen Gleichgewichtsbetrachtung folgen und den 
aus der elastischen Rechnung des Ersatzsystems resultierenden Normalkraften zu superponieren 
sind. Das gegebene Rahmensystem erleidet Verformungen infolge der Biegemomente, Quer- 
und Normalkrafte. Wie in der Rahmenrechnung tblich, vernachlassigen wir die Verformungen 
aus Quer- und Normalkraften. Dieselbe Vernachlassigung soll auch im Ersatzsystem Giltigkeit 
haben. Infolge der damit zugrunde gelegten Dehnungsstarrheit kénnen wir die oberen Riegel 
zusatzlich auBerlich gelagert annehmen. Die entstehenden oberen Lagerreaktionen sind dann 
z. B. fiir den gegebenen Fall, bei dem nur der unterste Riegel tatsdchlich gelagert ist, den er- 
mittelten Endpfostennormalkraften jeweils zuzuzihlen. Somit wird die Regelmafigkeit des 
inneren Systemaufbaus und des Lastangriffes noch durch die RegelmaBigkeit der duBeren 
Lagerungsbedingungen vervollstandigt. 


4, Ableitung des Naherungsverfahrens. Wegen der vorausgesetzten Dehnungsstarrheit der 
Lamellen sind die senkrecht zur unverformten Riegelachse auftretenden Verformungen (Riegel- 
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durchbiegungen) aller Riegel gleich. Dies bedeutet fiir die Lamellen gleiche Endverdrehungen 
an den Einspannstellen in die Riegel (Abb. 4, a). Es tritt auBerdem eine Relativverschiebung 
der Lamelleneinspannstellen in Achsrichtung der Riegel ein (Abb. 4, b). Bei seitlich verschieb- 
lichen Systemen kénnen sich die Riegel auch in Richtung ihrer eigenen Achsen verschieben, 
so daB auf diese Weise in jedem GeschoB zusatzliche Relativverschiebungen der Lamellen- 
einspannstellen auftreten (Abb. 4, c). Die Verformungslinien der Lamellen haben fiir die Ver- 
formungsanteile a, b und ¢ in halber lichter Lamellenlange einen Wendepunkt. Die Lamellen 
haben also in der Halfte ihrer Lichtweite einen Momentennullpunkt. Um ein statisch bestimmtes 
(oder auch statisch unbestimmtes) Grundsystem fiir unser Ersatzsystem zu erhalten, fiihren 
wir in jedem Gescho8 in halber lichter Lamellenlinge einen Schnitt. Entlang dieser Schnitte 
treten als unbekannte Funktionen die Querkrafte und Normalkrafte der Lamellen auf. Wir 
nennen die unbekannte Querkraftfunktion des m-ten Geschosses q,,(x) dx und die unbekannte 
Normalkraftfunktion des m-ten Geschosses n,,(x) dx (Abb. 5). Die Querkrafte des m-ten Ge- 
schosses seien positiv, wenn sie an dem zum Riegel m + 1 gehdrigen Lamellenschnittufer in 
Richtung der positiven x-Achse wirken. Die Normalkrafte seien als Zugkrafte positiv. 

Wir bezeichnen das Moment des Riegels m in- 
Mm -a(t)-das + $f hat tf folge der ihm zugeordneten duBeren Last p,,(x) bei 
Cs deca oe ero e ETT ANTI ATOT TD durchschnittenen Lamellen mit 


ca 2 Mo, m(x) = MG, m(x) te Aone (xp — x) + Boras (5) 


Riegel m 


ATV TNTTITTTTN Das Riegelmoment des Riegels m sei positiv, wenn 


es auf der dem Riegel m + 1 zugewandten Seite 


iE 


SEE ie ir satin aie: ER des Riegels m Zugspannungen erzeugt. M¢.,, ist 
or GeschokR m as Moment des Riegels m fiir irgendeine gedachte 
Mm(dc tt Pht tet t ait ea statisch bestimmte oder unbestimmte Riegellage- 
Gn (2-20 ee ee el rung. Ist die gedachte Riegellagerung mit der ge- 


gebenen Riegellagerung identisch, so sind Ao, , = 0 
und Bo, m = 0. Ist diese Identitaét nicht vorhanden, 
so wird die Differenz zwischen den Riegelmomen- 
ten fiir die gedachte und die gegebene Lagerung 
durch eine in x lineare Funktion erfaBt, die mit 
den Konstanten Ap,» und Bo,» festgelegt ist. Der 
Ansatz dieser linearen Funktion ist nach (5) so ge- 
wahlt, daB Ao, die lotrechte Lagerreaktion und 
Bo, m das Kinspannmoment des Riegels an der Stelle 
Xx = Xp infolge des durch Mj, ,,(x) nicht erfaBten Anteiles der Riegelmomente aus auBerer Last 
bedeuten. Ao,,, und Bo,,, sind positiv, wenn sie im Riegel m positive Momente erzeugen. 

Da in dem gewahlten Ersatzsystem die Endtangenten aller Riegel tibereinstimmen, miissen 
die Momentennullpunkte der Ersatzendpfosten wie die der Lamellen in halber lichter Pfosten- 
lange liegen. Die Endlamelle mit dem Tragheitsmoment I ,, dx/a, hat eine Lamellenquerkraft 
Ym(%4) dx baw. q,,(*g) dx; es folgt daher aus der Gleichheit der Verformung als Querkraft der 


Ersatzendpfosten 
Qf, m = UF) ig Im(% A) ? QF, m = US ig Im(*B) : (6) 


Bei der mathematischen Formulierung geht die Existenz der Ersatzendpfosten in die Rand- 
bedingungen ein: am Riegelende « = xg wirkt auf den Riegel m (Abb. 6) infolge der Querkrafte 
der Ersatzendpfosten das Moment 


Riege/ m+ 


Snot (Z)-an = ~<— = = <= | we ee 
Mm G2VVVV tty end 


Abb. 5. Unbekannte Funktionen Im(*) + dx und Ny (*) » dx 
als positive GréGen. 


aa [QB m— 1 Det =z OF m cI (Qo = 0 b) Q:, es 0) (7a) 
und am Riegelende x = x, das Moment 
ae [O04 m—1 bape a = Of | (04.0 4 0, Qt, = 0) D (7b) 


Es verhalten sich bei dem unserem Rechenverfahren zugrunde gelegten mechanischen Gedanken- 
modell (= Ersatazsystem) die Riegelmomente Mg,» proportional zu ihren Tragheitsmomenten, 
d. h. es ist 

My, (x) _ My,(2) 


I ae Qa rw S734 be pens 


I, (8) 


Bei allgemeinen Abmessungsverhaltnissen widersprechen sich (7a), (7b) und (8). Da Gleichung (8) 
die Voraussetzung zur Anwendbarkeit unseres Rechenverfahrens darstellt, fahren wir Hilfs- 


l,m 
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momente M# ,, und M®, » derart ein, daB die aus den Querkraften der Ersatzendpfosten her- 
riihrenden Momente nach (7a), (7b) an den Stellen « = x4 und x = x, so auf die vorhandenen 
r Riegel verteilt werden, daB die endgiiltigen Momente der Riegel des Ersatzsystems der Glei- 
chung (8) gehorchen (Abb. 7). An der Summe der durch (7a) bzw. (7b) gegebenen Momente 
iiber alle r Riegel tritt dadurch keine Anderung ein; es ist demnach 


DI Mg, =0; > MZ, =0: (9) 
y=] 


Die Momente M¥,,, bilden also an der Stelle x =x, und die Momente M%,,, an der Stelle 
x =x, Gleichgewichtsgruppen, die in dem gegebenen System einen Eigenspannungszustand 
hervorrufen. Die Einfithrung der Hilfsmomente ist daher qualitativ mit der Annahme zu ver- 
gleichen, da8 in der Balkenstatik die Auflagerreaktion durch parabelférmig iiber die Querschnitts- 
hohe verteilte Tangentialspannungen eingeleitet wird. Auch diese Annahme ruft:im Balken 
einen Kigenspannungszustand hervor, der nach dem Prinzip von de St. Venant etwa in einer 
Entfernung vom Auflager gleich der Balkenhéhe 
abgeklungen ist; dieser Eigenspannungszustand 
ist den aus einer Balkenrechnung ermittelten Be- 


TN al TT] 


i a ee GesohoB m-1 


& 
S 

7 Riegel m. 
iS 


t 
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Qe m-1 


La 
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| 
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Abb. 6. Querkrafte der Ersatzendpfosten. Abb. 7. Hilfsmomente am Ersatzsystem. 


anspruchungen zu superponieren. Da bei den von uns behandelten gegliederten Scheiben der 
Angriffsbereich der Gleichgewichtsgruppe bei hohen Scheiben die gleiche GréBenordnung wie 
die Spannweite erreichen kann, wird der Eigenspannungszustand nicht als rasch abklingend 
wie beim Balken vernachlassigbar sein; jedoch werden wir fiir ihn eine einfache Abschatzung 
durchfiithren kénnen. 

Wir behandeln zunachst das Ersatzsystem unter der Voraussetzung der Existenz der oben 
beschriebenen Hilfsmomente. Der Riegel m erhalt aus den Lamellenquerkraften und den Quer- 
kraften der Ersatzendpfosten der Geschosse m— 1 und m Momente. Wir bezeichnen das auf 
den Riegel m wirkende Moment aus diesen Kraften mit M,,,,(x): 


“B. 
My, m(x) = — [ [6m—1 dm—1(Q) + Cm Im())] 4 
aie Aa ™ (tp ae x) + Bow _ (O38, m—l1 bin —J] =| QF am c) e (10) 


(1 i, Go) —— 05 g,(7)) = 05 QF: 0 =0, QF, — 0) 
Der Riegel m erhalt aus den Lamellennormalkraften der Geschosse m — 1 und m Momente. 
Wir bezeichnen das auf den Riegel m aus diesen Kraften wirkende Moment mit M,, ,,(x): 
Mn, m(X) = ffm at) — rinlt)] ( — %) dy + (An,m + Nina — Nin) (%— — #) + Baym | (yy 
mn re (4) 05) a(n) = 0; IN = Oe INP =10);, 
Die angeschriebene Form der Konstanten wird aus ZweckmaBigkeitsgriinden gewahlt. Die 


Konstanten 4, , und B,,m lassen sich in die von n,,—; und n,, herriihrende Anteile zerlegen: 


7 ae — AC SA D) ispe =— Be ws os it eee . (12) 
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; A B 
AuBerdem erhalten die Riegel noch ein Moment My,,,(x) aus den Hilfsmomenten My, m und My, m: 


(13) 


B B A xp—x 
Mu, m(X) == Mz, j= (Mu, it oa Mya, ) a, . 
Damit lauten die Riegelmomente Mg,,,(x) im Ersatzsystem 


Mg, m(%) = Mo, m(%) + Mg, m(x) + Mn, m(x) + Mu, m(x) (lS msn). (14) 
GemaB (11) und (12) sowie (9) und (13) wird 


Me te) =0, > Mz,,(x) = 0, 
= 


v=1 v 

so daB als Summe aller Riegelmomente nach (14) 
3 Max) = 3 Mola) + 3 Mast) (15) 
vs eal y= 


folgt. Aufgrund der gemaB (4) vorgenommenen Aufteilung der gegebenen Last p(x) auf die 
Riegel wird 
T r 


Mo, m(x) =——=— 3) Mo,.(x) . (16a) 
a I, v oo 
y= es 
Wir fiihren als Abktirzung 
I 
in = (17) 
= I, v 
ein und erhalten aus (16a) 
Moca(*) =m orke) (16b) 
v=1 


Aus (8) folgt 
Mg, n(x) Le a Mg,,(x) ee : 


Die Summation dieser Gleichung iiber y liefert mit (17) 


Mr, me) == 11, m 2 Mg,,(x) . (18) 
Mit (16b) und (18) folgt aus (15) yi 
Mr, (x) = Mo, n(x) + ti, m a M,, (&) Si ee (19) 


Damit ist es gelungen, die Riegelmomente Mg,,,(x) des Ersatzsystems allein durch die unbe- 
kannten Querkraftfunktionen der Lamellen und die unbekannten Querkrafte der Ersatzend- 
pfosten auszudriicken. Wir fiihren eine Bezugssteifigkeit I, ein und definieren mit ihrer Hilfe 
die von m unabhangigen Momente 


Mg(x) 3 M,, m(*) Mo(x) ee Mo, m(X) = 
Le Lae ee ests: e) 
Als weitere Abkiirzung wahlen wir 
ig time (21) 
a 2 ie ” 


Aus den r Gleichungen (19) folgt damit eine einzige Beziehung: 
Mp(x) = M(x) + i, 2 Ms,(2) - (22) 


Wir setzen in (22) M,,,(x) nach (10) ein und erhalten unter Beachtung von a;,,, = b,, + ¢,, 


_ a r—1 *B Rol) 
M,(x) = M,(x) —i, a Bano i + t, Aq(%_— x) + 1, By —i, Dd) a1,,Q%,,~. (23) 
r= x a 


Darin bedeuten 


A, = on lap = 2 Ba» ‘ (24) 
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Die fir alle Riegel gleiche Durchbiegung y(x) geniigt der Balkendifferentialgleichung 


EI, y"(s) = — Ma,a(x) -< = — Mla). (25) 


Durch die Verformung y(x) der Riegel, an der sich die durchschnittenen Lamellen widerstandslos 
beteiligen, erfahren die Schnittufer eine Klaffung 6, (x) von der GréBe (Abb. 8a) 


01, m(x) con 1, ¥ (%) . (26) 

AuBerdem erfahren die Lamellen unter ihrer Querkraft q,,(x) dx eine Biegeverformung, die 
eine Klaffung 62, ,,(x) hervorruft (Abb. 8b) 
02, mite) 


a, IS m 
= a5 eT ie G(x) « (27) 
SchlieBlich kénnen sich bei seitlich verschieblichen Systemen die Riegel in Richtung der Riegel- 
achsen im m-ten GeschoB um die Strecke 03, m gegeneinander verschieben (Abb. 8c). 


| 
b,m(a)— bn (a) or 
a b Cc 


Abb. 8. Endverformungen einer in halber lichter Lange getrennten Lamelle. 
a) aus Riegeldurchbiegung; b) aus Lamellenquerkraft; c) aus seitlicher Riegelverschiebung. 


eo 


Die Kontinuitét des gegebenen Systems verlangt das SchlieBen der kiinstlich getrennten 
Lamellenschnittufer. Die Forderungen der Kontinuitaét lauten daher 


OfnAd) =F 02, nie) =e 03, m =a 0 ql Ss m = as 1) 0 (28) 


Nach Multiplikation mit E I,/a,,,, erhalten wir unter Verwendung von (26) und (27) 


“lim Ie Po pl Sm Sr 2 
Ely) +99 7, -im) +ELS™=0 (LSmSr—). (29) 
Wir definieren 
on B mie np 
Im(x) = ne oe m(*) » (30) 
= PAIS, ih, 
03, m oe x2, ine: 03, m . (31) 
Hiermit folgt aus (29) 
ae = =. 
E I, y'(*) + 45 [Gnd*) + 63m] =0 (LSmSr—1). (32) 


Die eckige Klammer in (32) mu von m unabhangig sein, da FE I, y’(x) und der Vorfaktor der 
Klammer von m unabhangig sind. Wir definieren daher eine von m unabhangige Hilfsfunktion 


G2) = 9%) +63, (Smsr—1) (33) 


und erhalten an Stelle der r — 1 Kontinuitatsbedingungen gem4B (32) als einzige Kontinuitats- 
gleichung 


2 
Bb Oly 


EI,y'(x) + G(x) =0. (34) 
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Wir differentiieren diese Gleichung einmal nach x und erhalten unter Verwendung von (25) 


= Dk a 
! (p\ te eae y= 0. 35 
d'(s) — yg, Male) = 0 (35) 
M,,(x) ist nach (23) in Abhangigkeit aller Querkrafte gegeben. Aus (30) und (33) folgt 
oe ioe 


g(t) =e?” (ax) -45,.] (1 SvSr—))- (36) 
Paks 


Die Gleichung (36) sagt aus, daB zwischen den Lamellenquerkraften verschiedener Geschosse 
lineare Beziehungen bestehen. Mit ihrer Hilfe formen wir die Summe 


*B 
eee xB r—1 
Dy E | an) iy = | 1E [41,5 anh dy 


x 
um; es wird 


=i} ES pele el 
fee pa 


= x I, y voy xt 2, ry 
2 a)» q,(y) = q(7) ron 2 I3 ; ihe Ba es if 3,” ° 

Wir fiihren ein: 

Stet, 

EB eal ie Ba 
a = 4 ey (37) 
: > qT, v 
y=1 
natal eee 
Dy iB 2,” 03, 
y =a (38) 
Le, 
v=) 
Aus der zweiten Gleichung (6) und aus (36) folgt mit (37) und (38) 
r—l 2 2 B . 
i Ds, a1,» Oz, = “Ey q(Xp) as is y : 
v=1 
Dann ergibt sich mit (23) 
<B 

aie a OP flee : : 4 aS 4 
M,(x) — M,(x) = — 208 | Gn) dn + (i, 4q—y) (%_— 2) + i, By— $20? iP G(xp) — a, By. (39) 


x 


My, nach (39) setzen wir in (35) ein und erhalten als Differentialgleichung fiir das Integral der 
gesuchten Funktion q(x) 


2B 


ay Oe fl 1s | F : oi : 

(2) + | eden = pg] Mole) + Ap —7) (to —2) +, By —Ty0PB Hea) —aaify]. (40 
Durch nochmalige Differentiation nach x finden wir als Differentialgleichung fiir die gesuchte 
Funktion q(x) 

" ff CEP 12 


q (x) TK q(x) = [M,(x)—¢,4,—y)] - (41) 


2 
xR a, 


Dies ist eine gewéhnliche Differentialgleichung zweiter Ordnung mit konstanten Koeffizienten, 
deren allgemeine Lésung mit der dimensionslosen Variablen 


sea (42) 
die Form 
: GE) = 4, (6) + C, Gof mE + Cy Sinwé (43) 
besitzt. 


Das Partikularintegral q)(x) hangt von M'(x) und der Konstanten (i, A, —y) ab. In y sind 
als Konstanten die unbekannten Seitenverschiebungen 63, mit 1 < »< r—1J1 enthalten. Die 


Konstanten A,, C,, C, und 63,, folgen aus den Randbedingungen. Fiir den Spezialfall, daB M,(x) 
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héchstens quadratisch von x abhangig ist, folgt aus (41) 


Gp(*) = —-—[ Mix) — (i, 4, —y)] (44) 


| al 

Damit haben wir die unbekannten Lameilenquerkrafte bestimmt, ohne die ebenfalls un- 
bekannten Lamellennormalkrafte zu kennen, die wir im Folgenden ermitteln werden. Aus (39) 
und der mit J,/J; , multiplizierten Gleichung (14) folgt durch Gleichsetzen 


xB 
Oy Ohh Sloe 8 
iM m(X) = SS mt} — My, mtx) ee : ince | q(7)) dij 
if Ti 5 | an ie Ty fe I m 
+ (im dq — Py) (tp =A) + 41,m By — gat Hea) — a BEB. (45) 
My, m(x) ist nach (10) gegeben; in (10) setzen wir gm— 1 und q,, nach (36) ein und erhalten 
xB 
= y — 5 M&, m™m a M4 m 
Mn, m(X) = Om “| q(y) dy + lism A,—Agm 03,m 4 = an (xp — 4) 
Bid 
si [t1, m Be poaae By, i iy 03, m Os: Vp A, is q(Xp) ae Mi, m] (46) 
mit 
=a ree a m—1 ~B | aaa Cm 941,m*~B Ta in ae a, Lie 
Om ; Cease Te a be ee | ees b Xp Ie ie 
7 
bo 4,9 Ly, 9 cra, ,J1,, 
er ee oe ee ae Naas 
(is Mf SIP i De B ) 
und 
Ue m Ue a), m—l1 x2, | Eee eee ry Cn a im x Ts im fat 
03, m = ip, y as [an Te 03, m—l } = ne Ts 03, m . (48) 
Wir fiihren als Hilfsfunktion n,(x) die auf den Riegel m wirkenden Normalkrafte ein: 
Tm(X) = n,,(%) — Mm —1(%) (n(x) = 0, n,(x) = 0), (49a) 
No Nea (Ni =0, Nf =0). (49) 
Damit lautet (11) 
x*B 
2 eB 
M,, n(x) = — J ta() (9 — 2) dn + (Anm — Non) (2p — %) + Brym (50) 


Gleichsetzen der rechten Seiten von (46) und (50) liefert nach zweimaliger Differentiation 
Nyn(X) pee G. XB q (x) : (51) 
Nach Einsetzen von (51) in (50) erhalten wir bei Beriicksichtigung der Gleichheit von 
xB ly *B Dera 
J 7,0) 9 — *) €n + On Xe J He) dy = Xp (%B— *) Om I(%x) 
aus der zu fordernden Identitat von (50) und (46) 
Anes rae Nz, ay (xp - x) =i Gre 
i i © of | Me a : 
= x5 (3 —*) Gn Wen) + [itm ty — Aarm Sim te |e) 
+ [i,m By — Bg, m — @ if 03, m + 2 12 XB Om I(%z) — MGs, ml] - 
Bei dem nun durchzufiihrenden K oeffizientenvergleich nehmen wir gleichzeitig eine Aufspaltung 


vor, und zwar definieren wir A,,m bzw. By,m als jeweils den Anteil, der aus den auBeren Reak- 
tionen der q,,-Krafte herriihrt, d. h. es gilt 


AG os — Uy, m As —— aw ’ (52a) 
Bim J ime B, 7 13 yes . (52b) 
24 
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Jedoch verbleibt auch dann, wenn die inneren Krafte q,, keine Reaktionen A,» erzeugen (also 
auch dann, wenn iiberhaupt keine Lager an der Stelle « = xg vorhanden sind), eine Konstante, 


é ‘ ; —=B 2 , 
die mechanisch als singulare Endnormalkraft N,, zu erklaren ist: 


B A 
Hm — Min m (52c) 
ds | 


Ne = — Om Xp Q(%p) + OF, m 


Weiterhin resultiert aus diesem Koeffizientenvergleich die Bestimmungsgleichung fiir das Hilfs- 
moment MF, n 


Mit,m = 4 1 [Gm *p Gn) — 08, ml « (52d) 
Aus Gleichgewichtsbetrachtungen folgt schlieBlich 
se ME — Ma 
Na = Gn Xp G4) — 08 m + — (G20) 
B A 
Mit, m = — 2 12 [0m %p O44) — 98m] - (52f) 


Gleichung (52d) und (52f) folgen auch aus der Summation von (7a) und (7b) tiber m von 1 bis r 
und anschlieBender Aufteilung dieser Summe im Verhaltnis der Riegeltragheitsmomente. Es 
wird gemaB (52d) und (52f) 


Miz, m — Mit, m = 42 Om Xp [iF G(x) + if G(x4)] — 68, m (12 + if)} - (52g) 
Nach (49a) gilt 


Da n,(x) = 0 ist, folgt 


Mit der Definition 
on 2: ey (53) 


ergibt sich fiir die Lamellennormalkraft explizit 
Mm(X) = Om Xp q (%) - (54a) 


Analog folgen als singulare Endlamellennormalkrafte an den Stellen x = x, und x = xp 


Ai = a pa fit a : 5 Ze 
Nin = Om 0 fala) + te UE Mes) + aleol| —[L + E+ if)| Soh. Ah) 
Beams = ay Bi tes a : : ie 
Nin = — Onn fle) + 8 (8 den) + i alea)l} + [1 + OB + is] Sot, Ae) 
mit 
te * __ xB Bl ah dose y ue Ci, 9, m Lom = 
2 OB» I, B By Oa x2, 2 Ty i ane m °s, | i Oy 


Fir das Ersatzsystem haben wir vorstehend die unbekannten Quer- und Normalkrafte der 
Lamellen ermittelt. Die eigentlich gesuchten unbekannten Quer- und Normalkrafte der Pfosten 
des gegebenen Systems erhalten wir hieraus durch Integration. Wir bezeichnen die Querkraft 
eines Pfostens an der Stelle x, im m-ten GeschoB mit Qp, »(x,), die Normalkraft dieses Pfostens 
mit Np, ,,(%,). Entsprechend den mit einem Vorfaktor nach (30) multiplizierten Lamellenquer- 
kraften q,,(x) fiihren wir auch Pfostenquerkrafte Qp, m(X;,) eins 


= I om Le 

Op, m(%) = x2, ane 1s Qp, m(*,) . (55) 
Es wird mit (33) fiir die Innenpfosten 

Qr, m(%,) = q(x) dx — ay 53, m (56a) 
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und fiir die Endpfosten 


x4+ = 

Qrmbsig) =f Ge) de + a if Gg) —$25s,m (1-42 if), (56b) 
a7 | 

— . 

Qrm(s) = { Gla) de + a5 iB Gp) —$°Bo,m (1 +218). (560) 
ae 


Fir ein System mit AuBenpfosten, die nicht am Riegelende angeordnet sind (Abb. 2), erfahren 


die Integrationsgrenzen eine entsprechende Anderung; der Vorfaktor von 63 » ist gleich der 
Lange des Integrationsbereiches. 


Mit der dimensionslosen Abkiirzung 


A = 2 Gin = (57) 
folgt mit (43) aus (56a) bis (56c) 
xy +2 
Qp, m(X;,) = i; G p(x) dx — a, 5s, m + aA (Cy Coj a &, + C, Sina E,) » (58a) 
x, — 
zat 


a, 2s 2 ow 2 aie p 
ES A\G, (1+ ; J+, 393% Coats +[C,%g 2% + G,(1 + ; Sin aE, (58b) 
*B 
De.nlin) = | Gz) de + 05 Gen) — $88s,m (1 +28) 
eee 
a, 24" & Ay &% Ay oe 4 
+#alla (a rey C,Ig%%| Cof a OV. Ge 7 cane (58c) 


Die tatsadchlichen Pfostenquerkrafte Qp, ,,, folgen aus den bezogenen Querkraften Qp,m nach (55). 

Die Pfostennormalkrafte des gegebenen Systems, an dem zusatzlich die in (52d) und (52f) 
angeschriebenen Hilfsmomente wirken, nennen wir N},,,(x); sie lauten, wenn die AuBenpfosten 
am Riegelende angeordnet sind (Abb. 1), fiir die Innenpfosten 


xp + rt 
NE n(x,) = | n(x) dx (59a) 
—— 
und fiir die Endpfosten 
x4t < 
Nb (x4) = f  mp(x) dx + Ne, (59b) 
“A 
eB 
Nt, (xz) = | n,(x) dx + Nn. (59c) 
ee 


Fiir ein System mit AuBenpfosten, die nicht am Riegelende angeordnet sind (Abb. 2), er- 
fahren die Integrationsgrenzen von (59b) und (59c) eine entsprechende Anderung. AuBerdem 
wirken die singularen Endlamellennormalkrafte N4 und N® nicht an den Stellen der AuBen- 
pfosten, so daB sie nicht der Normalkraft dieser AuBenpfosten superponiert werden dirfen. 


24* 
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Vielmehr sind sie mit umgekehrter Wirkungsrichtung an dem gegebenen Rahmen anzusetzen. 
Sie bewirken in allen Pfosten Quer- und Normalkrafte, die den bereits ermittelten cae a4 
superponieren sind. Mit guter Naherung kénnen die Pfostenquerkrafte a tea fs t 
vernachlassigt und die Pfostennormalkrafte als Reaktionen eines an den Ste eh ie a en 
fest gestiitzten und frei drehbaren Durchlauftragers errechnet werden, der an den Randern 
x =x, und x = xg wie die Riegel des gegebenen Systems gelagert ist. 


Mit (54a) bis (54c) folgt aus (59a) bis (59c) 
ilo) al —4) sou 


N?, (Xp) = Om ~B 


NP, m(24) = Om Xp i (x ue a +t Lid dla) + if ale) 
=" eee ee ec ‘| 8, (59e) 
Xp—_%X4 v=l 
N8, m(%z) = — Om *e i (*0— St) + spc, LP Wen) + if ate} 
tee ee ae (59f) 
*Bp—_ 4 oad 


Mit (43) ergibt sich dann 


Nb, m(%x) = Om ~B ? (60a) 


Gp (% + $)— 9 (%#— 7) + oA e(Cy Sina & + Cy Gof a &) 


NB. nla) = Om ®a [Ip (¥4 + $') + yg tw, Li dln) + i Gpleal 
[Xa(Ea tat) > + spot e, OF Coho + if Gof a ey 
B " oi | 


7 Gy ce Oe 


PA 
‘ 
+ C, gin (a (Es eG 219 + ie (i2 Gin « + i4 Gin « 4) 


xB x4 


| 


a 


ame pal 
NB, nls) = — Ome flip (*o— 9) yg 8 LE Glen) + 18 GC] 


4G; (Coy Ge —o » +" (iP Cof a + if Cof a €4) 


Gr jen (a(1 aL |} — “2 _ (2 Cina + i4 Gina £0) 
XK “4 | 


a, 


m 

ae [1 ee i) Dba (60c) 

Es wurde vorstehend das Ersatzsystem unter der Voraussetzung von auBen angebrachter 
Hilfsmomente M4, und M¥,,, behandelt. Die GréBe dieser Hilfsmomente folgt aus (52d) und 
(52f). Da im gegebenen System die Hilfsmomente nicht existieren, superponieren wir dem von 
uns am Ersatzsystem ermittelten Beanspruchungszustand einen weiteren Zustand, den wir durch 
Anbringen der Hilfsmomente in umgekehrter Wirkungsrichtung am gegebenen System erhalten 
(Abb. 9). Wir haben bei der Einfiihrung der Hilfsmomente erlautert, da® sie Gleichgewichts- 
gruppen bilden. Fir diesen Kigenspannungszustand vernachlassigen wir die Biegesteifigkeit 
der Pfosten. Dann erfahren die Riegel an den Kreuzungspunkten ihrer Systemlinien mit denen 
der Pfosten keine Vertikalverschiebung und wir kénnen sie als Durchlauftrager behandeln, die 
an ihren Enden x = x4 und x = xg wie die gegebenen Riegel gelagert und mit den Momenten 
— M4, und — M¥,,, belastet sind. Die aus dieser Durchlauftragerrechnung resultierenden 
Auflagerkrafte nennen wir Ny, m(%,)3 sie sind in Abb. 10 als positiv wirkende GréBen einge- 
tragen. Die daraus sich ergebenden Pfostennormalkrafte bezeichnen wir mit Nu, m(%,)3 es gilt 


Nu, m(%) = Nu, m(%) — Nu m—i(%)  (LSm<r, Nyo=0, Nu,=0). (61) 


Hieraus folgt 


Nu, n(x) = >' Ng, (x) - (62) 
y=1 
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Die endgiiltigen Normalkrafte der Pfosten des gegebenen Systems lauten 
Np, m(%y) = NB, m(%q) + Nur, m(%,) « (63) 


Kine eventuell vorhandene geringere Genauigkeit der auf diese Weise ermittelten ZusatzgréBen 
gegentiber den aus dem Ersatzsystem ermittelten BeanspruchungsgréBen ist ohne Bedeutung, 
da dieser Zusatzlastfall an sich nur einen geringen Anteil gegenitber dem im Ersatzsystem 
behandelten Fall hat und daher die méglicherweise vorhandenen héheren Fehlerprozente auf 
einen kleineren prozentualen Anteil an der Gesamtlésung zu beziehen sind. 

Bei Systemen mit zwei Riegeln, deren Tragheits- 
momente und Systemlinienabstande von der in halber 
: Lichtweite gedachten Schnittlinie gleich sind, d. h. 


ht — Li-2 5 


D=Iy D=Ipg 


Abb. 9. Umgekehrte Hilfsmomente am gegebenen 


Abb. 10. Durchlauftraiger zur Ermittlung der Beanspruchungen aus den 
System. 


Hilfsmomenten. 


gilt, verschwinden die Hilfsmomente Mie, und Mie die Lésung mit Hilfe des Ersatzsystems 
bedarf in diesem Falle keiner Zusatzlésung am Durchlauftrager.! 

Die Biegelinie y(x) der Riegel des Ersatzsystems erhalten wir durch Integration der Kon- 
tinuitatsbedingung (34) 


2 
x, Ae 


E Iya) + 2 [fan dy + Cs] = 0. (64) 
Setzen wir in (64) q(x) nach (43) ein, so ergibt sich 
E Ty(x) =— "4%" | (4,69) dy + Cy + 72(C, Ginak + Cy Coja8)]- (65) 


5. Systeme mit gleichen Abmessungen in allen Geschossen als Sonderfall. Die Lésung fiir diese 
Systeme erhalten wir durch Spezialisieren der oben abgeleiteten allgemeinen Lésung: 


(DF ore et Si Ss), 5 

Inm = I, , 

n= iy, 

i lee 

ta les (66) 
ie ,\lsmsr—})). 

Q1,m — 4, 

Ce 

bb, 

oon = les 


Es wird 


r 
eli = Tels 
v=1 


und es folgt aus (4) als gegebene Last je Riegel 


1H, Beck, Bauing. 31 (1956), S. 436. 
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Fiihren wir die Bezugssteifigkeit I, = I, ein, so folgt aus (20) 


M,(x) = Mz(x) = Mz,n(x) ; (67) 


Nach (30) und (31) wird 


Gee de ee 

Im(X) vie a a, is In(X) ? ( ) 

bs, i = ee z Os, m* (69) 
x a, a, 


Die Kontinuitatsgleichung ist unverandert nach (34) gegeben. Nach (37) und (38) wird 


~r—1 12 x3, a? I, 


2 70 
% r Ba, I, ee, 

1 *Raj I, 7—!~ 71 
ie ee 2 


Mit diesen Konstanten ist die Differentialgleichung unverandert durch (41) gegeben. Das 
gleiche gilt fiir die Lésung q(x) (43) und die Lamellennormalkrafte n,,(x) (54). Es lautet speziell 
die Abkiirzung 0,, (47) jetzt 


2 E 
ee sr = By (72a) 
2a 
2) 1 b) apf, 
Qa (2 -) claws (72b) 
= 1 c \ ax, I, 
Qe; = : <) I 16 (72c) 
und damit @,, nach (53) 
A er 1h 
&=(2—2) 7 Gamer) (73) 
B 1 
und nach (54d) 
Le CBE I Ne eel = 
Do Soe >: —— > Bs, ¢ + — Oa . (74) 
y=1 I, = yall ay 


6. Systeme ohne Ersatzendpfosten als Sonderfall. Fiir diese Systeme gilt, wenn sich die 
AuBenpfosten am Riegelende befinden, 


Lon = Lim => Tam (75) 
Daraus folgt gemaB (1) 
Tf olen (76) 
sowie gemaB (3) 
id — iP =0. (77) 
Es entfallen somit auch die Ersatzendpfosten-Querkrafte (6): 
Qi, m = QE,m = 0. (78) 


Folglich sind auch die Hilfsmomente fiir diesen Sonderfall nicht existent, was sich formal auch 
aus (52d) und (52f) ergibt, namlich 
Mi, m = Mim =0. (79) 


Wahrend die Formeln fiir die Riegelmomente bei durchschnittenen Lamellen infolge auBerer 
Last gemaf (5) und Lamellennormalkrafte gemaB (11) unverandert Giltigkeit haben, lautet das 
Riegelmoment aus den Lamellenquerkraften nach (10) jetzt 


Mym(*) = — i! [Bm —1 Gm—1() + Cm Gnl(9)] 49 + Agsm (Xp —%) -+ Bam 


(80) 
SmSr, qq) =0, 9,(y) =). 
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Konsequente Durchfiihrung dieser Spezialisierung liefert fiir die bezogenen Pfostenquerkrafte 


nach (58a) bis (58c) nunmehr 


ay +S 
Op, ma) = [  Gplx) du — a, 5s,m + 0,4 (C, Cof x £, + Cy Gina), (81) 
aks 
wa+2 
Or, n(x4) = | Gp(2) dz — 225s, y 
| 
Peay) (C a 6 Tg T)Cof og Hi Cpige = C,\ ome (81b 
5 1 22975, A 1 Sax, 2 Al? ) 
_ am 
Op, m(Xz) = | Gp(*) du — 255, m 
ee 
+ $4|(C —G Tate) Gofa +(— G, Tape + Ca) Sina. (81¢) 


Die tatsachlichen Pfostenquerkrafte Qp, » ergeben sich aus (55). Wegen des Wegfalls der Hilfs- 
momente entfallt auch die Zusatzlésung am Durchlauftrager. Folglich ist Ny, m(x,) =0 zu setzen, 
so daB (63) nun lautet 


Np, m(%,) = NP, m(%) - (82) 
Damit ergibt sich fiir die Normalkrafte gemaB (60a) bis (60c) jetzt 


Gp (% + 3) —Go(* 2 + oA (C, Sin o &, + Cy Cof a &;) 
B 


Np, m(X;) = Om *B , (83 a) 


Np, m(%4) = Qm*B 


— a a : m 

dp (x ae =| + C,Coj (« (E, 4 =). + C,Gin (a (E, ae eb) — 3 9t,,, (83) 
zie Up go ; a * 

Go| %2—3) + C, Cj (a ea) 4+C,Gin (a(1 = Ab + 3.68,. 83¢) 

7. Zusammenfassung. In der vorliegenden Arbeit ist ein Verfahren zur angeniherten Be- 
rechnung regelmaBig gegliederter Scheiben entwickelt. Dabei wird von der RegelmaBigkeit der 
gegebenen Systeme und der dadurch bedingten Vielzahl gleichartiger statisch unbestimmter 
Groen ausgegangen und auf Grund mechanischer Uberlegungen das gegebene Rahmensystem 
durch ein Kontinuum ersetzt. Die inneren Krafte des Ersatzkontinuums sind mit den tiblichen 
mathematischen Methoden ermittelt; es ist gelungen, alle im Ersatzkontinuum unbekannten 
GréBen auf eine einzige unbekannte Funktion zuriickzufiihren, die sich aus einer gewéhnlichen 
linearen Differentialgleichung zweiter Ordnung ergibt. Die gesuchten statisch Uberzahligen des 
gegebenen Rahmensystems erhalt man durch Integration der entsprechenden GroBen des Er- 
satzkontinuums. Wahrend die herkémmlichen Methoden der Rahmenstatik die Lésung linearer 
Gleichungssysteme mit hoher Gleichungsanzahl erfordern, bedarf das vorliegende Naherungs- 
verfahren nur der Lésung einer Differentialgleichung von bekanntem Typus. 

Bei der im heutigen Bauen zunehmenden Verwendung regelmafig gegliederter flachenhafter 
Fassaden liegt die Bedeutung des hier entwickelten Verfahrens darin, daf fiir derartige Systeme 
mit gleichen aueren Lagerungs- und Belastungsbedingungen die statisch tiberzahligen GréBen 
mittels Formeln angegeben werden kénnen, worin die speziellen Systemabmessungen als Para- 
meter eingehen. Fiir die praktische Handhabung werden Gebrauchsformeln fiir die am haufig- 
sten vorkommenden Lagerungs- und Belastungsfalle an anderer Stelle mitgeteilt. 


Np, m(Xp) = — Om*B 


(Eingegangen am 20. Dezember 1957.) 


Anschrift des Verfassers: Privatdozent Dr.-Ing. Hubert Beck, Frankfurt am Main, Hermannstr. 31. 
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Turbulente Ausbreitung runder Heifluftstrahlen in bewegter Luft 
(Erster Teil: Kernbereich) 


Von W. Szablewski 


1. Einleitung. Grundlage der Untersuchung bildet das im Jahre 1947 aufgestellte Gleichungs- 
system! zur Berechnung turbulenter Strémungen von Gasen stark verdnderlicher Dichte. In einer 
kiirzlich erschienenen Veréffentlichung? habe ich als Vorstufe der vorliegenden Arbeit die turbulente 
Vermischung ebener Heifluftstrahlen, die zunachst durch eine Scheidewand getrennt sind und sich 
dann mischen, fiir Unterschallstrémung berechnet. Dieses Modell wird durch das in unmittelbarer 
Nahe der Diisenkante befindliche Mischungsfeld eines Diisenstrahls der Machzahl kleiner als Eins 
mit dem umgebenden Medium realisiert und liefert mathematisch die Ausgangsposition (Anfangs- 
profile) fiir das in der Arbeit anzuwendende Fortsetzungsverfahren zur Lésung des Systems nicht- 
linearer partieller Differentialgleichungen, die das Mischungsfeld des Diisenstrahls in Strahlen- 
koordinaten, die dem Ausbreitungsvorgang angepaBt sind, beschreiben. 

Historisch ist anzumerken, da® mit Beschrankung auf Felder konstanter Dichte die turbulente 
Ausbreitung eines runden Diisenstrahls in ruhender Luft erstmalig von A. M. Kuethe® berechnet 
wurde. Es folgten dann H. B. Squire und J. Trouncer*, die in die Untersuchung auch den Fall be- 
wegter AuBenluft miteinbezogen. FuBend auf der Grundlage des Prandtlschen Mischungsweg- 
ansatzes, wurde von den genannten Autoren methodisch so vorgegangen, da eine plausible Funk- 
tion der Geschwindigkeitsverteilung iiber der Strahlbreite angenommen wurde, mit der man dann 
in den Impulsintegralsatz einging. Auf diese Weise wird das Problem zuriickgefiihrt auf die Lésung 
von gewohnlichen Differentialgleichungen fiir die das Feld festlegenden Funktionen der Misch- 
breite und der Kontur des Strahlkerns®. Es handelt sich also um ein Verfahren, das im Prinzip 
analog ist dem bekannten Pohlhausen- Verfahren in der Theorie der laminaren Grenzschichten. 

Demgegeniiber ist das in der vorliegenden Arbeit zur Anwendung kommende Fortsetzungs- 
verfahren, bei dem wir, ausgehend von den Anfangsprofilen, mittels des Systems von Differential- 
gleichungen das Feld in Differenzenschritten approximativ aufbauen, als exaktes Verfahren anzu- 
sprechen. Kin Differenzenverfahren ist in jiingster Zeit auch von S. J. Pai® fiir die Berechnung der 
laminaren Ausbreitung von heiBen Diisenstrahlen in Beispielen angewandt worden, wobei als eine der 
unabhangigen Variablen die Stromfunktion benutzt wird. Zum Rechengang selbst, z. B. iiber die 
GréBe der Differenzenschritte, wurden dabei keine Angaben gemacht. Das von Pai formulierte 
Gleichungssystem fiir die turbulenten Mischungsfelder von Gasen stark veranderlicher Dichte, 
das von unserem Gleichungssystem wesentlich verschieden ist (vgl. Arbeit II), ergibt fiir das Prob- 
lem der turbulenten Ausbreitung formal dieselben Differentialgleichungen wie bei der laminaren 
Ausbreitung, so daB die dort gewonnenen Ergebnisse auf die turbulente Ausbreitung tibertragen 
werden kénnen. Ein Beispiel einer derartigen Ubertragung ist jedoch noch nicht mitgeteilt worden. 

In der Durchfiihrung des Fortsetzungsverfahrens haben wir uns fiir den Kernbereich mit drei 
Schritten begniigt. Diese Beschrankung wird nahegelegt durch die sich aus bisherigen Forschungen? 
ergebende Feststellung, daB innerhalb des hier untersuchten Kernbereiches die Isotachen und Iso- 
thermen nur schwach gekriimmt verlaufen baw. die Profile der Geschwindigkeit und Temperatur 
in Strahlenkoordinaten nur eine geringe Verschiebung vom Ausgangsprofil erfahren und fast inzi- 
dieren. Unsere Rechnung ergibt denn auch, daf fiir die Erfassung des Kernbereiches praktisch 
sogar schon zwei Schritte als ausreichend anzusehen sind. 

Die Untersuchung beschrankt sich auf den Kernbereich, d. h. auf den Bereich des Mischungs- 
feldes, der durch das Vorhandensein eines Strahlkerns gekennzeichnet ‘ist, der durch die Vermi- 
schung noch nicht beriihrt ist. Wir werden dabei gemaB dem Umstand, da die Temperatur des 


1 W. Szablewski, Ing.-Arch. 20 (1952), S. 67. 

° W. Szablewski, Ing.-Archiv. 25 (1957), S. 10; diese Arbeit wird im folgenden als Arbeit II zitiert. 

° A. M. Kuethe, Journ. Appl. Mech. 2 (1935), S. 87. 

* H. B. Squire and J. Trouncer, British ARC R and M 1974 (1944). 

y Hierher gehéren auch die alteren, unter stark vereinfachenden Annahmen durchgefiihrten Arbeiten des 
Verf. tiber Felder konstanter und veranderlicher Dichte: W. Szablewski, U. u. M. 8003 (1944) baw. NACA TM 
1311 (1951); sowie Gottinger Dissertation 1947 bzw. NACA TM 1288 (1950). 

i S. J. Pai, Journ. Aeron. Sci. 16 (1949), S. 463; Quarterly Appl. Mathem. 10 (1952), S. 141. 

Vgl. z. B. A. M. Kuethe a. a. O.; ferner O. Pabst, U. u. M. 8004 (1944). 
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Diisenstrahls nach den Seiten starker ausgebreitet wird als seine Geschwindigkeit, fiir das Tempera- 
tur- und Geschwindigkeitsfeld voneinander verschiedene Kernbereiche zu unterscheiden haben. 
Grenzschichteinfliisse der Diisenwandung bleiben auBer acht, und die Profile der Temperatur und 
Geschwindigkeit in der Ebene der Diisenmiindung werden als rechteckig angenommen. Die Mach- 
zahl des Diisenstrahls wird, wie schon erwahnt, kleiner als Eins vorausgesetzt. 

SchlieBlich bemerken wir noch, daB wir hier nur den Fall bewegter AuBenluft betrachten. Der 
Fall ruhender AuBenluft, der die ganze Untersuchung abschlieBen wird, soll in einer nichsten Ab- 
handlung folgen. Er erfordert gegeniiber den anderen hier behandelten Fallen eine wesentlich gréBere 
numerische Genauigkeit und scheint auch sonst Abweichungen vom Verhalten dieser Falle zu zeigen. 
Diese Griinde lieBen es geboten erscheinen, den Fall ruhender AuSenluft besonders zu behandeln. 


2. Gleichungssystem. Das turbulente rotationssymmetrische, im zeitlichen Mittel stationare, 
Mischungsfeld (Abb. 1) wird unter der Voraussetzung konstanten Drucks durch folgende Gleichun- 
gen! beschrieben: 


Impuls: 

Ot Or Oey -* ef = du — 00 

gn ul ae 8 = ol) 5 ("05 + Eri at (1) 
Energie: 

OrQucpT erovcpT ~ Kee facocp T 9 

ae a ae)? @) x 
Masse: cee 
Aye 7 si Aye ra a = : a 0 Abb. 1. Geschwindigkeitsfeld. 
[oe or os eens oe a al (3) 


Hierbei bedeuten x, r Zylinderkoordinaten, u bzw. v die Geschwindigkeitskomponenten in der x- 
bzw. r-Richtung, 0 die Dichte, T die absolute Temperatur, c, die spezifische Warme bei konstantem 
Druck; FE ist eine empirische Zahl, die das Ubertragungsverhaltnis von substantiellen EKigenschaften 
und der (Druckschwankungen wahrend des Transports ausgesetzten) Geschwindigkeit beim turbu- 
lenten Austausch angibt. Das Uberstreichen bedeutet zeitliche Mittelwerte; im folgenden werden 
wir die Querstriche fortlassen. 

Die scheinbare kinematische Zahigkeit wird durch 

€(x) = x, b,(x) (Ug — uy) (Ug > Uy) (4) 

beschrieben, wo b,(x) die Mischbreite des Geschwindigkeitsfeldes, x, ein empirischer Koeffizient, 
Uy baw. u, die Geschwindigkeit des Diisenstrahls bzw. der AuBenluft bedeuten. Entsprechend be- 
zeichnen im folgenden T, baw. T, die Temperatur des Diisenstrahls bzw. der AuBenluft. 

Die Gré8e c, wird als Konstante angesehen. Mit der Zustandsgleichung fiir ideale Gase bei kon- 
stantem Druck 


o: T = konst. (5) 
ersetzen wir 0 durch 3, wo ? = T— T, : 
__ konst. konst. es 
= eee. ITs: 


Fiihren wir die dimensionslosen Funktionen g = (u — u,)/(ug — u,) sowie y = 8/99, wo Vy = 
T, — T,, ein und machen die Langen mit dem Disenradius ry dimensionslos, so kénnen wir 
unter Beachtung von (3) das Gleichungssystem (1) bis (3) auch schreiben: 


Impuls: 
. r/Tro [Yor uy) op , uv op 
1 + (8o/T,) % ( Ug a) Ax/ry ' Up Or/ro 
= Ae a) oye rlro Op r 0p H/T, ax (6) 
—= by Or/rg\1 + (8o/T,) x arira} Tea Orino (1 + (9/T,) x)? Gniroie 
Knergie: 
r/Tro Uy —_U, mea OF ae ok 
1 + (80/7) % Ug = Uy] Ax/ry | Up Or/ro 
i a r/ro LES a ey ee the oy 7 
= Eee) (re wnyg ai ete © 


1 Vel. Arbeit II. 
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Masse: 
ae ) r[ro (Up : U\ lig, eae r/ro v 
ax/ro Feaciaar: Up Up ' @r [tg Nl ae (8o/T) % Uo 
r) r Onde ax 
— E &,(x) Or/ro i ro (1 -++ (8/T 3) Fa ie aks) (8) 
mit 
b, up —u 
&(x) = % = : Us : 


Da wir es mit einem Ausbreitungsvorgang zu tun haben, ist es zweckmafig, auf Strahlen- 
koordinaten 


x 
mit 
(2/@r), =~ (A/en)e, (2/2), = (02) — (6/6) 
zu transformieren. 
Fithren wir noch mit dem Streckungsfaktor 
1 
0 Yim e; (lg — te)” i 


: ¢ Lx : : v : 
wo ¢, =6,/x, die Variablen & = >— sowie y =o7 ein und schreiben y = a, so erhalten wir 
inte 0 


schlieBlich folgendes Gleichungssystem : 


Impuls: 
ap op Up — Uy Uy; (8 /T,) x] Oy = § Ug— Uy , Uy ap 
ay? T ay |( Up ae a \y—y| ee NT IT) x Cues, 2/ Uo ae Uo 5 0g oe 
(11) 
Energie: 
OY 0% [2 [/uo— wy Uy (Fo/ Ty) 0/0 g _ 2 (uy Uy\ 2 OX _ 
dy? aE | Uo ey ayy ¥| * T+ @alTy) x ! dae y B| Ug eo Ug S96 a a 
Masse: 
Uy—U, /,. OP op ow oa 
Uo (= vx) Pap aaa Us 


Gegeniiber dem Feld konstanter Dichte treten die unterstrichenen Glieder zusatzlich auf. 
Dabei haben wir vorausgesetzt, daB 
o(x) = konst. (14) 
In Anbetracht dessen, da man iiber den nach (10) ino enthaltenen empirischen Koeffizienten x, 
derzeit noch nicht viel sagen kann, ob es sich etwa um eine Konstante universellen Charakters 
handelt oder auch nur um eine Modellkonstante, erscheint diese Annahme hier als die verniinftigste. 
Man kann sie so lesen, da fiir das Mischungsfeld des Kernbereiches x, und b,/x als konstant ange- 
sehen werden. Beides ist fiir den Kernbereich als in guter Naherung zutreffend anzusehen; denn 
beide Annahmen werden durch die experimentelle! und theoretische (vgl. Abb. 16) Feststellung ge- 
rechtfertigt, daB die Geschwindigkeitsprofile in Strahlenkoordinaten fast inzidieren, so da8B also 
tatsachlich o(x) im Kernbereich bis auf einen kleinen Fehler als konstant angesehen werden kann?. 
Die Randbedingungen fiir eine Lésung unseres Gleichungssystems (11) bis (13) im Kernbereich 


Jauten bei Beachtung, daB die Strahlachse die Strahlenkoordinate 1 = — 1/(x/r9) baw. y = — 1/é 
hat, 
1 
1 ee 
o> fiir <5 (15) 
0 aA 
1 > 
i fiir Je (16) 
0 ves 


1 O. Pabst, U. u. M. 8004 (1944). 
* Strenge Inzidenz wiirde besagen, daB erstens b, ~ x und da zweitens die Profile sich ahnlich verhalten 
bzw. daB x, = const. Denselben Schlu8 kann man auch bei einer Parallelversetzung der Profile ziehen; als 


Ahnlichkeitskoordinate fungiert dann 1* = Coe dito _ L—d/ro . . ; 
des Strahlkerns feackicibee d b/ro vs brane wo die Funktion d(x) die Kontur 
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Die jeweils erstere Annahme besagt, daB ein Strahlkern noch vorhanden ist. Wie schon elngangs 
erwahnt, werden wir zwei voneinander verschiedene Kernbereiche des Geschwindigkeits- und 
Temperaturfeldes zu unterscheiden haben. Das Geschwindigkeitsfeld wird dabei eine groBere 
Kernlange haben als das Temperaturfeld; wollen wir das Temperaturfeld auch noch iiber seinen 
Kernbereich hinaus, dabei aber immer noch innerhalb des Kernbereiches des Geschwindigkeits- 
feldes, berechnen, so ist die erstere Randbedingung in (16) zu modifizieren gemafB dem Umstand, daB 
hier bereits die Temperatur lings der Strahlachse abfallt; vgl. (48). 

Weiter ist auf der Strahlachse Verschwinden der Querkomponente der Geschwindigkeit anzu- 
nehmen, was die Randbedingung 


yo fiir Daa ae (17) 
ergibt. 
Fir die Diisenmiindung mit § = 0 ergibt (11) bis (13) die Gleichungen 

uw 1 \o1(4o__ 4 peel : (9o/T) Xo oa 
Po se Po P ( Uy Po 1 ae vo) (E 1) ie (9 /T,) Xo =i ? (18) 

ye {2 [/ug—uy, ma ee ta) yo 
Xo + Xo tz ( Up Po 1 abe vol 2 1+ (4/T,) A =0, (19) 

Uy — U t t 

ad eo Po + YY =9. (20) 


Das ist gerade das Gleichungssystem, das die turbulente Vermischung zweier ebener Heifluft- 
strahlen beschreibt, die zunachst durch eine Scheidewand getrennt sind und sich dann mischen, und 
das wir in der Arbeit II gelést haben. Dieser Arbeit kénnen wir also die dort berechneten Profile 
der Geschwindigkeit bzw. Temperatur @p, yy bzw. 7) als Anfangsprofile entnehmen: 
Anfangsbedingungen: 

P(E, ¥) > Poly) 

x(E¥) > Holy) dur oo 0) (21) 

plE, ¥) > Poly) 


Was die empirischen Koeffizienten der Theorie betrifft, so setzen wir fiir das Ubertragungs- 

verhaltnis, wie in Arbeit IT, 
Pa 2, (22) 

Nach (10) haben wir fiir den Streckungsfaktor o die Gleichung 

= 1 
2 x, (0 Cy) (Up —Uy)/Up * 

Gem48 unserer fiir den Kernbereich getroffenen Annahme (14) o() = konst. kénnen wir dann o(x) 
der Arbeit II (o — o(0)) entnehmen. Dort hatten wir die Breite o c, der Geschwindigkeitsverteilung 
im y-MaBstab zwischen den Grenzen 0,95 und 0,05 der berechneten Verteilungen abgelesen. Der 


Vergleich mit den Messungen von O. Pabst! ergab dann fiir den empirischen Koeffizienten x, den 
Zahlenwert x, = 0,0082. 


o} 


3. Fortsetzungsverfahren. Zur Lisung des Systems partieller Differentialgleichungen (11) bis 
(13) mit den Randbedingungen (15) bis (17) und den Anfangsbedingungen (21) sind wir auf Differen- 
zenrechnung angewiesen. Bei dem hier in Anwendung kommenden Fortsetzungsverfahren berech- 
nen wir, mit Start von den Anfangsprofilen (21), schrittweise das Mischungsfeld. Wir begniigen uns 
dabei, wie schon in der Einleitung erwahnt, fiir die Berechnung des Kernbereiches auf drei Schritte 
fiir das Geschwindigkeitsfeld, wahrend man fiir das Temperaturfeld bereits in zweien der hier voll- 
zogenen Schritte den Kernbereich durchmessen hat. 

Zur Charakteristik des Verfahrens ist zu sagen, daB der erste Schritt die Lésung eines Systems 
gewohnlicher Differentialgleichungen erfordert, die wir mittels Iteration gewinnen. Die nachsten 
Schritte erfordern dann nur noch Quadraturen. Allerdings steigert sich mit jedem weiteren Schritt 
der Formalismus, bedingt durch die notwendige Bildung immer héherer Ableitungen der Ausgangs- 
funktionen, ganz erheblich. 

Was die fiir den Erhalt glatt verlaufender Kurven der Verteilungsfunktionen in grafischer Auf- 
tragung erforderliche numerische Genauigkeit betrifft, so hat sich ergeben, daB die Anspriiche mit 
jedem weiteren Fortsetzungsschritt steigen; des weiteren wachsen die Anspriiche mit gréBer 


1Q. Pabst, a. a. O. 
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werdendem Geschwindigkeitsparameter (Uy—w)/Uo, bedingt durch die schwacher werdende Konver- 
genz} der Profile (21) fiir y > ov. 

Im folgenden berechnen wir die Geschwindigkeits- und Temperaturfelder im Kernbereich fiir die 
Geschwindigkeitsparameter (Uy—U,)/Uy © 0, = 0,25, 0,5, 0,75 und die Temperaturparameter 
OG T= Oaks 2 

a) Erster Fortsetzungsschritt. Fiir den ersten Schritt benétigen wir die Ableitungen 
(Ap/0E)e— 0 baw. (Az/05)¢— 0 baw. (Ap/0&)<— 0, bezeichnet mit Dy baw. Xq baw. oe 

Aus den Gleichungen (11), (12) lassen sich ®), Xq nicht direkt ermitteln, da diese Gleichungen 
sich fiir € = 0 auf die Gleichungen der Anfangsprofile (21) reduzieren. Durch Differentiieren von (11) 
nach £, wobei wir die Existenz und die Vertauschbarkeit der auftretenden Ableitungen als gesichert 
unterstellen, erhalten wir jedoch unter Beachtung von (18) folgende Bestimmungsgleichung fiir Do: 


” li 4 1 Ag pues 
US ®,(—F) I Po Pt ee ¥) cae D(i + sT) : ! 


2 (9904 ie] = 0. 3) 


Dabei bezeichnen im folgenden die Striche Ableitungen nach y. Gema® den Randbedingungen 
(15) ist von der Lésung asymptotisches Verschwinden zu verlangen, also 


@,—> 0 rapes a x (24) 
Entsprechend erhalten wir durch Differentiation von (12) nach unter Beachtung von (19) fiir X 
die Gleichung 


ASE Ge +4 F Gre Dy — ¥) 252i eee ) 4 7 


E\ wv 1+ (0/T;) xo 
ee HLin— U, 
= ( + a X, = Oe 
mit den aus (16) folgenden Randbedingungen 
x0 Oe Ges a at (26) 
SchlieBlich erhalten wir durch Differentiieren von (13) nach & 
Uy — U, 7 1 " 
a (D,— y D,) a ie tr Yo = 0, 
woraus unter Beachtung der sich aus (17) ergebenden Randbedingung 
Y->0 fir yr—oco (27) 
durch Integration 
y 
: = ¥ 
Po = — | yody + [y By —2 f By dy] (28) 
4 © — co 


folgt. 

Zur Lésung des Gleichungssystems (23), (25) unter den Randbedingungen (24), (26) haben wir 
den Weg der Iteration gewahlt. Wir bringen zu diesem Zweck (23), (25) zunachst auf eine zweck- 
maBige Form. 

Dividiert man Gleichung (23) durch gy), so erhalt man mit 


1 i ys ; Jane iWé Dp’ 7 
P, D, a) ae e 
Po Po 


durch zweimalige formale Integration von (23) unter Beachtung von (28) die Gleichung 


“0 


of & 
ny pee (F/T) Xo 
®o = [rs es ET + (O/T) % > 
—oco \ 
ay 
f p 
y Up —Uy, [ ‘Ug —U u,\ ® 
2 | voas ea | Py dx Py 10,4 a 7 \ Br ay + Fi (Yo —}), 29) 
0 eS. —0o 


1Vgl. Arbeit II, S. 14, FuBnote. 
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die in ihrer Formulierung bereits die erste der Randbedingungen (24) erfiillt. Die Integrations- 


konstante F’, ist so zu bestimmen, daB noch die zweite der Randbedingungen (24) Dy > 0 fiir y + oo 
erfiillt wird. 


Entsprechend erhalt man aus (25) die Gleichung 


af 
f O/T.) X, 
ae (9o/ Ty) Xo 
3 fo “T+ oT) Xe 
Sauces p 
a: Uy — Uy 2) [Sosa TEN OG 
al | Wo dx + 2 ‘A | D, dx ie tp ot ae dB\ dy + G, (yy — 1) (30) 


0 le, 2} == OX 


unter formaler Erfiillung der ersten der Randbedingungen (26). Die Integrationskonstante G, ist 
wiederum so zu bestimmen, dab X,— 0 fiir y— oo. 


Ww 


Abb. 2. Abb, 3. 


Die in (29) und (30) auftretenden Funktionen gy, ~, und ¥, 79 sowie y, haben wir, wie bereits 
erwahbnt, in der Arbeit IIT berechnet. 

Drr Erfolg eines iterativen Rechenganges — Konvergenz und Arbeitsaufwand — ist natiirlich 
weitgehend abhangig von der Wahl der Eingangsfunktionen OY, X(. 

In dieser Hinsicht bietet einen wichtigen Hinweis der Fall (ug — u,)/uy © 0. Hier haben wir, 
wie man am einfachsten durch Einsetzen in die Gleichungen (23), (25) verifiziert, die exakten L6- 
sungen+ 


a 5 (31) 
xX,= are fiir ae beliebig. 


1 
Das fiihrte uns allgemein zu dem Ansatz fiir die Kingangsfunktionen 
PY =A q> XY = By, 
wobei die Koeffizienten A, B durch Eingang mit obigem Ansatz in (29) baw. (30) und Kollokation 
an der Po See DZWe x0 max entsprechenden Stelle y bestimmt wurden. 

Mit den so bestimmten Eingangsfunktionen fihrte der Iterationsgang nach (29), (30), wobei 
zwischen den einzelnen Gangen noch gemittelt wurde, in den Fallen (uy — u,)/u, 0, = 0,25, 0,50 
nach 3—4 Iterationen zum Ziel, d. h. die Rechnung kam in der grafischen Auftragung eines DIN A 4- 
Blattes zum Stillstand. Fiir (uj — u,)/uy = 9,75 wurden die Eingangsfunktionen durch Extrapo- 
lation gewonnen; hier benétigte man 4 bis 5 Iterationen. 

Wir kénnen an dieser Stelle die berechneten Funktionen nur im Auszug wiedergeben. Abb. 2 
und 3 zeigen fiir ),/T, = 2 als Beispiel den Gang der Funktionen ®,, X, mit dem Parameter 


Any Mk co 


Vn Je (80%) + BP 


mit 


1 fs 
1 Hier ist (vgl. Ing.-Arch. 20 (1952), S. 73) ~ = — Vx (Pol T1) x0 =— 
yo y/VE und der Fehlerintegralfunktion ®(y*). 
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(uy) — u,)/uy. Die in Erscheinung tretende Rechtsversetzung der Kurven war zu erwarten, da sie 
der Rechtsversetzung der Anfangsprofile q, 7p mit dem Parameter (uy — u,)/u) entspricht (vgl. Ar- 
beit IT). 

Abb. 4 und 5 wiederum fiir (u, — u,)/u) = 0,5 als Beispiel zeigen eine mit dem Parameter @)/ T: 

0 V)/Uo iid § ; of 4 
gehende Linksversetzung der Kurven, die dem analogen Verhalten der Anfangsprofile entspricht 
vel. Arbeit IT). Des weiteren weist der mit 0,/T, wachsende Betrag der Funktionen auf die mit 
S o/ 4 8 

gréBerer Temperaturdifferenz zu erwartende starkere Kriimmung der Isotachen und Isothermen ~ 
hin. 

Wir kénnen nunmehr, ausgehend von den Anfangsprofilen gy und 7», den ersten Fortsetzungs- 
schritt durchfiihren und erhalten mit der Schrittweite a approximativ die Profile 


G1 =% +aD, Xy =o +aXo-. (32) 


Das Profil der Querkomponente der Geschwindigkeit p, an der Stelle € = a ermitteln wir nun nicht 
etwa in entsprechender Weise mittels Y. Vielmehr ist, wie wir unten zeigen, bei vorliegenden Pro- 
filen y,, 7, die Funktion y, eindeutig durch diese bestimmt und mittels einer expliziten Formel be- 
rechenbar, siehe (38). Wir kénnen uns also in bezug auf die approximative Konstruktion auf die 
Profile p, und y, beschranken. 


LX, 


Abb. 4. Abb. 5. 


In der vorliegenden Rechnung haben wir die Schrittweite a in Hinsicht auf iiberschlaglich er- 
mittelte Kernlangen (§ ~ 0,45) mit a = 0,2 bemessen. 


+b) Weitere Fortsetzungsschritte. Nach Vollzug der ersten Fortsetzung lassen sich die 
weiteren Fortsetzungen mittels expliziter Formeln, die bei der numerischen Auswertung nur 
Quadraturen erfordern, gewinnen. 


Zunachst stellen wir fest, da® wir mit der aus (13) durch Integration zu gewinnenden Formel 


bh 
up I i i 0p op 
Me Up ee Se | A+Enlyy—8 Ely, (33) 


—i/s 


welche die Randbedingung (17) erfiillt, aus Gleichung (11) die Funktion y eliminieren kénnen. 


Liegen nun an einer Stelle € += 0 die Profile vor, so erhal : F 
: > , t durch E 
in (11) und Umordnung die Gleichung oe QUE a SGD Er aan nr 


itd , Bo T, 4 
(1 +n p a ges eee ! TF 


; Op\ /Uy—u u Uy — 
+2fal bso9(oy £32) nee a Seay [a Healey eBay =9 
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oder 
1 ee (| o(E (A/T) x E 
C+eo)e’ +0'|-@ +s try 
ap 
ere dy 
Tay | Homme ies ye et Veo f a(vw ee a) ee 
Re! yas opel ; =0, 
ere ae = 


aus der durch Integration, bei Wiedereinsetzen von y gemaB (33) und Erfiillung der Randbedingung 
(17) folgt 


y 
Up—_ Uy ! uy [ , (o/ T;) Va é 
si —(E+1 + = 
ae glrugsc uy Ye Ug ( oe 1 | Sa T+ GIT YA rae 
y oti, l+éy = z,\3 dy. (34) 
Hg 4) 
Uo Uo 


—Yé 
Die den Profilen y, y entsprechende Funktion @g/0& gewinnen wir dann fiir & ++ 0 leicht aus (11) 


Zu 
wr ‘jo Commer uy E (9o/T;) x’ | S 
Beles Pee i OF UTE wityx * IFFY 


op 


ag ine - 26 (es Pe ss) 5) 
*—? 
Entsprechend erhalten wir aus (12) 
a , 2 Up Uy Wy (Fo/ T;) vis & 
2 x = 9 
ox Z : 4 Fall Ug etitly | 1 + (8)/T,) x Fy 36 
0a- =. 2 i fa 2 (to , Uy 5 (36) 
(=) (See ta 
Die Randbedingungen lauten gemaf (15), (16) 

ay ay 0 »  (y—>—I/é 

ME? OE | an es ie 5 (37) 


' Die Erfiillung der ersteren Bedingung erfordert beim Fortsetzungsverfahren besonderes Augen- 
merk; siehe 2. und besonders 3. Schritt. 

Wir merken an, dafs fiir den Fall ruhender AuBenluft (u, = 0) in (35), (36) fiir y > oo nicht nur 
die Zahler, sondern auch die Nenner Null werden. Das erklart schon, da dieser Fall eine besonders 
groBe Rechengenauigkeit erfordert. 

Es bleibt noch zu sagen, daf infolge der in (34) bis (36) vorkommenden Ableitungen g’, p’’ sowie 
7x, 7’, die im Fortsetzungsverfahren explizit mit den Ausgangsfunktionen qo, 7) und Do, Xo (erster 
Schritt) verbunden sind, jeder weitere Schritt (vgl. zweiter und dritter Schritt) die Bildung immer 
héherer Ableitungen dieser Funktionen erfordert, so daB der zu erbringende Rechenaufwand sich 
mit jedem weiteren Schritt ganz erheblich steigert. 

c) Zweiter Fortsetzungsschritt. Die fiir den zweiten Schritt benétigten Funktionen (dp/cé), 
bzw. (67/0&),, bezeichnet mit ®, bzw. X,, berechnen wir, ausgehend von den im ersten Schritt 
konstruierten Funktionen q,, 7,, nach den Gleichungen (34) bis (36). 

Gleichung (34) liefert 


eee 

Ug Uy oy 1k , een (9,/T,) val \ a 

_ 1 pans uo mt i] (+ ay) |e | 7 | = T+ @lT) x lee ay d (38) 
aes 2 Up Le ay. j Up — Uy e 7 y: 
mn 
p ; 
—I1fa 
Mit der Geweinnung von y, haben wir zunachst den ersten Schritt abgeschlossen. 
Die Gleichungen (35) und (36) ergeben dann 
" / Bhi ey 7 (Po/ Ts) x4 g 
v1 +0 ( on eisin ae v, (E+ oh Pa 
aio ts)x, la ay 
i= - ° _ Co ag (39) 
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und ; 
re 2 (F/T) x1, ent 
tat za 1+ (JT) x 1+ ay a 
xy a i / 2 u, The my) mes : ( ) 
oy Uy 
(eet a) 
Hierbei ist nach (32) 
P1 a Po ar D, ? Py aoe Po ar PD, ? (41) 
ta =o + 4Xy; ti =X + aX. 


Hinsichtlich der Randbedingungen (37) ist anzumerken, daf} die erstere Bedingung infolge des 
asymptotischen Charakters der Funktionen yp, @y und 7, Xo nicht exakt erfiillt ist; doch ist hier 
bei einem Zahlenwert von —l/a = —- 5 der Fehler noch sehr klein. rs 

Die Funktionen qp, gy, sowie Yo, %4 sind aus Arbeit IT bekannt. yy bzw. 7% gewinnen wir dann 
aus (18) bzw. (19). Die Funktionen ®,, ®) sowie Xo, Xj sind aus der im 1. Schritt erfolgten Be- 
rechnung nach (29) und (30) bekannt. ®f baw. X{ berechnen wir dann aus (23) baw. (25). 


Abb. 6. Abb. 7. 


Die fiir die schon genannten Parameter (uy — u,)/U9, O/T, berechneten Funktionen @,, X, 
kénnen wir wieder nur im Auszug bringen. Abb. 6, 7 fiir 0,/T,=2 baw. Abb. 8, 9 fiir (uy—u,)/uy 
= 0,5 zeigen denselben Gang mit (uy — u,)/uy baw. J/T, wie die entsprechenden Abb. 2, 3 
bzw. Abb. 4, 5 des ersten Schrittes. 

Die Werte des zweiten Schrittes haben wir in Anbetracht der zu erwartenden Verkiirzungen 
der Kernlange mit wachsenden Temperaturparameter %)/T, zu 


b=0,15 fir 9,/T,=0, 
0,125 
0,10 2 


2 


bemessen. 


Damit befinden wir uns zwar, wie man aus der Auftragung des sich aus Por Xo (Arbeit IT), y,, vy; 
nach (32) und %,, 7, nach (42) ergebenden Isotachen- und Isothermenfeldes ersieht, zwar noch inner- 
halb des Kernbereiches des Geschwindigkeitsfeldes, wenn auch schon in der Nahe des Kernendes ; 
gelangen aber mit diesem Schritt bereits itber den Kernbereich des Temperaturfeldes hinaus (vel. 
Abschn. IV, Tabelle 2 u. 3), so daB fiir das Temperaturfeld bereits mit dem zweiten Schritt die 
Aufgabe beendet ist. 


Wir erhalten in Fortsetzung nunmehr die Profile 


Me | 


i (42) 
te = 1 +b6X, =y% +aX,+bX,. 
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a) Dritter Fortsetzungsschritt. Der dritte Schritt verliuft ganz analog dem zweiten 
Schritt. 

Obwohl wir mit dem zweiten Schritt bereits iiber das Kernende des Temperaturfeldes hinaus 
gelangt sind, uns also schon in dem Bereich befinden, der durch standigen Abfall der Temperatur 
langs der Strahlachse gekennzeichnet ist, kinnen wir die Konstruktion des Temperaturfeldes noch 
weiter fortsetzen, solange wir uns nur im Kernbereich des Geschwindigkeitsfeldes befinden; denn 
erst nach Erreichen des Kernendes des Geschwindigkeitsfeldes tritt eine Anderung in unserem 
Rechenschema dadurch ein, da infolge des nun einsetzenden Abfalls der Geschwindigkeit langs 
der Strahlachse die scheinbare kinematische Zahigkeit nunmehr beschrieben wird durch 


e(x) = 4, by (x) (w4(x) — uy) , (43) 
wo u,(x) Geschwindigkeit auf der Strahlachse. 
Fiir den dritten Schritt miissen wir die Funktion (Cp/CE)a +» bezeichnet mit ®,, und, wenn 


wir das Temperaturfeld noch weiter fortsetzen wollen, die Funktion (07/0&)a 4.5, bezeichnet mit X,, 
berechnen. 


+2, eo Ht 5 Xx; 
0 
Lay 
Abb. 8. Abb. 9. 


Die sich aus (34) bis (36) ergebenden Gleichungen lauten jetzt 


0 
= ald ln Ge a+b | 
Uo pel ea " A= Bs ] ( Oo} 1) A2 | 
‘ ied jis Po + =a + (0+ d)x1fob + 94] (Ee TE GIT) ee rFeroHI| . 
= —"5 Up l1+(a+b)y re + By 
Ug Uf) 
—1f(a +b) (44) 
womit wir zunachst den zweiten Schritt abgeschlossen hatten. 
" , Up Uy Uy | (/T;) Xe ; a+b 
a (ee mt ity v. OF OTE OT | 1+ @+ Oy ne 
a Cy iL ul ; 
2(a4 6) te a 
" , 2 39 (80/7) Xe a Te b 
Peale T+ (O/T) to | 1+ (@ + B)y na 
aa : 


pot) 9, +2) 


Uo 


Hierbei ist nach (42) 
Poe r Ol Bie), fag t bor = of 400% +b OF: bs 
Hite XM Hy teXt+bX, Hm tox t+aG +bXy. 
Gegeniiber dem zweiten Schritt bendtigen wir also zusatzlich die Funktion OO, und Aja 
25 
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Die Funktion @', OY baw. Xj, X7 sind durch sukzessives Differentieren aus (39) baw. (40) 
zu erhalten, was wiederum die Berechnung der Ableitungen 9’, gf) bzw. 71> xy ) und damit 
nach (41) die Berechnung der Ableitungen Po > PAY) und Ves xh), die saben near ere 
tieren von (18), (19) erhaltbar sind, sowie die Berechnung der Ableitungen Dy, YY) und Xy, X} 
erfordert, die durch sukzessives Differentieren von (23) baw. (25) gewonnen werden kénnen. 

Man sieht schon, daB® der Rechenaufwand gegeniiber dem zweiten Schritt ganz erheblich an- 
gewachsen ist. ; 

Eine Bemerkung erfordert noch die erstere der Randbedingungen (37). Mit nunmehr 
—1/é =—1/(a + b) (w — 2,9 fir B)/T, = 0 2. B.) ist fir @, die Abweichung von der au er- 
fiillenden Randbedingung schon merkbar, so daB die Anbringung einer Korrektur erforderlich wird. 
Fiir X, andererseits ist gemaB dem Umstand, daf die Temperatur langs der Strahlachse bereits 
abfallt, die Bedingung (37) zu ersetzen durch die Randbedingung 


a o 
Cle mires STUN A/a) ay i 


Man erhalt diese Bedingung aus (12) durch Grenziibergang, wenn man beachtet, da auf der 
Strahlachse 0y/dy = 0. 


Abb. 10. Abb, 11. 


Nach dem angegebenen Schema haben wir fiir die schon genannten Parameterwerte die Funk- 
tionen ®, und in einigen Fallen auch die Funktion X, berechnet. Es ist zu bemerken, daB die durch 
den ersten Schritt festgelegte Rechengenauigkeit (glatter Kurvenverlauf von ®, baw. X, in der 
Auftragung eines DIN A4-Blattes) nunmehr in der dritten Fortsetzung nicht mehr ganz aus- 
reichend war. Die sich im unregelmafigen Verlauf der Kurven — insbes. fiir X, — zu y— oo hin 
zeigende Ungenauigkeit der Rechnung (vgl. z.B. Abb.12) trat bei gréBerem (uy — u,)/uy und 
kleinerem #,/T, starker hervor. Abb. 10, 11 fiir 0)/T,=2 und Abb. 12, 13 fiir (uy — u,)/uy = 0,5 
zeigen wieder die berechneten Funktionen ®,, X, in Auswahl, wobei die letzteren nunmehr die 
Randbedingung (48) erfiillen. 


Die Werte des dritten Schrittes haben wir mit c=0,1 bemessen. Wir erhalten dann in Fort- 
setzung die Profile 


Ps =%, +e®, =p +aG+bO,+cO, (49) 
und, falls wir die Konstruktion des Temperaturfeldes noch weiter fortsetzen wollen, 
Xs =o +eX,=H% +aX,+bX,+c¢O,. (50) 


Mit diesem Schritt gelangen wir bereits tiber das Kernende des Geschwindigkeitsfeldes hinaus. 
Da jedoch der Abfall der Geschwindigkeit, wie auch der Temperatur, langs der Strahlachse zunachst 
ganz allmahlich erfolgt (vgl. Abb.17), besteht hier noch keine Notwendigkeit einer Anderung 
unseres Rechenschemas gemaB (43). 

In Abb. 14 bzw. 15 haben wir schlieBlich noch fiir (uy — u,)/uy = 0,25, 8/T, ~ 0 als Muster- 
beispiel die Funktionen ®,, ®,, ®, baw. Xp, X,, X, nebeneinander dargestellt. 
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Abb. 12. Abb. 13. 


Abb. 14. Abb, 15. 
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4. Beschreibung des Geschwindigkeits- und Temperaturfeldes im Kernbereich ‘ Die in Abschn. 3 
durchgefiihrte Berechnung ergab die Verteilungsfunktionen der Geschwindigkeit bzw. Temperatur 
Por Pir Por Ps bzw. Yoo Xr» Xo Xs in den Schnitten & = 0, a, a :9 b, a+b +c des Kernbereiches; 
vgl. die in Abb. 16 dargestellten Profile im Fall des Musterbeispiels (uy — U,)/ Uo == (Li. Di Ty 0. 

Wir beschreiben nunmehr das Geschwindigkeits- bzw. Temperaturfeld im Kernbereich mittels 
des Verlaufs der Isotachen bzw. Isothermen: r/r, tiber x/r, fiir g = konst. baw. y = konst. Nach 
der in Abschn. II gegebenen Definition y = (r/ry — 1)/&, € = (1/0) 4/19 haben wir 

St =1t+éy, ~ =a. (51) 
Yo Yo 
GemaB unserer in Abschn. 2 getroffenen Annahme o(x) = konst. fiir das einzelne Modell ent- 
nehmen wir o der Arbeit II (o = a(0)): Tabelle 1. 


Tabelle 1 
iam ea) Gages Re aah. pee 2050 i ae Te 
Uo 
9,[T, = 0 | 1 | 2 
Come 26,28 26,28 | 26,75 | o | 26,23 | 25,93 
Uo 
Geschwindigkelt 
------ Jemperatur 
PX ih 


aot = 05 “ot 60 Upn tly _ pas 


Vot/g 


| | | 
0 70 20 30 40 50 
Abb. 17. Geschwindigkeit und Temperatur lings Strahlachse. 


Dabei hat, wie in Abschn. 2 erértert wurde, der empirische Koeffizient x, den Zahlenwert 
x, = 0,0082. Fiir den Fall, daB %, den Charakter einer Modellkonstanten haben sollte, waren 
auch andere Zahlenwerte denkbar; die Werte der Tabelle 1 waren dann gemaB der Relation 


* 
Cage: 


o ut 
zu modifizieren.? 
Das Ergebnis unserer Rechnung enthalten die Tabellen 2 und 3. Die effektiv vorhandenen 
Geschwindigkeiten bzw. Temperaturen ergeben sich gema den Formeln 


u=(uy—u,)y+u, baw. T=(T,—T,)7+T,. 

Auf eine zahlenmaBige Angabe der Kernlangen haben wir verzichtet, da bei dem ganz allmahlich 
einsetzenden Abfall der Geschwindigkeit und Temperatur langs der Strahlenachse (vgl. Abb.17), 
derartige Angaben wenig sinnvoll erscheinen. Wir haben statt dessen den Verlauf der Geschwindig- 
keit und Temperatur langs der Strahlachse in Abb. 17 dargestellt. 

Abb. 18% zeigt die sich ergebenden Isotachen- und Isothermenfelder im Beispiel der Faille 


(uj — u,)/a5 = 0,5; 3/7, =O, 1) 2. 
Aus Abb. 17 und 18 ergibt sich als Temperatureffekt eine longitudinale und laterale Verkiirzung 


des Mischungsfeldes; diese Verkiirzung entspricht der aus Arbeit II folgenden Drehung des Winkel- 
raumes der Mischung, mit Scheitel in der Diisenkante, zur Strahlachse hin. 


1 Die Durchfiihrung der schwierigen numerischen Rechnungen verdanke ich meinen Mitarbeitern Herrn 
Herbert Moch und Frau Lilli Schulze. 


® Das ergibt die aus (10) folgende Gleichung: 
1 
— : 
2 (0 cy) (Up — Uy)/Uty 


8 In Abb. 18 haben wir, um den Kontrast zwischen Geschwindigkeits- und Temperaturfeld besser hervor- 
zuheben, nur den Kernbereich des Temperaturfeldes wiedergegeben. 
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Doch ist der Temperatureffekt hier nicht sehr erheblich. Von groéBerer Bedeutung erscheint 
dagegen die wesentlich andere Struktur des Temperaturfeldes gegeniiber dem Geschwindigkeitsfeld, 
die durch den sich aus dem Experiment ergebenden Wert E=2 [vgl. (22)] fiir das Ubertragungs- 
verhaltnis bedingt ist und sich in einer gréReren Breite sowie geringeren Tiefe des Temperaturfeldes 
pragnant dartut. 
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Abb. 18. Feld der Isotachen und Isothermen. 
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v 
In Tabelle 4 haben wir die Zahlenwerte der Querkomponente der Strémung y = one angegeben, 
aus der v=u,y/o folgt. Abb. 19 zeigt die erhaltene Querkomponentenverteilung im Beispiel der 


Falle (uj — u,)/Up = 0,5; Bo/T, = 9, 1, 2. Das im Falle bewegter AuBenluft bei gréBeren Tem- 
peraturdifferenzen in Erscheinung tretende Phanomen des Abstroms vom Mischungsfeld haben wir 


bereits in Arbeit II diskutiert. 


Tabelle 4. Querkomponente der Geschwindigkeit 


Mis 1 | 0,01 | 0,03 | 0,04 | 0,05 | 0,07 0,00 | 0,00 | 0,01 | 0,01 | 0,01 | 0,01 
Eso 108 0.08 | 0,17 | 0,21 | 0,31 | 0,38 0.01 | 0,02 | 0,03 | 0,04 | 0,06 | 0,07 
M6025 1.0.30) |.0.67 |.0,78 | 7, || heal 0.05 | 0,07 | 0,14 | 0,16 | 0,23 | 0,26 
12} 076 | 085 | 1.77 | 191 | 2,80 | 2,94 | 0,17] 0,19 | 0,37 | 0,42 | 0,57 | 0,61 
—o'3| 161 | 171 | 3:24 | 3,28 | 4,62 | 4,56 | 0,36] 0,38 | 0,70 | 0,73 | 0,98 | 0,98 
—o'4} 2’s0 | 2252 | 433 | 416 | 5,66 | 5,35 | 0,56] 0,57] 0,96 | 0,96 | 1,24 | 1,19 
o | 2'31 | 2274 | 450 | 4:16 | 5,67 | 5,18 | 0,63 | 0,63 | 1,02 | 0,97 | 1,27 | 1,18 

0,4] 2:30 | 2.20 | 3,97 | 3,53 | 5,11 | 4,51 0.51 | 0,50 | 0,89 | 0,82 | 1,14 | 1,02 

03} 137 | 129 | 3.27 | 2:80 | 4,47 | 3,81 | 0,27] 0,27} 0,70 | 0,62 | 0,96 | 0,84 

121 059 | 0,56 | 2,74 | 2.25 | 3,98 | 3,29 |-+0,04 |+0,06 | 0,53 | 0,46 | 0,82 | 0,70 

16] 018 | 017 | 2.42 | 1.91 | 3,65 | 2,94 |—0,10 |—0,07 | 0,42 | 0,35 | 0,72 | 0,60 

20} 004 | 0,04 | 223 | Ll | 3,41 | 2,69 |—0,17 |—0,13 | 0,36 | 0,29 | 0,66 | 0,53 

241 0,01 | 0,01 | 210 | 1,58 | 3,22 | 2,49 |—0,18 |—0,14 | 0,33 | 0,26 | 0,62 | 0,49 

2.8 1:99 | 1.47 | 305 | 2,33 |—0,18 |—0,14 | 0,31 | 0,24 | 0,58 | 0,45 

3.9 139 | 137 | 2:90 | 2.19 |—0,17 |—0,13 | 0,29 | 0,22 | 0,55 | 0,42 

5G 31 | 1.29 | 2:77 | 2:06 |—0,16 |—0,12 | 0,28 | 0,21 | 0,52 | 0,40 

4.0 1.72 | 1.22 | 264 | 1,95 |—0,16 |—0,12 | 0,27 | 0,20 | 0,50 | 0,38 

4.4 1.65, | 1,15 | 2653 | 1,85 |220,15-1=0,11 | 0,26, |, 0,198 0,400 ee 

Lomas 
Uo 
O/T, = O/T, = 1 O/T; = 2 

A eta oe BAGH 2 S08 || e2a 

10 y, 10 yw, 10 y, 10 wy, 10 y, 10 w, 

2,4 0,00 | 0,01 | 0,01 | 0,01 | 0,02 0,00 | 0,00 | 0,01 | 0,01) 0,02 | 0,02 
—2'0! 0,02 | 0,03 | 0,05 | 0,07 | 0,09 | 0,11 0,03 | 0,04] 0,06] 0,08) 0,10} 0,12 
—16| 0,10] 0,12 | 0,22 | 0,27 | 0,35 | 0,41 0,13 | 0,16] 0,26 | 0,311 0,39 | 0,46 
1:21 0.29 | 0.33 | 0,62 | 0,70 | 0,91 | 1,00 | 0,39] 0,44) 0.731) 0,834) 1,041) 116 
—9.8| 0,63 | 0,68 | 1.19 -| 1,18 | 1,62 | 1,69 | 0:84 | =0,92)- 1,449) 155))eoasmumotna 
p41 10:99.) 1103-4) oti s4p eo 14 loa 1,33 | 1,39] 2,08| 2.141 2,651 2.69 
0 V13 > 1,34 Teo | 1:63 31-9.03 <2,14 1,53.| “1,55°| 2,994) 2.07 oe ead: 

0,41 0.92 1 0.01 |) 1,578 | 1.38 |-1,99. | 1:85 1,24. | ©1,24 \ 1,98 | 1.91 1) 9.51 Jeeoead 

0,8|-+0,44 |+10,46 | 1,17 |.0,98 | Le2 | 1,47 |+-0,54 | 20,59 |-.1,36) 1.30008 Loom aeine 

1,2 | —0,09 |—0,02 |-0,78 | 0,62. | 1,28 | 1,13 |—0,30 |—0,15 | 0.70) 0.6900 01059) 95 

1,6 | —0,49 |—0,37 | 0,50 | 0,37 | 1,03 | 0,88 |—1,02 |—0,77 |-L0,13 |-+0.17 | -0,85 | 0,78 

2,0|—0,71 |—0;57 | 0,33 | 0,23 | 0,87 | 0,72 |—1:53 |—1;22:/=-0.07 011900 0.48 omnes 

2,4|—0,79 |—0,64 | 0,25 | 0,16 | 0,77 | 0,62 |——1,81 )—1,49 |~0,52, |--20.450) 0.03) Og 

2.8 |—0,80 |—0,64 | 0,21 | 0,13 | 0,71 | 0,56 1,93 | —1,60 | —0,66 |—0,61 | +0,07 | --0,05 

3,2|—0,78 |—0,62 | 0,19 | 0,11 | 0,67 | 0,52 |—1,95 |—1,60 |—0,72 |—0,65 |—0,01 |—0.03 

3.6 |—0.75 |—0,59 | 0.17 | 0.10 | 0,63. | 0.49 1,92 |—1,56 | —0,74 | —0,66 |—0,06 | —0,07 

4,0]—0,72 |—0,56 | 0,16 | 0,10 | 0,60 | 0,46 |—1,86 |—1,49 |—0.73 | —0,64 |—0,08 |—0,09 

4,4|—0,69 |—0,53 | 0,16 | 0,09 | 0,58 | 0,44 |—1,79 |—1.42 |—0.71 |—0,62 | —0,08 | 0.10 

4,8|—0,66 |—0,50 | 0,15 | 0,09 | 0,55 | 0,42 |—1,72 |—1.35 |—0,69 |—0,60 |—0.08 | 0,09 

5,2 |—0,63 |—0,47 | 0,15 | 0,08 | 0.53 | 0.40 1,66 |—1,28 | —0,66 | —0,57 |—0,08 | 0,09 

5,6 |—0,61 |—0,45 0,08 | 0,51 | 0,38 1,60 |—1,22 | —0,63 | —0,54 |—0,08 | —0,09 


Wir bemerken noch, da8 die hier in drei Fortsetzungsschritten erhaltenen Verteilungen mit nur 
geringer Abweichung bereits mittels zweier Schritte erhalten werden kénnen, wie Vergleichs- 
rechnungen ergaben. Auf Grund dieser Feststellung haben wir in den Fallen, in denen die be- 


XXVI. Band 1958 = W. Szablewski: Turbulente Ausbreitung runder Hei®luftstrahlen in bewegter Luft 375 


rechnete Funktion X, einen sich durch gréfere Schwankungen erweisenden unbefriedigenden Ver- 
lauf zeigte, das Temperaturfeld in der gesamten Langserstreckung des Kernbereiches des Geschwin- 
digkeitsfeldes mittels zweier Schritte berechnet (Falle (uy — u,)/uy) = 0,5; %/T,=0 und 
eo UU 0.15; 9/1 0,1; 2). 


5. Vergleich mit Messungen. Zum Vergleich zichen wir, wie in Arbeit II, die von 0. Pabst! 
an HeiSluftstrahlen in bewegter Luft durchgefiihrten Messungen heran. Der aus einer Diise mit 
10 mm Durchmesser tretende Strahl hatte eine Geschwindigkeit von 2400 m/sec und eine Tem- 
peratur von ~ 400° Celsius. Dieser HeifSluftstrahl trat in einen Freistrahl von 80 mm Durchmesser 
mit Geschwindigkeiten von ~ 180, 100 und 26 m/sec sowie normaler Temperatur. Die dimensions- 
losen Daten betragen 


Wir betrachten hier nur die beiden ersteren 
Falle, wahrend der letzte Fall bei der Untersuchung 
der Ausbreitung in ruhender Luft (vgl. Einleitung) 
herangezogen werden soll. 

Die in Arbeit II erhaltenen Profile qo, v7) in 
Strahlenkoordinaten (die Anfangsprofile & =0 der 
vorliegenden Rechnung) hatten in den oben ge- 
nannten Fallen befriedigende Ubereinstimmung 
mit den Messungen sowohl hinsichtlich der Misch- 
breite als auch hinsichtlich des Verhaltnisses von Ay 
Auen- zu Innenwinkel, diese von der durch den 
Diisenrand in Richtung der Strahlachse gelegten 
Geraden gemessen, ergeben. Dabei hatten wir 
beziiglich der Geschwindigkeitsverteilungen dem 
durch Grenzschichteinfliisse an der Diisenwan- 
dung bedingten Impulsverlust durch Einfiithrung 
eines effektiven Diisenradius r* gemaB der Glei- 
chung 


* 
feouu—wu) r dr = Qo Uy (up — wy) “2 
6 

Rechnung getragen. 

Nach dem Gesagten ist, da sich sowohl nach 
den Messungen von O. Pabst wie nach unserer 
Theorie (vgl. Abb. 16) fiir den Kernbereich nur 
eine geringe Verschiebung der Profile in Strahlen- 
koordinaten vom Anfangsprofil ergibt, von vorn- 
herein gute Ubereinstimmung der Theorie mit dem 
Experiment hinsichtlich des Verlaufs der Grenzen 
der Mischfelder im Kernbereich zu erwarten. 

In Abb. 20 haben wir fiir die genannten Fille die theoretischen und experimentellen Felder 
der Isotachen und Isothermen fiir den Kernbereich aufgetragen. Die dem Bericht von O. Pabst 
entnommenen ,,experimentellen“’ Punkte wurden dort durch Schnitte der aus den einzelnen Meb- 
punkten konstruierten Profile erhalten und sind demgema48 mit einer gewissen Streuung behaftet. 
Punkte, die nach der im Bericht selbst-vorgenommenen Mittelung (mittlere Kurven) und im unten 
angegebenen Fall starker streuten, wurden hier fortgelassen. Eine Korrektur wurde an einer der 
angegebenen Isotachen vorgenommen.* 

Die diisennahen Schnitte, die noch merkbar den Einflu8B der Grenzschicht an der Diisenwandung 
zeigen (vgl. Arbeit II), wurden aufer acht gelassen. Entsprechend der Einfiihrung des effektiven 


1Q, Pabst, a.a. O. 

2Es wurde im Fall (up — u,)/Up = 0,53; F/T, = 1,3 der Punkt x/2r9 = 8, 1/27) = 0,58 der Isotache 
u/uy = 0,6 an Hand der MeBdiagramme zu r/2 ry = 0,54 korrigiert, was einen anderen Verlauf der Isotache 
ergibt. Die Nachpriifung wurde veranlaBt durch den unwahrscheinlichen Verlauf der angegebenen Isotache, 
der fiir das Geschwindigkeitsprofil y(7) bei x/2 rp) = 8 aufen einen steileren Gradienten ergeben wirde als innen. 
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Abb. 21. Geschwindigkeit und Temperatur langs Strahlachse. 


XXVI. Band 1958 = W. Szablewski: Turbulente Ausbreitung runder Heifluftstrahlen in bewegter Luft ayn i 


Radius ry wurden die Punkte der experimentellen Isotachen in die Diagramme mit den Koordinaten 
x/2rf und r/2 rf eingetragen.! 

Der Vergleich zeigt die erwartete gute Ubereinstimmung der Theorie mit dem Experiment hin- 
sichtlich der Begrenzung der Mischfelder. Des weiteren wird die von der Theorie angezeigte wesent- 
lich verschiedene Struktur der Geschwindigkeitsfelder und Temperaturmischfelder — gréBere Breite 
sowie geringere Tiefe des Temperaturfeldes — bestatigt. 

Dagegen ist trotz der Ubereinstimmung der Grenzen der Mischfelder eine z. T. starkere Ab- 
weichung der experimentellen Isotachen und Isothermen im auBeren Teil der Mischfelder von den 
theoretischen Kurven festzustellen. Welche Griinde dafiir maBgebend sind, lassen wir dahingestellt. 


oi 
ii a 


Theorie 


Ug ~U, Ls Vy = e 
09 Uo = thes ) L, =13 e 
08 e Messung 0.Pabst U. u.M. 8004 (1944) 
4 
[ee ! 1 aa ee I | | i l l Wiha N uo 
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Abb, 22, Geschwindigkeit langs Strahlachse. 


In Abb. 21 ist noch der Verlauf der Geschwindigkeit und Temperatur langs der Strahlachse 
angegeben. Dabei wurden die MeBpunkte der Geschwindigkeit wieder iiber x/r* eingetragen. Die 
Ubereinstimmung von Theorie und Experiment ist hier wieder als befriedigend anzusehen. Ins- 
besondere wird die theoretische Aussage von der geringeren Tiefe des Temperaturfeldes gut be- 
statigt. 

Wiinschenswert ware noch ein Vergleich mit Messungen, die den Effekt der Temperaturdifferenz 
zwischen Strahl und umgebendem Medium herausstellen. Solche Messungen scheinen fiir bewegte 
AuBenluft nicht vorzuliegen. Wir haben jedoch in Abb. 22 den Geschwindigkeitsverlauf langs der 
Strahlachse des Falles (uy — u,)/uy = 0,53; %)/T, = 1,3 mit der des Falles (uy — u,)/uy) = 0,50, 
},/T, = 9, wobei die Messung einer Alteren Arbeit des Verf.2 entnommen wurde, verglichen. Die 
effektiven Radien betrugen in beiden Fallen rf /ryg =0,88. Der Temperatureffekt der longitudinalen 
Verkiirzung wird experimentell befriedigend wiedergegeben. 


6. Zusammenfassung. Die Mischfelder der Geschwindigkeit und Temperatur bei der turbulenten 
Ausbreitung von runden HeiSluftstrahlen in bewegter Luft werden in der Zone des Kernbereiches 
fiir verschiedene Werte der Geschwindigkeits- und Temperaturdifferenz berechnet. Grundlage 
bildet das vom Verf. in einer vorangehenden Arbeit (Ing.-Arch. 20) aufgestellte Gleichungssystem 
fiir turbulente Stromungen von Gasen stark veranderlicher Dichte. Das sich hier ergebende System 
nichtlinearer partieller Differentialgleichungen wird mittels eines Fortsetzungsverfahrens gelést. 
Ein Vergleich mit Messungen zeigt befriedigende Ubereinstimmung. 


(Eingegangen am 21. Dezember 1957.) 
Anschrift des Verfassers: Dr. W. Szablewski, Berlin W 8, Otto-Nuschke-Str. 22/23 


1 Aus denselben Griinden, wie in der vorigen FuBnote angegeben, wurde im Fall (uy — u,)/u9 = 0,74 die 


Isotherme vy = 0,2 statt 0,3 herangezogen. 
2 W. Szablewski, U. u. M. 8003 (1944). 
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Kuppeln gleicher Festigkeit 
Von H. Ziegler 


1. Problemstellung. Unterliegt eine Rotationsschale einer rotationssymmetrischen Belastung, 
gegeben durch die Radialkomponente p, und die Meridionalkomponente p, der Last je Flachen- 
einheit, dann sind die Schnittkrafte N,, N» ihres Membranspannungszustandes! durch 

d N, N, 
_- (r N,) —1r, Ng cos 9 = — Pot » Ae Aa eaap (1.1) 
1 
bestimmt. Dabei sind q, r,, Tr, und 
r =f, sing (1.2) 
geometrische Bestimmungsstiicke der Schale, die der Abb. 1 entnommen werden kénnen. Der 


Spannungszustand in der Schale ist eben und durch 
N, N. 
p= > 05 = = (1.3) 


gegeben, wenn t die Schalenstarke ist. 
Schalen gleicher Festigkeit sind durch einen Membranspannungszustand gekennzeichnet, der 
iiberall an der zulassigen Grenze liegt. Darunter wird iiblicherweise der Spannungszustand 


Oy = 95 = 7 Oo (1.4) 
verstanden, wobei o, die zulassige Grenze fiir beide Spannungen bezeichnet. Auf dem Ansatz (1.4) 
beruht beispielsweise der Tropfenbehalter? fiir hydrostatischen Innendruck p, = yz, py =0 und die 
erstmals von Milankovié? angegebene, mit 
Pp, =—ytcosg, Po =Ytsng, (1.5) 
d. h. nur durch ihr Eigengewicht belastete Kuppel gleicher Festigkeit. 


Abb. 1. Schalenelement. Abb. 2. FlieBfigur (oder Bruchfigur). 


Der Ansatz (1.4) hat den Vorzug, an keine FlieBbedingung (bzw. Bruchhypothese) gebunden 
zu sein; er beschraénkt aber andererseits die konstruktiven Méglichkeiten. Im Falle der Kuppel 
gleicher Festigkeit legen die Beziehungen (1.1) bis (1.5) sowohl die Mittelflache wie den Verlauf der 
Wandstarke vollstandig fest. Es unterliegt keinem Zweifel, daB die so erhaltene Kuppelform weder 
vom dkonomischen noch vom Asthetischen Standpunkt zu iiberbieten ist; insbesondere diirfte die 
Kuppel von Milankovié bei gegebener Schalenstarke t, im Scheitel mit einem Minimum an Material 
auskommen, um eine gegebene (und nicht zu grofe) kreisférmige Offnung zu tiberspannen. Fiir 
Falle, die (etwa aus Griinden einfacherer Herstellung) eine andere Mittelflache verlangen, bietet der 
Ansatz (1.4) aber keine Ausweichméglichkeit. 

In der Spannungsebene (¢,, dy) wird die FlieBbedingung von Tresca (bzw. die Bruchhypothese 
konstanter Schubspannung) durch das Sechseck von Abb. 2 dargestellt, die FlieBbedingung von 

1 W. Fliigge, Stati 

ae rie eee Ee der Schalen, 2. Aufl., S. 26. Berlin 1957. 


°M. Milankovié, Arbeiten der Jugoslaw. Akad. der Wissensch., Agram 175 (1908) S. 140. Man vergleiche 


auch W. Fliigge, a. a. O., S. 36, sowie de nf: i J i 
Foeuakoit: Bee 30 (1939) 8. 366 n zusammenfassenden Artikel von K. Federhofer, Uber Schalen gleicher 
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v. Mises (baw. die, Bruchhypothese konstanter Gestaltainderungsarbeit) durch die umgeschriebene 
Ellipse. Dabei bedeutet og die FlieBgrenze (bzw. Bruchgrenze) bei reinem Zug oder Druck. Die 
Bedingung (1.4) verlangt, daB der Spannungszustand in der ganzen Schale durch den Punkt A 
oder C in der Spannungsebene dargestellt werde. Das ist eine unndtige Einschrankung, denn je 
nachdem man der Betrachtung die Bedingung von Tresca oder diejenige von v. Mises zugrunde legt, 
kommt fiir die verschiedenen Schalenelemente als Bild in der Spannungsebene die Peripherie des 
ganzen Sechsecks bzw. der ganzen Ellipse in Frage. 

Es ist klar, dafS man bei dieser letzten Betrachtungsweise, die meines Wissens bisher nur von 
Biezeno’ auf Kesselbéden angewandt worden ist, iiber einen Freiheitsgrad mehr verfiigt. So kann 
man etwa wie Biezeno die Wandstarke vorschreiben und die zugehérige Mittelflache bestimmen. 
Oder man kann umgekehrt (innerhalb gewisser Grenzen) die Mittelflache vorschreiben und den zu- 
gehérigen Verlauf der Wandstarke berechnen. 

Wenn (1.4) nicht mehr gilt, besteht die Verformung freilich nicht 
mehr in einer Ahnlichkeitstransformation. Es ist daher infolge der 
Deformation ein sekundarer Biegespannungszustand zu erwarten, 
der aber bei diinnen Schalen bekanntlich vernachlassigt werden darf. 
Auch an den Randern hat man ja in jedem Fall durch eine geeignete 
Lagerung allfallige Biegespannungen zu vermeiden. 

Im Folgenden soll als Beispiel die Kuppel unter ihrem Eigen- 
gewicht betrachtet und gezeigt werden, da sich eine spharische 
Kuppel im Sinne der Bedingung von Tresca (sogar auf zwei Arten) als Kuppel gleicher Festigkeit 
ausbilden ]aBt. 


Abb. 3. Spharische Kuppel. 


2. Die spharische Kuppel; Lésung 1. Bei einer Kugelschale (Abb. 3) sind die Radien 
fT. r=asing. (2.1) 
Setzt man (2.1) zusammen mit (1.3) und (1.5) in (1.1) ein, so erhalt man fiir die spharische Kuppel 
unter ihrem Higengewicht die Beziehungen 
5 (to, sin g) —t oy 005 9 + ve sin go 0 Op + 0g + ay cosy —0. (2.2) 
Im Scheitel p = 0 gilt o, =o, < 0; man hat also bei vélliger Ausnutzung des Materials 


Oo ad) OGG — 0 = OG, (2.3) 
Der Bildpunkt P des Scheitels in der Spannungsebene ist A, und aus der zweiten Beziehung (2.2) 
sowie (2.3) ergibt sich der Radius der Mittelflache zu 


a= 2%. (2.4) 


Geht man einem Meridan entlang nach unten, dann kann nach den Feststellungen von Abschn. 1 
der Bildpunkt P nicht in A bleiben. Er muf langs einer der Seiten J oder 2 des Sechsecks abwandern, 
und diese beiden Méglichkeiten stellen die zwei bereits erwahnten Falle dar. 

Auf der Seite 1 des Sechsecks gilt 


J, = — 06 (2.5) 
und wegen der zweiten Beziehung (2.2) sowie (2.4) 
Oy =O (1—2 cos). (2.6) 


Setzt man dies, (2.5) und (2.4) in die erste Beziehung (2.2) ein, so erhalt man fir t die Differential- 
gleichung 


4 [t (1—2 cos g) sin g] + ¢ (008 p + 2 sing) = 0 (2.7) 
mit dem allgemeinen Integral 
t =c(1 + cos g)—7/3 (1 —2 cos g)— 78. (2.8) 
Mit der Randbedingung 
ti(p =0) = (2.9) 
wird hieraus 
t = — 27/3 4, (1 + cos g)— 23 (1 — 2 cos g)— 4, (2.10) 


1C, B. Biezeno, Der Ingenieur 39 (1922) S. 781. 
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und damit ist der Verlauf der Schalenstarke gegeben, sobald ihr Kleinstwert t) im Scheitel gewahlt 


ist. 


Auf der Seite 1 des Sechsecks (Abb. 2) gilt neben (2.5) 
—oS0o,59. (2.11) 


Setzt man hier (2.6) ein, so kommt nach kurzer Zwischenrechnung 


1 
Bote =a 2.12 
5 -cosgy Sl. ( ) 


8 


0° a 2° go 4° ag? Glenn 0) amines mame 
Abb. 4. Verlauf der Schalenstirke. 


Der Punkt P wandert also von A bis B, wahrend von 0 bis 60° wachst, und auf dieses Intervall 
ist somit die Lésung (2.10) beschrankt. In Wirklichkeit gilt sie sogar nicht ganz so weit, wie die 
Tabelle 1 bzw. die Kurve I in Abb. 4 zeigen, welche das Verhaltnis t/t, als Funktion von ¢ darstellen. 
Fiir p— 60° geht t/t) > co; die Lésung 1 kommt demnach nur fir verhaltnismaBig flache Kugel- 
schalen (gy < 60°) in Frage. 


Tabelle 1. Verlauf der Schalenstarke fiir Lésung 1 


¢ t/t | ~° t/ty 
5 1,019 35 3,037 
10 1,080 | 40 4,736 
15 1,193 45 8,649 
20 1,378 | 50 PINS 
25 1,676 | 55 102.6 
30 2,168 60 0° 


3. Die sphirische Kuppel; Lésung 2. Auf Seite 2 des Sechsecks (Abb. 2) gilt 
Oy = — Dp (3.1) 
und wegen der zweiten Beziehung (2.2) sowie (2.4) 
Og = 0) (1—2 cosq). (3.2) 


Setzt man dies, (3.1) und (2.4) in die erste Bezichung (2.2) ein, so erhalt man fir t die Differential- 
gleichung 


ee 
a (t sin p) — t (2 — cos py) = 0 (3.3) 
mit dem allgemeinen Integral 


t=c(1+cosqg)?. (3.4) 
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Mit der Randbedingung (2.9) wird hieraus 
t=4t(1+cosg)?*=t, cos 4 ; (3.5) 


und damit ist ein zweiter Verlauf von ¢ auf t) zuriickgefiihrt. 
Auf der Seite 2 des Sechsecks (Abb. 2) gilt neben (3.1) auch 
=Go- 6, = 0, (3.6) 


und wenn hier (3.2) eingesetzt wird, dann kommt man wieder auf (2.12). Der Punkt P wandert also 
von A nach D, wahrend ¢ von 0 bis 60° anwachst, und auf dieses Intervall bleibt die Lésung (3.5) 
beschrankt. Tab. 2 und die Kurve 2 in Abb. 4 geben den Verlauf von t/t, als Funktion von : 


Tabelle 2. Verlauf der Schalenstarke fiir Lésung 2 


e t/to | gy t/t, 

5 1,004 35 1,209 
10 1,015 40 1,282 
15 1,035 | 45 Mode 
20 1,063 50 1,482 
25 1,101 a) 1,615 
30 1,149 | 60 1,778 


Da fiir p = 60° erst t/t) = 16/9 wird, ist die Lésung 2 durch einen viel schwacheren Anstieg der 
Schalenstarke charakterisiert als die Lésung I. Sie ist somit wesentlich ékonomischer. Zudem laBt 
sie sich tiber den Winkel 60° hinaus fortsetzen. Wenn namlich P in D auf die Seite 3 des Sechsecks 
(Abb. 2) iibertritt, gilt 

O59 =O) + Oy (3.7) 


und wegen der zweiten Beziehung (2.2) sowie (2.4) 
Op =— 9p 5 ++ cos ®) : (3.8) 


Setzt man dies, (3.7) und (2.4) in die erste Beziehung (2.2) ein, so erhalt man fir t die Differential- 
gleichung 
alt (1 + 2 cos p) sing] —t (2 — cos y + 2 sin? y) = 0 (3.9) 
mit dem allgemeinen Integral 
t = c (1 + cos g) (1 — cos y)~ 43 (1 + 2 cos g)— 3, (3.10) 
Der Anschlu8 an (3.5) wird durch die Randbedingung 


t(p = 60°) =, (3.11) 
hergestellt, und diese liefert 
t= PIF ty (1 + cos g) (1 — cos g)— 49 (1 + 2 cos g)— MP. (3.12) 
Auf der Seite 3 des Sechsecks (Abb. 2) gilt 
Jer pay (3.13) 
Setzt man hier (3.8) ein, so folgt 
—+Scsp ss. (3.14) 


Der Punkt P wandert demnach von D nach E, wahrend @ von 60° bis 120° wachst, und tiber diesen 
letzten Wert hinaus laBt sich die Lésung (3.12) nicht fortsetzen. In der Tat wachst t/t), wie die 
Tabelle 3 bzw. die Kurve 2 in Abb. 4 zeigt, im Giiltigkeitsbereich von (3.12) sehr stark an und geht 


fiir py — 120° gegen ov. 
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Tabelle 3. Verlauf der Schalenstérke bei der Fortsetzung 
von Lésung 2 


gv t/t ~ t/t 

60 1,778 95 21,40 
65 2,199 100 43,81 
70 2,727 | 105 118,8 
US 3,600 110 486,6 
80 5,008 amis) wPO195 
85 7,428 120 oo 

90 11,95 


Fir p = 90° folgt aus (3.12) t/t) = 11,95. Bei einer Halbkugel wachst also die Schalenstarke 
vom Scheitel bis zam Rand rund auf das 12 eae: und zwar entfallt der sa Teil des Anstiegs 
auf die letzten 30°. 

Der Breitenkreis ¢ = 60° entspricht dem Ubergang des Punktes P von der Seite 2 im Sechseck 
(Abb. 2) auf die Seite 3. Er stellt die Bruchfuge dar, an welcher oy das Vorzeichen wechselt und zur 
Zugspannung wird. Die Funktion ¢(y) hat hier einen Sprung in der Ableitung. 


4, Diskussion. Die Lésungen J und 2 sind an die Voraussetzung gebunden, da das Material der 
FlieBbedingung von Tresca (bzw. der Bruchhypothese konstanter Schubspannung) geniige. 

Der Radius a der spharischen Kuppel ist fiir beide Lésungen (J und 2) derselbe. Er ist nach (2.4) 
durch das Material bestimmt, sofern man dieses voll ausniitzen will. Fir Beton mit y = 2,4 g/cm® 
und oy = 150 kg/cm? wird z. B. a = 1250 m, fiir einen Stahl mit y = 7,85 g/cm? und 0, = 1500 kg/ 
cm? sogar a = 3820 m. 

Man kann natiirlich fiir oy nur einen Bruchteil der zulassigen Spannung ansetzen und erhalt so 
eine Kuppel von kleinerem Radius, die noch immer iberall die gleiche Festigkeit besitzt. Ferner 
diirfte es nicht schwierig sein, die Ergebnisse fiir den Fall zu verallgemeinern, daB neben dem Eigen- 
gewicht noch eine konstante Last je Flacheneinheit der Schale (Verputz, Isolation) oder ihrer 
Horizontalprojektion (Schneelast) vorhanden ist. 

Nach (2.4) stimmt der Radius a mit dem Kriimmungsradius tiberein, welchen die Kuppel von 
Milankovié im Scheitel aufweist. Die Mittelflache der hier diskutierten Kuppel ist daher die oskulie- 
rende Kugel im Scheitel derjenigen von Milankovié. Bei dieser nimmt mit wachsendem Abstand 
vom Scheitel der Radius r, zu, wahrend bei der Kugel r, = a bleibt.. Mit der Kuppel von Milankovié 
lassen sich also gréBere Flachen tiberspannen als mit der spharischen. Andererseits geht t/t, bei 
der Kuppel von Milankovié fiir py — 90° nach oo, wahrend bei der spharischen Kuppel (Lésung 2) 
der Winkel y = 90° leicht iiberschritten werden kann. 

Die hier behandelten Lésungen stellen nur zwei von unendlich vielen Méglichkeiten dar, eine 
Kuppel mit gleicher Festigkeit auszustatten. Man beachte in diesem Zusammenhang neben der 
Freiheit in der Wahl der Mittelflache die Méglichkeit anderer FlieBfiguren (oder Bruchfiguren), 
wie etwa der v. Misesschen Ellipse in Abb. 2. Hinsichtlich der Form der Mittelflache bestehen in- 
dessen gewisse Schranken. So laBt sich beispielsweise das Problem fiir eine Kegelschale ebenso- 
wenig lésen wie fiir eine oben offene (d.h. mit einer kreisformigen Offnung um den Scheitel versehene) 
spharische Schale. 


(Eingegangen am 24, Dezember 1957.) 
Anschrift des Verfassers: Professor Dr. Hans Ziegler, Riischlikon bei Ziirich, Weiherweg 6 
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: = : % Dreizehnte neubearbeitete Auflage. Von Dr.-Ing. Wilhelm Rai®, o. Professor an der Tech- 

eft oe nischen Universitat Berlin, Direktor des Institutes fiir Heizung und Liiftung. Unter Mit- 
Ke - arbeit von Dr. phil. habil. F. Bradtke }. Mit 524 Abbildungen, 43 Zahlen- und 2 Bildtafeln 
ae sowie den Axbeitsblattern 1—15, XII, 568 Seiten 4°, ee Ganzleinen DM 79,50 
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Adie dharanr: I. Teil: Es warcmaen's zur Heizung und ie ate Einzelheizung: Ofen fiir feste 
Brennstoffe. Der Schornstein, Gasheizéfen. Elektrische Raumheizung. — Zentralheizung: Allgemei- 
nes. Bauelemente der Warmwasser- und Dampfheizungen. Warmwasserheizungen. Niederdruckdampf- 
_ heizungen. Heizzentralen. Sonderbauarten der Dampf- und Wasserheizung. Luftheizung. Gesichtspunkte 
zur Wahl des Warmetragers und der Heizungsart. — Zentrale Warmwasserbereitung: Warmwasser- 
bereitung mit Gebrauchswasserspeicherung. Warmwasserbereitung mit HeiBwasserspeicherung. Heiz- 
-flachen- und Speicherbemessung. Verbrauch und Wirtschaftlichkeit. Steinbildung und Korrosion. — 
Fernheizung: Allgemeines. Fernleitungen. Dampffernheizung. Hei8wasserfernheizung. Warmwasser- 
_ fernheizung. Stadtheizung. — Heizkraftanlagen. Warmepumpe. — Liiftungs- und Klimatechnik: 
% Allgemeines. Freie Liiftung. Liiftungsanlagen. Klimaanlagen. — II. Teil: Berechnung von Heiz- und 
Liiftungsanlagen. Hygienische und warmephysiologische Grundlagen: Einleitung. Warme- 
regelung des menschlichen Kérpers. Durch die Haut vermittelte Einfliisse von Umgebungsluft auf den 
‘menschlichen ‘KG6rper. Behaglichkeitsmafstabe. Durch die Atmung vermittelte Einfliisse der Umgebungs- 
luft fiir den menschlichen Kérper. — Meteorologisch-klimatische Grundlagen: Hinleitung. Die 
Temperatur der AufSenluft. Die Feuchtigkeit der Auenluft. Der Wind. Die Sonnenstrahlung. Physika- 
i lische Grundlagen fiir das Rechnen mit feuchter Luft.— Warmeiibertragung: Allgemeines. Warme- 
= perce --Jeitung. Warmeiibergang. Warmedurchgang. ‘Warmeiibertragung durch Strahlung. — W armetechnische 
_ _Berechnungen: Der Warmebedarf. Kesselheizflachen. Raumheizkorper. Decken- und Fufbodenheizung. 
_ Warmeaustauscher. Rohrisolierungen. — Strémungsfragen: Die Gesetze fiir die Strémung in Lei- 
_ tungen. Der Druckverlust in Leitungsnetzen. Ventilatoren und Kreiselpumpen. — Berechnung von 
_ Rohrnetzen: Allgemeines. Berechnung von Fernleitungen. Berechnung der Strangnetze von Warm- 
wasserheizungen. Berechnung der Strangnetze von Niederdruckdampfheizungen. Die Berechnung von 
Hochdruck- und Unterdruckdampfheizungen. — Betrieb von Heizanlagen: Begriffe zur Warmebilanz. 
Der Heizwarmebedarf und seine Veranderlichkeit. Warmebedarf und Warmeverbrauch. Heizbetrieb. — 
Berechnung von Liftungs- und Klimaanlagen: Berechnung der Luftverteilungsleitungen. Be- 
rechnung einer Liftungsanlage. Klimatechnische Berechnungen. — Il. Teil: Zahlen- und Bildtafeln. 
_ Zahlentafeln: ‘Stoffwerte und warmetechnische Tabellen. Warmebedarfsrechnung. Berechnung von 
; Heizflachen -und Isolierungen. Rohrnetzberechnung. Allgemeines. Bildtafeln, — Regeln, Richtlinien, 
Normen. -- Namenverzeichnis. — Sachverzeichnis. 


Die asfamercicke Bautatigkeit der letzten : zehn Jahre hat die technische Entwicklung auf dem Fachgebiet 
Heizung und Liiftung nachhaltig geférdert. Neue Bauformen von Kesseln und Heizkérpern sowie die zu- 
- nehmende Verwendung von Deckenheizungen haben das Gesicht der Heizanlagen teilweise geindert. Bei 
der Warmeversorgung gewinnt fiir dichtbesiedelte Stadtgebiete die Fernheizung und fiir Einzelanlagen die 
Olfeuerung immer mehr an Bedeutung. Ahnlich einschneidende Fortschritte sind i im Bau von Luftungs- 
und Klimaanlagen zu erkennen. — 
Damit ist eine vollige Neubearbeitung ser Lehrbuches notwendig geworden, die sich nicht nur auf die 
beschreibenden Abschnitte, sondern auch auf den gesamten Berechnungsteil erstreckt. Eingehender als 
_ geither werden in der 13. Auflage die Gesetze der Warmeiibertragung behandelt. Zur Vereinfachung von 
-Warmeaustauschrechnungen werden die Losungen schwieriger oder umfangreicher Formeln weitgehend in 
- Diagrammform gebracht. Dem Abschnitt iiber die Bestimmung des Warmebedarfs von Gebauden liegt die 
_ jetzt fertiggestellte Neufassung der DIN 4701 zugrunde. Experimentelle und theoretische Untersuchungen 
aus den letzten Jahren erméglichen vielfach eine genauere Ermittlung der Wiarmeleistung verschiedener 
Heizkérperbauarten bzw. Heizsysteme. U. a. werden fiir die Berechnung von Deckenheizungen neue Ver- 
fahren angegeben, deren Anwendung an Beispielrechnungen erlautert wird. 


re Die Rietschelschen Hilfstafeln zur Bemessung der Leitungsnetze von Heizanlagen sind nach dem derzeitigen 

Stand unseres Wissens tiber die Rohrreibung neu berechnet und durch weitere, zum Teil graphische Tafeln 

erganzt worden. Fiir die Berechnung von Fernleitungen und von Schwerkraftheizungen mit oberer Vertei- 

- lung werden gegeniiber den fritheren Auflagen vereinfachte Verfahren angegeben. Ausfihrlich wird auf 

die Berechnung von Klimaanlagen eingegangen. Ein Sonderabschnitt ist den Fragen des Betriebes und der 
Wirtschaftlichkeit von Heizanlagen gewidmet. 
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Materialpriifung 
mit Rénteenstrahlen 


unter besonderer Beriicksichtigung der Réntgenmetallkunde 


Von 
Dr. phil. Dr. med. h. c. Richard Glocker 
Professor fiir Réntgentechnik an der Technischen Hochschule Stuttgart 


Vierte erweiterte Auflage. Mit 455 Abbildungen. VII, 530 Seiten Gr.-8°. 1958, Canzleinen DM 61,50 


Inhaltstibersicht: Einteilung: Die Natur der Réntgenstrahlen und die verschiedenen Verfahren der 
Werkstoffpriifung mittels Réntgenstrahlen. — Erzeugung der Réntgenstrahlen: Allgemeines itiber Rontgen- 
réhren. Ausfiihrungsformen von Réntgenrdbren und Glihventilréhren. Réntgenapparate. Strahlenschutz. 
— Eigenschaften der Réntgenstrahlen: Absorption und Sekundarstrahlung. Beugung und Brechung. 
Ionisation und photographische Wirkung. — Grobstrukturuntersuchung: Allgemeine Grundlagen der Grob- 
strukturuntersuchung. |Praktische Anwendung der Grobstrukturuntersuchung. — Spektralanalyse: Rént- 
genspektroskopische Apparate. GesetzmaBigkeiten des Rontgenspektrums. Qualitative Réntgenspektral- 
analyse. Quantitative Réntgenspektralanalyse. — Feinstrukturuntersuchung: Uberblick iiber die verschie- 
denen Verfahren der Feinstrukturuntersuchung und ihre Anwendungsgebiete. Kristallographische Grund- 
lagen. I. Kristallographische Grundlagen. II]. DEBYE-SCHERRER-Aufnahmen (Pulverdiagramme). 
Diffraktometer (Zahlrohrgoniometer). LAUE-Aufnahmen und Drebkristallverfahren. Photographische 
Réntgengoniometerverfahren. Intensitatsgesetze der Réntgeninterferenzen. Uberblick iiber den Gang 
einer Strukturbestimmung. Beschreibung von Kristallstrukturen anorganischer und organischer chemischer 
Stoffe und Grundziige der Kristallchemie. Struktur von Legierungen. Verbreiterung der Réntgeninter- 
ferenzen und Bestimmung der Kristallgré8e. Messung von elastischen Spannungen. Kristalltexturen. Nicht- 
kristalline Stoffe und Flissigkeiten. Kleinwinkelstreuung. — Mathematischer Anhang. — Schrifttum — 


Namenverzeichnis — Sachverzeichnis. 


Die Anwendung der Réntgenverfahren erfordert gewisse Vorkenntnisse. In zwei einleitenden Abschnitten 
werden die Réntgenanlagen und die Eigenschaften der Réntgenstrahlen einschlieBlich des Strahlenschutzes 
besprochen. Bei der Grobstrukturuntersuchung zum Nachweis von Fehlstellen z. B. bei SchweiBungen ge- 
winnt die Anwendung der Strahlungen von radioaktiven Isotopen steigende Bedeutung. Die Réntgenspek- 
tralanalyse ist seit der Einfiihrung des Zahlrohres eine Routinemethode der Industrie z. B. zur laufeaden 
Bestimmung des Molybdangehaltes von Stahlproben. Bei der Feinstrukturuntersuchung ist die Metall- 
kunde besonders beriicksichtigt; technisch wichtige Vorgange, wie z. B. die Aushartung von Legierungen, 
verdanken ihre endgiiltige Aufklarung der Réntgenstrukturuntersuchung. Die rontgenographische Messung 
von Spannungen erméglicht nicht nur die Ermittlung der Eigenspannungen, die als Folge von Bearbeitungs- 
vorgangen oder Warmebehandlungen in einem Werkstiick entstehen, sondern verspricht auch ein wichtiges 
Hilfsmittel zum Studium der Zerriittung bei Wechselbeanspruchung zu werden. Der Bestimmung der Atom- 
anordnung in amorphen festen Kérpern und in Metallschmelzen ist ein besonderer Abschnitt gewidmet. 
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